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A certain man had the power of grunting exactly like a pig, so much so that whenever he 
grunted where pigs were grazing, they would all turn to see if any new member had come 
into their fold. Slowly he became popular and started presenting a show of his art for a 
price. Once when he was showing his skills in a village pandal, a sage with his disciples 
chanced to pass through the village and was surprised to hear about the show. The sage 
immediately got another pandal erected and advertised that there was going to be a bigger 
grunting show and that too without any tickets. Naturally a big crowd assembled to watch 
the feat. But, to their, utter disappointment, the sage was just making an actual pig grunt 
loudly. The crowd left the pandal in sheer disgust claiming that they heard this everyday 
and that the sage was merely wasting their time. No amount advertising could get the crowd 
back, whereas, the show of the man imitating the pig’s grunt was still pulling huge crowds. 
Now the sage addressed his disciples, “Here is a splendid lesson for us. Men seldom care 
for reality, but always go in for imitation. That is why this world exists, which is a mere 
imitation, a reflection in the distorting mirror of Maya of the great Atman. No external 
help is needed to see the self but very few want it and even if you eagerly advertise it, none 
will go to you except those who love Truth for Truth’s sake. Reflect on this”. 
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The ’’head” of amphiphilic molecules (also called surfactants) consist of a polar or ionic 
group. The ’’tail” of many amphiphiles consist of a single hydrocarbon chain whereas that 
of some other amphiphiles, e.g., phospholipids, are made of two hydrocarbon chains both of 
which are connected to the same head. Yet another class of surfactants, called gemini surfac- 
tants, consist of two single-chain amphiphiles whose heads are connected by a ’’spacer” chain 
and, hence, these ’’double-headed” surfactants are sometimes also referred to as ’’dimeric sur- 
factants”. When put into an aqueous medium, the ’’heads” of the amphiphiles like to get 
immersed in water and, hence, called ’’hydrophilic” while the tails tend to minimize contact 
with water and, hence, called ’’hydrophobic” . The spacer in gemini surfactants is usually 
hydrophobic but gemini surfactants with hydrophilic spacers have also been synthesized. 
In aqueous medium amphiphiles optimize their contact with water forming various types 
of self-assemblies (i.e., supramolecular aggregates), e.g., monolayer and bilayer membranes, 


micelles, vesicles, etc. 



The aim of the thesis is to develop, generalize and/or extend realistic, but simple micro- 
)pic models of self assembling amphiphilic systems and to study their bulk and interfacial 
)perties. We check the reliability of the models by comparing the results of computer 
lulation with known experimental results and, then, predict new phenomena by carrying 
fc novel computer experiments. Since Monte Carlo (MC) simulation is more efficient than 
)lecular Dynamics (MD) method for the equilibration of such systems, we compute most 
the quantities of interest by MC simulation. 

In chapter 1 we give a brief overview of the physical properties of the amphiphilic 
items, the nature of the seif assemblies formed by the amphiphiles, and also describe the 
joretical and computational techniques to model these systems. Chapters 2,3 and 4 deal 
:h almost exclusively the interfacial properties, chapter 5 deals with both the interface 
d the bulk while chapters 6 and 7 deal with self-assemblies in bulk and chapter 8 deals 
:h self-assemblies in porous media. Our conclusions are summarized in chapter 9. 

It is known from continuum models of membranes/interfaces that when bending rigidity 
ne governs the out-of-plane thermal fluctuations, amphiphilic membranes can crumple 
3 ve a critical temperature. In the chapter 2 of the thesis we show that the oil-water 
erface is rough, but not crumpled, in the droplet phase of the Widom model which is the 
iplest microscopic lattice model of ternary microemulsions. 

In chapter 3 we construct a class of discrete lattice models of membranes/interfaces 
ere we can make a clear distinction between the interfacial tension and bending rigidity 
ich can be switched OS'’ independently of each other. More speciflcally, in chapter 3, we 
pose a dH- 1-dimensional restricted discrete solid-on-solid (RDSOS) model as well as un- 
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restricted discrete solid-on-solid (UDSOS) model for d-dimensional amphiphilic membranes 
where the bending stiffness of the membrane is incorporated following the same prescrip- 
tion as in the Widom model. We also point out some analogies between the relaxation to 
rough equilibrium states of these SOS models and nonequilibrium growth of some kinetically 
roughened model surfaces. 

In all the above models bending energy of the amphiphilic monolayer membranes is 
introduced as a parameter. In the remaining chapters we use the Larson prescription for 
modelling amphiphilic systems, as it explicitly takes into account the conformation^ degrees 
of freedom of the amphiphilic chain molecules. Jan, Stauffer and collaborators reformulated 
the Larson model in terms of Ising-like variables. We have used Jan and Stauffer’s prescripi- 
tion for our models in chapters 4 r 8. 

In chapter 4 we study the roughening of the oil-water interface in the presence of single- 
chain surfactants and investigate the physical origin of the lowering of the oil-water interfacial 
tension. 

After studying the properties of the amphiphilic membrane at the oil-water interface 
in chapters 2,3 and 4, we investigate the air-water interface in the presence of single-chain 
amphiphiles in chapter 5. In this chapter we introduce a microscopic lattice model of a 
binary mixture of amphiphilic molecules, of two different lengths, in a system where water 
is separated from the air above it by a sharp well defined interface. We demonstrate an 
entropy-driven phase segregation in a direction perpendicular to the air-water interface when 
the initial total surface density of the amphiphiles is sufficiently high. We also investigate 
(a) the conformations of the amphiphiles, (b) the distribution of the sizes of the clusters 



of monomers belonging to the long amphiphiles as well as of those belonging to the short 
amphiphiles in planes parallel to the interface, (c) the effects of varying the lengths, total 
concentration and the ratio of the numbers of the two types of the amphiphilic molecules, 
:he strength of the inter-monomer interactions and temperature. 

In the chapters 4 and 5 only single-chain amphiphiles have been considered although the 
general conclusions drawn from our MC studies are valid also for double-chain amphiphiles. 
n chapter 6 we compare and contrast the self-assemblies of single-chain and double-chain 
imphiphiles, respectively, in water. In this chapter we study the self-assemblies of binary 
nixtures of single-chain and double-chain amphiphiles and compare these with those formed 
ly the corresponding single species of amphiphiles. 

In chapter 7 we propose a microscopic lattice model of gemini surfactants. We show 
tat these model gemini surfactants form long, thread-like and entangled micellar aggregates 
\?en at low concentrations where single-chain model surfactants form spherical micelles. We 
jmpare the morphologies of the micellar aggregates formed in our MC simulations with 
lose observed in laboratory experiments and in earlier molecular dynamics simulations. 
h also study the variation of the critical micellar concentration of these model gemini 
irfactants with the variation of the (a) length of the hydrocarbon spacer connecting the 
ro hydrophihc heads, (b) length of the hydrophobic tail and (c) the bending rigidity of the 
'drocarbon chains forming the spacer and the tail. Some of the trends of variation are 
unter-intuitive but are in excellent agreement with the available experimental results. 

In chapter 8 we propose a scheme for MC simulation of ternary micro emulsions in porous 
jdia (such as vycor glass). Finally, in chapter 9 we summarize our conclusions and outline 



directions of future work. 
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Chapter 1 


Introduction 


1.1 The Aim 

Soap molecules are common examples of amphiphilic molecules; these not only find wide 
ranging applications in detergent and pharmaceutical industries, food technology, petroleum re- 
covery, etc. but are also one of the most important constituents of ceils in living systems. There- 
fore, physics, chemistry, biology and technology meet at the frontier area of interdisciplinary 
research on association colloids formed by surfactants [1-5]. The amphiphiles give rise to many 
interesting properties of complex (multi-component) fluids containing water and amphiphiles. For 
example, we know that at ordinary room temperature oil and water do not mix because of the high 
oil-water interfacial tension. However, the effective oil-water interfacial tension can be reduced 
to almost vanishingly small value by adding soap to the otherwise immiscible binary mixture of 
oil and water. That is why amphilhilic molecules are often referred to as surfactants (surface 
active agents). 
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ipite of decades of research efforts, a clear microscopic understanding of the behaviour of 
: fluids containing amphiphilic molecules has not yet emerged. In most practical applica- 
ie amphiphilic systems are so complex that it is difficult to calculate any physical property 
system analytically starting from a microscopic (molecular) description which includes all 
sails of these systems without making drastic approximations. Computer simulation of 
:opic models can provide valuable insights into the bulk and interfacial properties of the 
)hilic systems, such as, the physics of the formation of supramolecular aggregates and their 
ologies as well as the influence of these on the macroscopic physical properties. In order 
n insight into the physical origin of some of these bulk and interfacial properties, in this 
we develop, generalize and/or extend realistic, but simple microscopic models of complex 
-component) fluids containing amphiphilic molecules and water. We check the reliability of 
odels by comparing the results of computer simulation with known experimental results and, 
predict new phenomena by carrying out novel computer experiments. Since Monte Carlo 
simulation is more efficient than Molecular Dynamics (MD) method for the equilibration 
:h systems, we compute most of the quantities of interest by MC simulation. 

The Amphiphilic Molecules and Their Supramolec- 
ular Aggregates 

The "head" part of amphiphilic molecules consist of a polar or ionic group. The “tail” of 
/ amphiphiles consist of a single hydrocarbon chain whereas that of some other amphiphiles, 
phospholipids, are made of two hydrocarbon chains both of which are connected to the 
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Figure 1.1: Different types of amphiphiles. 


same head [6, 7, 9, 10]. Yet another class of amphiphiles, called gemini amphiphiles [11-14], 
consist of two single-chain amphiphiles whose heads are connected by a "spacer” chain and, 
hence, these "double-headed" amphiphiles are sometimes also referred to as “dimeric surfactants” 
[15, 16]. Some common example of the different amphiphiles are as follows: Single chain: 
Sodiumdodecylsulphate {Na'^SO^ ■ {CH 2 )nCH 3 }, 

Double chain: Dihexadecyl dimethylammonium bromide {Cie-ffas • N'^{CHz) 2 Br~ ■ Cig-ffas}, 
and Gemini surfactants: {Br", m — - NMe^ — {CH2)n — N'^Me2 — m — C 16 H 33 , Br~} 

where n is the length of the spacer. Different types of model amphiphiles are shown in figure 1.1 

When put into an aqueous medium, the “heads” of the amphiphiles like to get immersed 
in water and, hence, called “hydrophilic” while the tails tend to minimize contact with water 
and, hence, called “hydrophobic" [10]. The spacer in gemini amphiphiles is usually hydrophobic. 
However gemini amphiphiles with hydrophilic spacers have also been synthesized [17]. Therefore, 
in aqueous medium, amphiphiles form a variety of "self assemblies" (supramolecular aggregates). 
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.gure 1.2: Different structures formed by amphiphiles in aqueous medium (after[6]). 

; monolayer and bilayer membranes, micelles, inverted-micelles, vesicles, etc. (figure 1.2). 
^ varying the chemical structure, size, shape and flexibility of the amphiphiles, the relative 
tration of water, surfactant (and other components in the fluid), temperature and external 
•e the system moves from one region of the phase diagram to another. A brief overview of 
perimentally observed phase diagram and the associated physical properties of binary and 
/ fluid mixtures, which contain amphiphiles, is given in the next section. Many attempts 
)een made to develop simple models to describe the generic behaviour of these systems, 
tion 1.4 a short summary of the theoretical and computational techniques used to model 
systems is presented. 


Phase Diagrams and Physical Properties 

/arious experimental techniques have been used to study these systems; these include freeze- 
ire electron micrograph, X-ray, light scattering and neutron scattering experiments, nuclear 
etic resonance, transmission electron microscopy, measurement of conductivity and rheo- 
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logical measurements, e.g. measurement of viscosity. 

1.3.1 Interfacial properties 

One of the most spectacular effects of the amphiphiles on a mixture of oil and water is that a 
small amount of amphiphiles can reduce the effective interfacial tension between water and oil by 
several orders of magnitude. A particularly clear example is described below. Experiments were 
performed on a mixture of ionic amphiphile (Sodiumdodecylsulphate), water, salt and organic 
solvents [21]. The system was observed as the concentration of salt was increased and the 
system moved from (1) a mixture of organic phase and oil/water microemulsion to (2) a mixture 
of organic phase, a middle phase microemulsion and an aqueous phase and finally to (3) a mixture 
of water/oil microemulsion and an aqueous phase. It was observed that at the critical end-point 
at which the middle phase and the water-rich become identical, the interfacial tension between 
them goes to zero. The interfacial tension between the middle and the oil-rich phase also vanishes 
at another critical endpoint, where the middle and oil-rich phase are identical. Since the oil-water 
interfacial tension (which is only defined in the region of three phase coexistence) must always be 
less than the sum of the interfacial tensions between oil-rich and middle phases and the aqueous 
and middle phases, the oil-water tension was dramatically reduced, by 3 orders of magnitude in 
this case. 


1.3.2 The Phase Diagram 


• Phases of ternary systems: 
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re 1.3: Schematic diagram illustrating the phase behaviour of real surfactant-oil-water sys- 
. (after[43]). 
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A schematic phase diagram of ternary system consisting of oil, water and surfactant has 
been shown in figure 1.3. Ternary mixtures of oil, water and surfactants, where the concentration 
of amphiphiles is low, are called tnicroemulsions [9]. Suppose the volume fraction occupied 
by oil is very small compared to that occupied by water. Then at the microscopic level, droplets 
of oil are dispersed in water [10, 18] and the amphiphiles form a monolayer at the oil-water 
interface. The amphiphiles arrange themselves, at the oil-water interface in such a manner that 
the hydrophilic heads are in contact with water and the hydrophobic tails are in oil. Conversely, 
at low amphiphile and low water concentration, one observes ‘droplets’ of water dispersed in oil, 
coated by a monolayer of amphiphiles with the heads pointing into the centre (water) and the tails 
pointing outwards (oil) [10, 18]. These systems are often called globular microemulsions 
[18]. Schematic diagrams of both structures are shown in figures 1.4(a) and (b). 

In contrast to emulsions, which are dispersions of particles of typical size lOA to lOOOOA 
in water globular microemulsions are dispersion of droplets of typical size lOOA [26] in water. 
Consequently, at the macroscopic level microemulsions are transparent, in contrast to turbid 
emulsions, as light does not get scattered by the droplets whose size is much smaller than the 
wavelength of light. 

When the concentrations of oil and water are comparable and much higher than that of the 
amphiphiles, an interconnected network of oil and water domains is formed with a monolayer of 
amphiphiles distributed at the oil-water interface [18]. This phase is usually called the bicon- 
tinuous micro emulsion [25]. A bicontinuous microemulsion is shown schematically in Fig. 
1.4(c); the monolayer at the interface between oil and water in such a phase is an example of 
amphiphilic monolayer membranes. 
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Figure 1.4: Different structures formed by amphiphiles in a mixture of oil-water-amphiphile 
system. Figure (d) and (e) after[2] and figure (c) after[19]. 
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As the concentration of amphiphiies is increased, liquid-crystalline phases with long- 
range periodic order are observed [28, 29]. Similar phases have been observed in liquids [27] 
consisting of small non-spherical molecules and hence the name. The phases observed are the 
lamellar phase (Fig. 1.4(d)), the hexagonal phase (Fig. 1.4(e)), the cubic phase and the 
inverted hexagonal phase [28]. In the lamellar phase, fiat or rippled layers of water and oil are 
separated by monolayers of amphiphiies. In the hexagonal phase, cylinders of oil domains, coated 
with amphiphiies, are so arranged that the centers of the circular cross sections of the cylinders 
form a (quasi-) two dimensional triangular lattice. The inverse hexagonal phase is obtained from 
the hexagonal phase by interchanging water and oil. Finally, in the cubic phase the centers of 
the spheres of amphiphile-coated oil droplets are arranged so as to form a three dimensional 
body-centered cubic lattice. 

It is worth mentioning that the ternary mixture exhibits a disordered fluid phase at sufficiently 
high temperature. 

•Phases of binary systems: 

Binary mixtures of water and amphiphiies also have a rich and varied phase behaviour. At 
very low amphiphile concentrations, hardly any aggregation of the amphiphiies takes place; iso- 
lated single amphiphiies are uniformly distributed in the mixture. Above a particular concentration, 
called the critical micelle concentration (CMC) [10], the amphiphiies form supramolecular 
aggregates which are called micelles. (Actually CMC is not a single concentration but a very 
narrow range of concentration and, therefore, it may be more appropriate to call it the character- 
istic micellar concentration). Just above the CMC one finds micelles dispersed in water. There 
is also a ‘transition’ from a micellar to a condensed phase at a particular temperature called the 
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idensation transition temperature. Below this temperature, the amphiphiles are in a 
densed phase where they form one large cluster and above this temperature, they from stable 
:elles of varying sizes [78], 

Ordered liquid crystalline phases have also been observed [30, 31]. With the increase of am- 
phile concentration, the progression is usually hexagonal, cubic, lamellar and inverted hexago- 
I in that order. It should be noted that all of these phases are not always present for particular 
iphiphilic systems. Disordered phases similar to the bicontinuous microemulsion have also 
en found. This disordered phase is often called the L 3 or sponge phase [32], This L 3 
lase is therefore made up of sheets of amphiphiles that are examples of amphiphilic bilayer 
Lembranes(see Fig.1.2). 

Although the generic behaviour described above is common to most types of amphiphiles, 
le shape, size, flexibility, chemical structure, etc. of the amphiphiles are known to play a part 
determining the details of the phase diagram as well as the morphology of the supramolec- 
lar aggregates. Surfactants can also have their behaviour affected by the addition of a fourth 
amponent: salt. Ionic amphiphiles have an electrostatic repulsion between the heads which be- 
omes screened in the presence of counterions when salt dissolves in water. Therefore one finds 
lat a quaternery system transforms from oil-in-water to bicontinuous and finally water-in-oil 
licroemulsion as the concentration of salt is increased [ 20 ]. 
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Figure 1.5: Theoretical description of amphiphile molecules on different length scales (after[19]). 


1.4 Models describing Amphiphilic Systems 


The theoretical description of amphiphilic systems depends on the length scale on which 
the phenomena are to be explained and predictions to be made, and on how much detail of 
information is expected to be extracted from the model. This has led to three different classes 
of models (see fig 1.5), which have been used in the last few years to study self-assembling 
amphiphilic systems [2]. 
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1.4.1 Microscopic models 

The natural starting point for a study of complex fluids is a microscopic model, in which the 
position of all the atoms of different molecules and their microscopic interactions are taken into 
account, as shown schematically in Fig 1.5. Such models have been used to study the structure 
of single micelles [95] , but unfortunately they are too complex to calculate the macroscopic 
properties of microemulsions etc. Thus, the models have to be simplified considerably to make 
them tractable for both approximate analytical calculation and Monte Carlo simulations. These 
simplified models are most often formulated on a lattice in the spirit of lattice gas models. The 
models which have been studied range from the Ising model with next-nearest-neighbour and 
multispin interactions [35, 101, 104, 105], and the Blume- Emery-Griffiths (BEG) model [36] 
with three spin interaction [97, 37] to models in which the amphiphile is described by a vector 
[64, 38] or by a pair of charges, which are equal in magnitude but opposite in sign [48, 49] or by 
3 short polymer chain [41-47, 60]. 

The Widom Model 

The simplest lattice model is the Widom model [35] where the amphiphile, water and oil 
nolecules all occupy bonds on a lattice, which is assumed to be square in two dimensions or 
I simple cubic in three dimensions. Each site of the lattice is occupied by classical Ising spins 
vhich can take one of two values. Si = ±1. Without loss of generality, water molecules occupy a 
)ond between two up spins, oil molecules occupy bonds between two down spins and amphiphile 
nolecules occupy bonds between one up and one down spin, the head obviously corresponding 
0 the up spin. The Widom Model is illustrated in Fig. 1.6. 
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Figure 1.6: A typical configuration of Widom model. 
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Figure 1.7: Origin of bending and different neighbors in Widom model 
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The effective Hamiltonian is, 

nn nnn nnnn 

H = -Jj2SiSj-M'£SiSk-2M SiSk ( 1 . 1 ) 

(ij) (ik) (ik) 

where nn refers to nearest neighbours (i.e. 1st neighbour), nnn next nearest neighbours (i.e. 2nd 
neighbour. Fig. 1.7(c)) and nnnn is the next to next nearest neighbour (4th neighbour in case 
of simple cubic lattice. Fig. 1.7(d), and 3rd neighbour in case of square lattice). J is positive 
(ferromagnetic) and M is negative (antiferromagnetic). The farther neighbour interaction was 
included to take into account the bending of the amphiphilic molecule at the interface. Fig 1.7(a) 
represent a flat interface. On the other hand whenever heads or tails of two amphiphilic molecules 
meet at the same lattice point bending of the interface occurs (Fig. 1.7(b)). 

The main advantage of this model is that it is simple; in fact it is essentially an Ising model 
with not only nn interactions but also farther-neighbour interactions of specific type. The Widom 
model can account for droplet and lamellar phases and a phase corresponding to bicontinuous 
microemulsion. More details of the model and results will be discussed in chapter 2. 

The Three Component Model 

Here oil, water and amphiphile all occupy sites on a lattice and are represented by the three 
states of a classical spin 1, where the spins can have a value Si = 0, ±1. The tendency of 
amphiphiles to go to the oil-water interface can be included by the addition of a three-body 
interaction in the Blume, Emery and Griffiths (BEG) model [36], [39]. The effective Hamiltonian 
is 


iJ = 


- Y. [75^5, • -f -F M{SfSj + 5i5j)] - L ^ 5^(1 - S‘^)Sk , 

{ij) (ii^> 


( 1 . 2 ) 
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ere {ij) is a sum over nearest neighbours and {ijk) is a sum over all sets of three connected 
;s. The three body interaction is antiferromagnetic between sites i and k, and favours lamellar 
nfiguration. Lamellar and disordered phases are observed in this model. 

The Alexander Model 

In the Alexander model [40], which is a mixed site-bond model, oil and water occupy sites 
a lattice. The “water-loving" head of each of the amphiphiles is assumed to be identical to a 
iter molecule and the “oil-loving” tail of the amphiphiles is assumed to be identical to an oil 
olecule. Therefore, every lattice site is occupied by either water or oil molecule. Those pair of 
dike nearest-neighbour sites (i.e., oil-water nearest-neighbour pairs) are identified as amphiphilic 
olecule which are connected by a rigid bond. Therefore this model is often referred to as the 
iecorated spin 1/2” model. The Hamiltonian of the system is taken to be [50, 51] 

H = -hJ2Si-JE3iSj-pY:'^is-v Y. (i-s) 

i (ij) {hi) {iiik} {{ijik}) {i,j) 

here Si = ±1 occupy sites i of a d-dimensional, hypercubic lattice and Tij = 0, 1 occupy 
s links < ij >; (z,j) are distinct, nearest-neighbor pairs of sites, and {ij,ik) and {{ij,ik)) 
-e, respectively, distinct pairs of links that meet at sites i at 90° and 180°. The field h is 
dated directly to the difference between the chemical potentials of oil and water; the exchange 
iteraction J is chosen to be positive so that it favours the phase separation of oil and water. 
j = 1(0) denotes presence (absence) of a surfactant molecule on link (ij), p is the chemical 
otential of surfactant molecules, and V and Lj parametrize the interactions of these molecules. 

1 is the strength of the surfactant-mediated interaction between nearest-neighbour spins. The 
base diagram of this model has also been investigated exhaustively. 



1.4 Models describing Amphiphilic Systems 


17 


Vector Models 

The microscopic models we have discussed so far do not include the orientational degrees of 
freedom of the amphiphile molecules and water and oil molecules are taken to be structureless. 
Without taking into account the orientational degrees of freedom of the amphiphile molecules 
and by treating water as structureless these models miss orientational bonding effects between 
water molecules and between water and amphiphile molecules. There are vector models where 
the amphiphile and/or water molecules are given orientations to more realistically simulate the 
directional nature of the interactions between the amphiphile and solution [64]. 

The Larson Model 

The Larson model, which was introduced originally by Larson and co-workers [41-44] and 
subsequently reformulated by Jan, Stauffer and others [45, 46], is an attempt towards more 
realistic modelling by taking into account the conformations of amphiphiles as polymer chains 
rather than merely single units. 

We first describe the original version of the model, as formulated by Larson and coworkers. 
The Larson model is simulated on a L x L square lattice (in 2d) or L x L x L cubic lattice (in 
3d). Each oil and water molecule can occupy one site of the lattice while the amphiphile, which 
is denoted [41] by T^Tlq, occupies (m + q) neighbouring sites connected by rigid bonds with the 
head containing q sites and the tail m sites. All sites of the lattice are occupied. Amphiphilic 
chains can not break but bend and move in the solvent. Each head unit is effectively identical to 
a water molecule so far as its interactions with other molecules are concerned, and similarly, from 
the point of view of inter-molecular interactions, each tail unit is effectively identical to an oil 
molecule. Each site of the lattice interacts with nearest neighbours and next nearest neighbours. 
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The effective Hamiltonian can therefore be given by 


H = -Y.^E^PiPi (1.4) 

ab ij 

ere {ij) is over nearest and next (diagonally) nearest neighbours and the sum {ah) is over the 
)e, ‘water’ or ‘oil’, of monomers. The strength of the interactions with nearest and next nearest 
ighbours are equal. The occupation variable, = 1 if a molecule of type a is on site i and is 
'0 otherwise. There are periodic boundary conditions in all directions. The energy parameters, 
vw, E-ujo 3nd Ego, are between water-water, water-oil and oil-oil respectively. One can define 
= {2Ewo — Eoo — Eyjyj)/2kBT which is the only interaction parameter in the model. It may 
considered a dimensionless inverse temperature. 

A typical amphiphilic molecule with q = l and m = 8 is shown in Fig. 1.8(a). The system 
allowed to change its configuration in three ways. First, a water and oil molecule can change 
isition irrespective of how far apart they are, second, an amphiphile chain can move by reptation 
;. it slithers along its length like a reptile [63], and the chain can twist by forming a kink when 
le of its sites moves to a diagonal neighbour; these ‘moves’ are shown in Fig. 1.8(b). The 
ains are self-avoiding, i.e. an amphiphile is not allowed to intersect itself or other amphiphiles. 

The system evolves using a Metropolis algorithm as follows: one of the allowed moves 
scribed above is ‘attempted’ and the corresponding change in energy, AE, that would be 
used by implementing the attempted move, is calculated using the Hamiltonian in equation 
.4). If AE < 0, then the move is accepted with probability unity and if AE > 0 then it is 
cepted with probability p = exp{-AE/kBT). This procedure was defined by Larson [41] as 
Monte Carlo Step. Larson performed simulations on lattices of size 40 x 40 and 100 x 100 
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Figure 1.8; (a) A model amphiphilic molecule, (b) Allowed moves for amphiphiles. 

lattices [41] in 2d and on 20 x 20 x 20 lattices [42] in 3d. 

The model has been used extensively by Larson [41-44] to study the phase behaviour of 
ternary mixtures for different amphiphiles, TiTii, 72712, The simulations were per- 

formed at a temperature corresponding to 0.6Tc where Tc is the critical temperature below which 
oil and water, in the absence of amphiphile, phase separate on this lattice. This choice of tem- 
perature corresponds to w = 0.5 in 2d and w = 0.1538 in 3d. 

Three-phase coexistence was observed in all cases except for the amphiphiles TiTii- Lamellar 
phases, hexagonal phases and cubic phases were observed. For example, in the case TiTi^, for 
amphiphile concentrations of = 30 — 40%, close-packed spherical structures are formed while 
for concentrations 4)a = 50 — 60% hexagonally packed cylinders were formed and at concentration 
(j)^ = 65—85% co-existence of cylinders and lamellae and above (f)a = 85% predominantly lamellar 
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ctures were observed. The microstructural transitions from less ordered to more ordered 
ses as the temperature was lowered was also studied. It was also observed for the 
cture that as the concentration of amphiphile was increased, the shape of micelles changed 
n spherical to cigar-shaped to 'worm-like'. It was observed that head-tail asymmetry such as 
Increased the tendency to form cigar-shaped or ‘worm-like’ micelles. 

Later Jan, Stauffer, and collaborators [45-47] reformulated the Larson model in terms of 
g-!ike spin variables, in the same spirit in which a large number of simpler lattice model had 
n formulated earlier for the convenience of calculation. This modified larson model will be 
:ussed in chapter 4 of the thesis. 

Moleculair models 

The Shell Amsterdam group [60] has developed a molecular model of the amphiphiles in 
jtion in which water and oil molecules are represented by spherical monomers interacting via 
icated Lennard-Jones potentials in the continuum. An amphiphile is modelled as a linear 
in of water and oil molecules. The forces are given by the gradient of the potentials and 
numerically integrating Newton’s second Law - the system evolves in real time towards its 
ilibrium state. Equilibrium configurations are calculated using the molecular Dynamics (MD) 
ulation technique. The advantage of this method is that the simulation results can be related 
the measured parameters for the real physical system. The disadvantage of this approach 
hat it is computationally very expensive, requiring long execution times. With present day 
iputers, one can normally only simulate nanoseconds making equilibrium configurations ex- 
nely difficult to calculate. Using this model, Smit and co workers have studied some particular 
sical realizations such as the clustering of amphiphiles at the oil-water interface. Because MD 
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is computationally expensive, they did not fully explore the phase diagram. 

1.4.2 Landau- Ginzburg models 

On a length scale somewhat longer than the microscopic one, a ternary amphiphilic system 
can be described by the local concentrations po{r), py,{r) and Pampi'!') of the three species of 
molecules, namely, oil, water, and amphiphile (see Fig. 1.5). The interaction between the 
molecules is now described by free energy functional F[po, Pw, Pamp]- This functional can either 
be obtained from a microscopic model by integrating out degrees of freedom on smaller length 
scales, or from symmetry considerations. The construction of Landau-Ginzburg model depends 
on the range of phenomena the model is expected to describe. If we want to study only ternary 
amphiphilic systems for which the water and oil concentration do not differ enormously, then a 
model with a single scalar order parameter is sufficient [65, 67]. On the other hand, if the whole 
range of concentrations between the binary and the ternary system is to be investigated, a model 
with two or three order parameter is required [68]. 

The Landau-Ginzburg free energy functional with a single scalar order parameter ©(r) is 
given by [65] 

F(<l>(r) = I (iV[c( -H gi^) ( -p /($) - //$] (1.5) 

with ^ < 0 and c > 0. The negative g tends to create the interface and positive c stabilize the 
system. Here the order parameter 'i'(r) describes the difference between the local densities of 
oil and water and p, is the chemical potential difference between oil and water. The amphiphile 
concentration does not appear explicitly in the free energy functional and should be considered 
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have been integrated out [96], In order to account for the three phase coexistence possible 
this ternary system, Gompper and Shick imposed the condition that the function / must be 
osen such a way as to have three minima. In the absence of amphiphile molecules /($) has 
'0 minima at $ = <^0 and # = which describe the oil rich and water rich phases, and 
a positive constant. When amphiphile is added, a third minimum of / appears at $ = 0, which 
scribes the disordered microemulsions phase. 

The free energy functional with two order parameter can be written as [66] 

F[#(r), p(r)] = J d^r[c(V$)^ — + ap^ - pp- + Fg (1.6) 

ere p(r) is the density of the surfactant molecules at the location r. c,ro,u, g,a and p are 
lenomenological coefficients, p is the chemical potential of the amphiphilic molecules, whereas 
.p is the difference of the chemical potentials of water and oil. Cubic terms have been omitted 
om above equation; however in more general situations such terms may play some role, g is the 
rength of the coupling between two fields p and The surfactant property of the amphiphiles 
taken into account through 

Fs = s J (frp{V^f ( 1 . 7 ) 

more detailed Landau-Ginzburg model with two order parameters was introduced by Chen et. 

. [69] for studying the equilibrium phase diagram. 

Finally models can be constructed in which both the density and the orientational degrees 
amphiphile are taken into account and in that case we have three order parameters in the 
mdau-Ginzburg free energy functional [68]. 
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1.4.3 Interface/Membrane models 

On a length scale which is even longer than that in the Landau-Ginzburg formulation, the 
detailed structure of the amphiphilic monolayer at the interface between oil and water (or that of 
a bilayer in water) can no longer be resolved. The interfaces (or “membranes") are now described 
as two-dimensional mathematical surfaces. By studying the "out-of-plane” thermal fluctuations 
of these amphiphilic membranes, the membrane theories provide a universal description for ternary 
and binary mixtures. At this level of description, the effects of the amphiphilic molecules enter 
only through the surface tension and the other elastic constants, e.g., the bending rigidity, of 
these surfaces. 

The nature of the out-of-plane thermal fluctuations of interfaces and amphiphilic membranes 
as well as their effects on the equilibrium conformations have been investigated extensively over 
the last ten years [106-108, 5, 109-111]. In the case of fluid-fluid interfaces in the absence of 
amphiphiles, these fluctuations are controlled by the interfacial tension. On the other hand, the 
effective interfacial tension of amphiphilic membranes can be vanishingly small and, therefore, 
their out-of-plane thermal fluctuations are controlled by the bending rigidity [111] . The interfaces 
can exhibit a transition from a “smooth” phase to a “rough" phase at a temperature known as the 
roughening temperature [112, 114] . Similarly, depending on the nature of the in-plane ordering, 
an amphiphilic membrane can exhibit a crumpled phase at sufficiently high temperatures. In the 
following we examine the different definitions of crumpling. 

For a two dimensional fluid membrane embedded in a three dimensional space, bending 
deformation are expressed in terms of curvature. We define a unit normal h{x, y) to the membrane 
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a point {x, y) with height z. The s-component and y-component of the normal being given 
Ux = —{%) and Uy = -(||)- The orientational correlation between the normals to the 
embrane at two different points decays exponentially and is given by 

< n(0).n(r) >'^ (1.8) 

lere Cp is the persistence length and is given by [26] 

(p ~ aexp{ATrK /SksT) . (1.9) 

t length smaller than Cp. the membrane is effectively flat, i.e. its normals at different points 
e well-correlated. Beyond the length Cp. the normals at different points become effectively 
Korrelated, and the membrane is said to be be crumpled. 

A second alternative, but equivalent, definition of crumpling is due to Helfrich [108, 111] 
hich we outline below. Having defined the normals one can define two principal curvature at 
iy point on the surface Ci = (^) and C 2 = (^) The curvature elastic energy density can be 
(pressed as a function of the two principal curvatures and is given by [108, 111] 

H = K/2{ci + C2)^ — Kco{ci -H C 2 ) -f KgCiC2, (1-10) 

here Cq is the spontaneous curvature, K{> 0) is the modulus of bending rigidity, and Kq is the 
astic modulus of Gaussian curvature C 1 C 2 . This Hamiltonian is rotationally invariant in thex-y 
ane. Let A = ht the basal area of the membrane in the x — y-plane. The instantaneous 
)nfiguration of the membrane is defined by the local height variable u{x,y) and u(^x,y) obey 
le periodic boundary condition. Fourier expansion of u{x,y) leads to 

^(^ =HKco5(y ■ f} + bgsin{q • f)] 

9' 


( 1 . 11 ) 
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Because of the periodic boundary condition, allowed values of g are given by (g^., qy) = (27r/X)(m, n) 
where the integers m and n cannot vanish simultaneously. Moreover, the sum is to be carried 
over only half the g- plane. For weakly deformed membrane (|Vw|) one has 


Ci+C 2 = Y^ q^[agCOs{g- f) + bgSin{q- r*)] (1-12) 

q' 

The total bending energy is given by 

+ (1.13) 

^ Q 

The absorbed membrane area per unit basal area is l/2(Vu)^ to lowest order in Vu. Integrating 
the local values gives the total absorbed area 

S,A = \Y.e{al + t:‘,)A (1.14) 

In the presence of non-vanishing surface tension <j the total energy of thermal undulation is given 
by 

= 7 E(<'«" + Kq')(a?, + 6j)/l (1.15) 

^ Q 

From equipartition theorem of energy the mean square amplitude is given by 


2 l 2 2kBT 

< o, >-< 5, >- 


(1.16) 


The sum in the above equation (1.15) can be replaced by integration in the g plane to obtain 
the average mean square amplitude of undulation and is given by 


.<-u2 > = y < >2 (117) 

= I (1/2) < a, >2 dq (1.18) 

2k bT r qdq 

Air J {aq^ -T Kq'^) 


(1.19) 
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r K = 0,a ^ 0 (only interfacial tension) 


<u^ > 


/ 


ksT 
2'Ka 
^ ksT 
2'Ka 


dq 

g 

ln{L/a) 


( 1 . 20 ) 


this case the membrane is logarithmically rough. For tr = 0, if 7 ^ 0 (only bending rigidity) 


<u^ > 


k^T r dq 

2kK J q^ 


keT 

4kK 


iHaf 


(1.21) 


The crumpling of the membrane can be seen from the excess area 5 A per unit basal area A. 
the absence of interfacial tension (a = 0) from the equations (1.14) and (1.16) we have 

6 A ksT , / -r I \ / \ 

T 4^ 

So in the thermodynamic limit L 00 , < 5A > /A 00 i.e. (excess area/basal area) 
verges and the membrane is said to be crumpled. Roughly the length scale at which < 5A > 
1 = 1 , the membrane looses its orientation and is said to be crumpled. This length can be 
entified with the persistence length Qp and solving eq. ( 1 . 22 ) we find 

. AkK, 

Cp = (1.23) 

lere K is the bending rigidity. This formula differs from that given by eq. (1.9) by a numerical 
:tor only. 


Next, we mention a third alternative definition of crumpling. If the interfacial width W 
ales as FF ~ with a > 0 , the membrane is rough. Moreover, if cr > 1 for i < g, where Cp 
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is the persistence length [26], then the membrane is crumpled on length scales much longer than 
^p\ thus, all crumpled membranes are rough but the converse is not necessarily true [109], 

All the above definitions for crumpling are known to be equivalent when applied to continuum 
models of amphiphilic membranes. 


1.5 Organization of the thesis: 

In this thesis we study interface/membrane models and microscopic lattice models but do 
not consider any Landau-Ginzburg model. So far as the membrane models are concerned, we have 
studied only the out-of-plane thermal fluctuations of some model membranes and investigated 
its consequences on the conformational properties of the membrane. We have not studied the 
self-assemblies and bulk properties of the amphiphilic systems using membrane model approach. 
Different bulk properties and phase behaviours of the amphiphilic systems have been studied with 
microscopic lattice models. With microscopic models we have investigated the self-assemblies of 
three different classes of amphiphilic molecules with increasing complexity; we begin with single- 
chain amphiphiles and end with gemini surfactant and, in between, present a brief discussion 
on double-chain surfactants to maintain the logical sequence of development of our ideas and 
methodology of modelling. 

The outline of the thesis is as follows: 

In the chapter 2 of the thesis we show that the oil-water interface is rough, but not crumpled, 
in the droplet phase of the Widom model which is the simplest microscopic lattice model of ternary 


microemulsions. 
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in chapter 3 we construct a class of discrete lattice models of membranes/interfaces where 
j can make a clear distinction between the interfacial tension and bending rigidity which can be 
witched off" independently of each other. More specifically, in chapter 3, we propose a d + 1- 
mensional restricted discrete solid-on-solid (RDSOS) model as well as unrestricted discrete 
did-on-solid (UDSOS) model for d-dimensional amphiphilic membranes where the bending 
iffness of the membrane is incorporated following the same prescription as in the Widom model, 
'hen the out-of-plane thermal fluctuations of the UDSOS model are controlled only by this 
ending stiffness we show that it is crumpled at all non-vanishing temperatures but the nature 
■ the crumpling is very different from that observed in continuum models. We also point out 
ime analogies between the relaxation to rough equilibrium states of these SOS models and 
jnequilibrium growth of some kinetically roughened model surfaces. 

In chapter 4 we study the roughening of the oil-water interface in the presence of single-chain 
irfactants and investigate the physical origin of the lowering of the oil-water interfacial tension. 

After studying the properties of the amphiphilic membrane at the oil-water interface in chap- 
rs 2,3 and 4, we investigate the air-water interface in the presence of single-chain amphiphiles 
chapter 5. In this chapter we introduce a microscopic lattice model of a binary mixture of 
nphiphilic molecules, of two different lengths, in a system where water is separated from the air 
lOve it by a sharp well defined interface. We demonstrate an entropy-driven phase segregation 
a direction perpendicular to the air-water interface when the initial total surface density of the 
nphiphiles is sufficiently high. We also investigate (a) the conformations of the amphiphiles, 

) the distribution of the sizes of the clusters of monomers belonging to the long amphiphiles 
well as of those belonging to the short amphiphiles in planes parallel to the interface, (c) the 
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effects of varying the lengths, total concentration and the ratio of the numbers of the two types 
of the amphiphilic molecules, the strength of the inter-monomer interactions and temperature. 

In the chapters 4 and 5 only single-chain amphiphiles have been considered although the gen- 
eral conclusions drawn from our MC studies are valid also for double-chain amphiphiles. In chapter 
6 we compare and contrast the self-assemblies of single-chain and double-chain amphiphiles, re- 
spectively, in water. In this chapter we study the self-assemblies of binary mixtures of single-chain 
and double-chain amphiphiles and compare these with those formed by the corresponding single 
species of amphiphiles. 

In chapter 7 we propose a microscopic lattice model of gemini surfactants. We show that 
these model gemini surfactants form long, thread-like and entangled micellar aggregates even at 
low concentrations where single-chain model surfactants form spherical micelles. We compare 
the morphologies of the micellar aggregates formed in our MC simulations with those observed in 
laboratory experiments and in earlier molecular dynamics simulations. We also study the variation 
of the critical micellar concentration of these model gemini surfactants with the variation of 
the (a) length of the hydrocarbon spacer (both hydrophobic and hydrophilic) connecting the 
two hydrophilic heads, (b) length of the hydrophobic tail and (c) the bending rigidity of the 
hydrocarbon chains forming the spacer and the tail. Some of the trends of variation are counter- 
intuitive but are in excellent agreement with the available experimental results. We also investigate 
the conformations and spatial organizations of the gemini surfactants at the air-water interface 
both in the dilute and semi-dilute regime. 

In chapter 8 we propose a scheme for MC simulation of ternary microemulsions in porous 
media (such as vycor glass). Finally, in chapter 9 we summarize our conclusions and outline 
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rections of future work. 



Chapter 2 


Oil- Water Interface in the Widom 

model is it smooth, rough or 

1 

crumpled? 


2.1 Introduction 

The Widom model (described in chapter 1) is. effectively, a spin-1/2 Ising model with not only 
nearest-neighbour ferromagnetic interactions but also with farther neighbour antiferromagnetic 
interactions. The Hamiltonian for the Widom model is given by equation (1.1). The interaction 
J is positive (ferromagnetic) whereas M is negative (antiferromagnetic); the latter arise from a 
particular prescription, suggested by Widom, for taking into account the bending energy of the 

paper based on this chapter has already been published by us ( Prabal K. Maiti and Debashish 
Chowdhury, J. de Physique. I, 5, 671 (1995)) 
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So far no attempt has been made to study the possibility of roughening/ crumpling of the 
nphiphilic membranes in the lattice models of complex fluids, e.g., monolayers at the oil-water 
iterface in ternary microemulsions. In this chapter we study the nature of interface in the Widom 
todel. The M = 0 limit of the Widom Hamiltonian (equation 1.1) is identical to the ordinary 
iing Model with nearest neighbour interactions on a simple cubic lattice; the interface between 
omains of up and down spins in this case is known to be rough at all T ^ 0. Therefore, our 
ivestigation of the interface in the Widom model in this chapter will show the effects of the 
pecific farther-neighbour antiferromagnetic interactions of the Widom model (non-vanishing M) 
m the roughened interface of the nearest-neighbour Ising model; we restrict our Monte Carlo 
MC) study to square lattice {d = 2) only. We have focussed attention only on that part of the 
)hase diagram on the j — m-plane where the droplet phase (i.e., the ferromagnetic phase in the 
pin language) is thermodynamically stable. 


2.2 Simulation 

We have used the standard techniques of Monte Carlo (MC) simulation. The system consists 
)f LxX{Ly+A) lattice where, by convention X and Y denote the horizontal and vertical directions, 
espectively. Each site of the lattice is occupied by a classical Ising spin (S,- = ±1). Periodic 
)oundary condition is applied in the X-direction whereas a fixed boundary condition is applied 
n the F-direction; the spins in the two uppermost layers remain "frozen” in the up (+) state 
/hile those in the two lowermost layers remain "frozen" in the down (-) state throughout the 
imulation. In some of our runs we used an alternative boundary condition where, instead of 
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fixed boundaries, anti-periodic boundary condition was applied in the y-direction. However, 
within the accuracy of our computation, the data for these two sets of boundary conditions were 
indistinguishable. Therefore, almost all of our production runs were made with fixed boundary 
condition in the K-direction. Application of these boundary conditions are used routinely in 
creating interfaces in the Ising model where majority of the spins in the upper half are in the up 
state whereas majority of the spins in the lower half are in the down state. For each set of values 
of j,m,Lx and Ly we began with an initial state in which all the spins in the upper half of the 
system were in the up state whereas those in the lower half were in the down state; this initial 
condition helps us in faster equilibration of the system. 


2.3 Results 

If the two nearest neighbour spins are antiparallel we call the bond connecting these two spins 
to be a "broken” one. After the equilibration of each configuration, we measured the number 
of broken bonds, Nb{y), as a function of the row-index y; in the droplet phase (ferromagnetic 
phase) one expects a maximum in Nb{y) near Ly/2 provided accurate data have been generated 
by averaging over sufficiently large number of configurations. In fig.2.2 we plot Nb{y) as a 
function of y for the parameter values j = 1.0 and m = 0; this corresponds to a two-dimensional 
Ising model with only nearest-neighbour interactions. The figures 2.2(a) and (b) correspond to 
Lx = 2000 and 20000, respectively. Similar plots were also obtained for j = 1.25, m = 0.0 and 
j = 1.6, m = 0.0. 

The full width at half maximum (FWHM) in the plot of Nb{y) can be taken as a measure 
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"igure 2.2: Two typical plots of A^(j/) for j = 1.0, m = 0.0 are shown. = 2000 in (a) and 
.0000 in (b). Ly — 200 in both (a) and (b). The data in (a) and (b) have been obtained by 
iveraging over 20 and 6 configurations, respectively. 
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Figure 2.3; The width of the interface of the Ising model with only nearest-neighbour ferromag- 
netic interactions (i.e., m = 0) are plotted against for three different values of j. 

of the width of the interface. Since the width W of the interface of two-dimensional Ising model 

with only nearest- neighbour interaction is known to vary as \/Lx with the length Lx we show in 

fig. 2. 3 that our corresponding data are, indeed, consistent with this behaviour. Moreover, note 

that for the same L, the smaller is the value of j the larger is the width of the interface, as it 

♦ 

should be. 

Next, in fig.2.4, we plot the width of the interface as a function of \/Xi for the same values 
of j as in fig.2.3, but now for non-zero values of m. Clearly, the interface is still rough in spite 
of the non-zero values of the farther-neighbour antiferromagnetic interactions. Moreover, note 
that, for given L and j , the interface is wider for non-zero values of m as compared to the Ising 
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interface for which m = 0. This arises from the fact that, for given j, the system is closer to the 
critical point for larger value of m as long as the system remains in the ferromagnetic phase. 

We conclude that although the farther-neighbour interactions in the Widom model give rise 
to a richer bulk phase diagram, as compared to that of the Ising model with nearest-neighbour 
interactions on the simple cubic lattice, these do not affect the interface roughening exponent, 
except for the fact that the interface becomes wider when these further-neighbour interactions 
are switched on. 

In contrast to our Monte Carlo simulation in 2d, Khang et. al [99] have studied the same 
Widom Hamiltonian (equation 1.1) by performing low-temperature series expansions in the three 
dimensional cubic lattice. In their study they used periodic boundary condition in the x and y 
directions and antiperiodic boundary conditions in the z direction. The interface at T = 0 was 
chosen to lie at z = 0 and in the zero temperature state, all spins are positive for z < 0 and 
negative for z > 0. Thus the ground state normalized layer density is p{z) = 1 for z < 0, 
associated with positive spins, and p(z) = 0 for z > 0, associated with negative spins. In their 
paper they also studied a solid-on-solid (SOS) model and a restricted solid-on-solid (RSOS) model 
(which will be discussed in next chapter) associated with the Widom Hamiltonian. To study the 
roughening transition they calculated moments of density gradient and the surface density. The 
moments of density gradient is defined as 

(z"”) = f:iP(^ - 1/2) - Piz + 1/2)1^" (2-1) 

oo 

which diverges as T ^ Tr. The other order parameter is defined as 

H = 

1 - p(l/2) - p(3/2) 


( 2 . 2 ) 
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3 hng et al obtained the low-temperature series of the the order parameters , {z^) and R in 
rms of the small parameters x = and y = for the Hamiltonian (1.1) as well 

, for both the SOS and RSOS models. They observed that roughening transition curve passes 
irough the region in which paramagnetic, ferromagnetic and modulated phases merge, and is 
arallel to the phase boundary of the ferromagnetic and modulated phases. It was also found 
lat crossover behaviour occurs in the roughening transition from the surface tension dominant 
;gime to the bending-energy dominant regime. 




Chapter 3 


Restricted and unrestricted discrete 
SOS models of interface with 
bending rigidity: equilibrium and 
non-equilibrium properties’ 


3.1 Introduction 


In the previous chapter we have studied the roughening properties of the oil-water interface 
in the Widom model in the droplet phase (ferromagnetic phase). Though the interface was 
found to be rough at all non-vanishing temperatures there was no evidence for crumpled phase 
^The following papers are based on this chapter: 

P. K. Maiti, D. Chowdhury and J. K. Bhattacharjee, Phy. Rev. E 54, 2670 (1996); P. K. Maiti and D. 
Chowdhury, IITK preprint (1997), submitted for publication 
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the interface. This might be due to fact that in the droplet phase of the Widom model the 
it-of-plane thermal fluctuations are controlled by surface tension as well as bending rigidity, 
'hen bending rigidity alone governs the out-of-plane thermal fluctuation, fluid like membranes 
e crumpled at all non-vanishing temperatures [108, 5, 109, 8]. One disadvantage of the Widom 
odel is that we can not switch off the surface tension J since that would destroy the droplet 
riase itself. Now the question arises whether it is possible to construct a discrete model of 
n interface where we can make a clear distinction between the interfacial tension and bending 
gidity and whether or not that model can exhibit a crumpled phase. 

The roughening and crumpling of isolated planar interfaces have been studied extensively in 
ie solid-on-solid (SOS) approximation [112, 114], The results obtained from the discrete and 
ontinuum models of interfaces have been found to be qualitatively consistent with each other 
L15, 113]. 

In this chapter we propose a rf-H 1-dimensional unrestricted discrete SOS (UDSOS) as well 
s restricted discrete SOS (RDSOS) model for d-dimensional amphiphilic membranes where the 
ending stiffness of the membrane is incorporated following a prescription suggested earlier by 
Vidom [35] in order to account for the bending rigidity of the membrane at the oil-water interface 
1 the lattice model of ternary microemulsions studied in chapter 2. 

In the case of the RDSOS model the restrictions imposed on the step sizes are such that a 
lear distinction between the interfacial tension and bending stiffness can be made. The motiva- 
on for imposing this restriction will be clarified later in the appropriate context. 
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Figure 3.1: A discrete model SOS interface 

3.2 The Model 

For simplicity, we formulate the model as a d-dimensional membrane embedded in a d + 1 
-dimensional space. For d = 1 we have a one dimensional membrane (curve) embedded in 
two-dimensional space and for d = 2 we have two-dimensional membrane embedded in three- 
dimensional space. Suppose the sites on a d-dimensional lattice are labelled by the integers i 
(f = 1, 2, 3, .., where L is the linear size of the lattice in each of the d-directions. The height 
of the membrane, measured from the reference line (in d = 1) or reference plane (in d = 2), is 
assumed to be given by a single-valued discrete height function; hi, the height at the z-th site, 
can take only integer values. A model discrete SOS interface has been shown in figure 3.1. If /ij is 
constant and independent of z then the membrane is flat everywhere. On the other hand, unequal 
heights at neighbouring sites gives rise to bending of the interface; these bends are identical to 
the bendings of the amphiphilic monolayer at the oil-water interface in the Widom model. 

The hamiltonian for the model membrane is given by [116] 

H = jY,\hi-hj\+K:Y,i^- (3.1) 

i>j i>j 

where 5hi,hj is Kronecker d-function and the constant JCis a measure of the bending stiffness of 
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e membrane whereas J is a measure of the interfacial tension. The summations in both the 
rms are to be carried out only over pairs of nearest neighbour sites on the d-dimensional lattice, 
ote that an energy cost JC is associated here with every bending of the interface at right angle 
/er a single lattice spacing. In the limit /C = 0, this model reduces to the standard SOS model 
• an interface with only interfacial tension. 

We denote the absolute value of the difference in the height variables at two neighbouring 
tes i and j, i.e., \hi — hj\, by the symbol |A/ij| which, in general, can be any positive integer 
, 1, N. If N is allowed to be infinitely large (i.e. the height variable at a particular lattice 
oint i can take any value irrespective of its neighbouring height variables) we call the model 
nrestricted discrete solid-on-solid (UDSOS) model. On the other hand if N is allowed to take 
inly a finite integer value we call the model as restricted discrete solid-on-solid (RDSOS) model. 

If \Ahi\ is allowed to take only the values 0 and 1 (i.e., N = 1), both J" term (surface 
ension term) and K term (bending stiffness term) are identical and the model may be regarded 
is a standard restricted discrete SOS model with effective interfacial tension J + fZ. Therefore, 
n order to make clear distinction between bending rigidity and interfacial tension the minimum 
illowed value of N should be 2. In all our computer simulations of RDSOS model 0, 1 and 2 
i/ere the allowed values of |A/ij|. 

We have applied periodic boundary condition in all the d-directions of the d-dimensional 
attice. We work with the parameters J = J/iksT) and K = IC/ikeT) where T is the tem- 
lerature and ke, the Boltzamann constant is chosen to be unity in our units. Note that if the 
lagnitudes of J and K are increased (decreased) simultaneously keeping the ratio .7/ K constant, 

; effectively mimics a decrease (increase) of the temperature for fixed J and )C. 
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In our MC simulation of the UDSOS model, the system is allowed to evolve, starting from a 
flat initial condition {hi = 0 for all i), following the standard Metropolis algorithm [117] : a site 
on the d-dimensional lattice is selected randomly and the corresponding height of the interface 
hi was allowed to change by letting d,- -> dj ± 1 certainly if AE < 0 and with a probability 
proportional to exp{—AE) if AE > 0, where AE is the change in energy that would be caused 
by the change in height variable hi. In the case of RDSOS model the algorithm for time evolution 
is identical to that for the UDSOS model except that none of those changes in the height variable 
which would violate the restriction \Ahi\ = 0, 1, ...,7/ is allowed at any step. 

The width , W, of the interface is defined through the relation W^{L,t) = [< {hi{t)— < 
h >)^ >]a„ where hi(t) (i = 1, 2, ..., is the height at the i-th site at time t, the symbol < X > 
implies average of X over all the lattice sites, i.e. < X >= Xi, whereas the symbol [T]a„ 
implies average of Y over a large number of MC runs. The quantity W{L,t) is monitored, for 
fixed L, as a function of time t] W{L,oo) is the width of the system in equilibrium. We define 
the exponents a and p as follows: W{L,t) ~ for fixed L and for very early times whereas 
W{L,oo) ~ in equilibrium. 

Another quantity that scales in the same way as the interface width is the height-height 
correlation function 

(7(0 = [< {hi^i - hif {I < L) (3.2) 

The roughness exponent can be determined from the relation 


C{1) cx P 


{I « L). 


(3.3) 
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We have also computed two other correlation functions in the case of UDSOS model in 
= 1. The correlation function G(x) is defined as 

G(x) = [(ho- <h>)(h^-<h>)U (3-4) 

the continuum formulation, De Gennes and Taupin [26] studied the correlation between the 
jrmals to the membrane at two-spatially separated points. We have introduced a novel discrete 
3 unterpart of this normal-normal correlation appropriate for our discrete model. Only three 
irections are possible for the normal to our discrete interface- it. can be directed either towards 
ft (— f) or right (x) or can be vertical. Therefore, we can define the "normal" as follows; 

Ui = -t-1 if hi > hi+i 

rii = 0 if hi = /ij+i 

m = -1 if hi < hi+i 

1 terms of these "normals” the normal-normal correlation function N{x) for our discrete model 
defined as 

N^{x') [t^O * ‘^xjav (3-5) 


1.3 UDSOS Model 

In order to test our algorithm and computer program we have calculated W{L,oo) as a 
inction of L first in special case d=l,K = Q,J= 1.5, i.e. for the standard SOS model. From 
? 3.2 we find that in this special case of J ^ 0, 1^(1, oo) oc VZ. in agreement with the 


)own result for the standard SOS model. 
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J=1.5,K=0.0 


j=1.5.K=1.0 
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Figure 3.2: The equilibrium width W of the interface plotted against for two different sets 
of values of the parameters J and K, namely, J = 1.5, FC = 0.0 and J = l.o,K — 1.0. The 


lines through the data points are the corresponding best fitted straight lines . 
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The advantage of the SOS model ( 3 . 1 ) is that one can "switch off J and thus, chosing non- 
;hing K, can have a discrete model of interface with purely bending rigidity only. Therefore, 
lext computed W{L,t) in the other special case d = I, J = 0,K ^ 0. We found that 
L,t) kept increasing for ever with the increasing time t. In fig.3.3 the solid curves represent 
L,t), plotted against for three different values oft, when K = 3.0; the corresponding 
a for JT = 1.0 are shown by the dashed curves in the same figure. Note that smaller K 
-esponds to an effectively higher temperature T. From Fig.3.3 it is clear that (a) for a fixed 
he interface is wider at a higher temperature than at a lower temperature after the elapse of 
; same time interval beginning from the flat conformation att = 0, and (b) at any given t, W 
ndependent of L but W{L,t) ^ oo as t — > oo for all L. 

Thus the actual area A of the membrane in equilibrium is infinitely large, if J = 0 but 
^ 0, at all r 7 ^ 0, even for all finite values of the basal area (length L \n d = 1 )! Does it 
cessarily imply that the membrane is crumpled at all non-vanishing temperatures? 

in order to answer this central question we shall now summarize the three different definitions 
crumpling. All these definitions yield identical results when applied to standard continuum 
Ddels of membranes and, therefore, usually regarded as equivalent. But, our discrete lattice 
adel, with J = 0 and K ^ 0, satisfies the criterion of crumpling set by one of these definitions 
it appears to be smooth when examined against the two other definitions, (f) If the ratio of 
e excess area 5A = A — Ap created by thermal fluctuations and the original basal area Ap of 
e membrane diverges in the thermodynamic limit the membrane is crumpled. This criterion of 
umpling is certainly satisfied by our lattice model ( 3 . 1 ) at all T 7 ^ 0 if J = 0 but AT 7 ^ 0 . (ii) 
the interfacial width W scales as W ~ 1 “ with a > 0. the membrane is rough. Moreover, if 



Figure 3.3: The width W of the interface plotted against for the parameter values J 0.0, 
K = 3.0 (full curves) and J = 0.0, K = 1.0 (dashed curves) at three different times, namely. 


t = 1000, t = 2000 and t = 5000 Monte Carlo steps. 
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> 1 for L < Cp. where Cp is the persistence length [26], then the membrane is crumpled on 
igth scales much longer than Cpi thus, all crumpled membranes are rough but the converse is 
)t necessarily true [109] . From fig.2, we conclude that, after the very early times (when no 
ear trends emerge due to the presence of transients), one can extract an effective roughening 
cponent a^ff for all longer times and one finds that cte// = 0. This, however, would imply 
lat at all times, except the very early stages, during its growth the membrane remains effectively 
nooth! [Hi) If the correlation between the normals to the membrane at two spatially-separated 
oints follows < n[f) •n(0) > ~ , then the membrane is defined as crumpled on length 

:ales much larger than Cp [26], But, we neither find any exponential decay in the ” normal-normal 
Drrelation function" N[x) nor observe any qualitative difference in N'[x) computed for K = 1.0 
nd iRT = 5.0 (see figs.3.4(a) and 3.4(b)), which correspond to an interface of fixed bending 
gidity K. at two different effective temperatures. 

In fig.3.5(a) and 3.5(b) we present the height-height correlation function C{x) for K = 
.0 and K = 1.0, respectively. The corresponding correlation functions G[x) are plotted in 
gs. 3.6(a) and 3.6(b), respectively. Noting that K = 1.0 corresponds to an effectively higher 
imperature than K = 5.0 for same fixed bending stiffness K, we conclude that the only effect of 
higher temperature in both these figures is the larger magnitude of the corresponding correlation 
inction arising from stronger thermal fluctuation. 

In order to understand these unusual features of the out-of-plane thermal fluctuations in the 
odel (3.1), let us compare schematic typical configurations for this model with those of standard 
mtinuum models. Since wiggles exist in the rough phase of the continuum model on all length 
ales (smaller wiggles are part of the bigger ones which, in turn, are parts of even bigger ones, 
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Figure 3.5: The height-height correlation function C(x) for (a) K = 5.0, (b) K - 1.0. 
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d so on) the width increases with increasing L. The absence of wiggles of all sizes leads to 
= 0 in our lattice SOS model. Moreover, it is obvious from these schematic representations 
typical configurations why the normal -normal correlation exhibits a periodicity in our lattice 
)S model. Note that the diverging area of the lattice membrane in equilibrium arises from the 
merging area (length in d = 1) of the vertical regions of the membrane. 

Next, we have computed the width of the interface as a function of L for the more general 
ise JT ^ 0, J yf 0 simultaneously. In order to test the effect of non-vanishing K on the 
'Ughness of the SOS interface we chose the same value of J, namely, J = 1 . 5 , as before, while 
' was chosen to be 1 . 0 . It is clear from fig. 3 . 2 , that the non-vanishing bending rigidity tends 
) suppress the out-of-plane thermal fluctuations thereby reducing the width of the interface but 
le roughening exponent remains the same as in the special case K = 0. In other words, in 
le general case, when both J and K have non-vanishing values, the interface is rough but not 
umpled. 

We shall now explain why "switching on" of the bending stiffness K, does not lead to any 
lange in the value of the roughening exponent of the standard UDSOS model with J" 7^ 0 , /C = 0 
Consider, for simplicity, d = 1. Then, J has the dimension of \Energy\l [Length] whereas /C 
js the dimension of [Energy], Therefore, there exists a crossover length = Kj J such that 
le equilibrium structure of the interface on length scales L » is governed by J . However, 
the special case J" = 0 , the equilibrium is governed solely by K. It is straightforward to see 
16 ] that if J" = 0 the interfacial width in equilibrium diverges. On the other hand, in the 
esence of non-vanishing J , the interface must have a finite width; for a given K, the larger is 
e J the smaller is the equilibrium width of the interface. Moreover, a larger value of /C, for 
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fixed J{> 0), leads to stronger suppression of the interfacial fluctuations. Nevertheless, if JT 7^ 0 
there is a very short early time regime where W{L,t) oc larger is the magnitude of J 

the quicker is the onset of the deviation from this law. These arguments can be easily generalized 
to d > 1. 

Now we study the dynamics of the UDSOS model interface analytically [118] in the special 
limit J = 0, but K ^ 0 by using a standard Master equation approach which is similar to that 
for the kinetic Ising models [120]. 

For simplicity, we consider only d= 1. 

hL;t) = (3.6) 

2=1 {/I'J 

where P{hx...hL,',t) is the probability of having the system in a state (fii, ..., /li) at time t. 
W{hi -4- h'i) is the transition probability, i.e., the probability of changing the state by changing 
the height variable at the z-th site from a value hi to a value /i[ = /ij ± d. The first sum on 
the right hand side goes over ail sites and the second sum goes over all the allowed values of h\. 
Expectation values of quantities A({/ii}) which depend on the height variables {/zj} are defined 
by 

< yl({/ii}) >= 2; yl({/n})F({/li}, 1) (3.7) 

(li) 

Therefore, we can write the equation of motion for the average value of < A({/ii}) > as [121] 

^ < ^(/Zi) >= - < - Mh'i)]W{hi h'i) > (3.8) 

dr {h-j 

To derive equation (3.8) form equation (3.7) we have used principle of detailed balance namely 


W{hi — )■ h'^P{h\...hi...hL\t) — W{h'i ^ hi)P{hi...h'i...hi]t) 


(3.9) 
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ice the transition probabilities W{hi -> /i-) depend on the neighbouring height variables hi+i 
d hi-i only we choose the probabilities W(hi h'^) to be of the following form [ 120 ] 


W{hi -> h[) = [1 - + 5k, - 5k-„h,^^ 




here 7 = tanh{PAE{hi,h[)/2) with AE{hi,h[) being the energy change of the system for 
— >• h[. For the case J = 0 and iiT 7 ^ 0 we have 


AE{hi,h[) = K{5hi,hi+i + 5hi,hi-, - 5f,>.,hi+i - 5h'.,hi-i) 


(3.11) 


ir simplicity, from now onwards, we shall consider only transitions between neighbouring states 
-7 hi ± 1. Some of the possible transition probabilities for different neighbouring height 
iriables have been tabulated in table 3.1. Associated changes in the energy also have been 
\/en in the table. 

The equation of motion for the average value < hi{t) > can now be written from equation 
1 . 8 ) as 

d 

~ < hi{t) > = — < [hi — {hi + 1)]1T^ {hi — > hi 4- 1) (3.12) 

+[hi-{hi-l)]W{hi-^hi-l)> 

= — < W{hi hi — 1) — W{hi hi + 1) > 

= P{hi...hi...hL; t)[W{hi hi - 1) - W{hi ^ hi + 1)] 

{A.} 

•om the table it is evident that for every possible transition there is a reverse transition. The 
lal configuration in one case is the initial configuration of the other case (see row 1 and row 12 
■table 3.1) whose energy contribution is opposite in sign. The contributions coming from these 
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Table 3.1: Some of the possible transitions and the corresponding transition probabilities for 
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cancels each other in the above sum. Summing over all possible transition we get 

= 0 

i.e. < lii(t) > = constant 


(3.13) 


Physically this means that as the system evolves with time the mean position remains fixed. This 
is in accordance with the results of our MC simulation reported earlier [116]. One can compare 
this result with the growth law of the standard SOS model with only interfacial tension [122] where 
< hi{t) >~ with an exponent p > 1/2 for d < 2 and < hi{t) >~ exp{-\{T)t), (A(T) > 
0; T < oo) for d > 2. 

The width of the interface is defined through 

W{t) =< < hiit) >f >=< hfit) >-< hi{t) >2 

To compute the growth law for the interface width we calculate < hj{t) >. From equation (3.8) 
we can write the equation of motion for < h'-{t) > as 

4:<hj(t)> = <[W{hi,hi + l) + W{hi,hi-l)]> (3.14) 

at 

+2 < hi[Wihi, hi + l)- Wihi, hi - 1)] > 


Summing over ail the possible transitions we find that the second term in the above equation 
vanishes and the first term gives a constant which is a function of temperature alone. So we have 

— < hi{t) > = constant (3.15) 

dt 

i.e. < hf(t) > cc t 

So, the width of the interface W{t) oc \/t; the results of our MC simulation (fig 3.7) [116] are in 
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Figure 3 . 7 : The square of the interfacial width plotted against time f for J = 0 . 0 , K = 3 . 0 . The 
inset shows the same data plotted against \/t to emphasize the fact that, when J = 0, JuT 7^ 0, 
the interfacial width diverges in the limit t —¥ 00. 

excellent agreement with this growth law. So for d = 1 we have /? ~ 1/2. Subsequently, we have 
also repeated the computation in d = 2 and have observed that ~ 1/2 also in d = 2 [ 119 ]. 

3.4 RDSOS Model 

We now motivate the imposition of restrictions on the value of N and the purpose of 
introducing the RDSOS model. If the discrete SOS models with bending stiffness are to be 
interpreted as models of amphiphilic membranes then it is useful to examine the physical meaning 
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of the unusual crumpling exhibited by the UDSOS model. For a fluctuating amphiphilic membrane 
there are two independent extensive variables, namely, the total area A of the membrane and 
the area Ap of the projection of the membrane on a reference base plane; for incompressible 
membranes A can change only if the membrane exchanges amphiphiles with a reservoir maintained 
at a fixed chemical potential ja whereas Ap can be altered by applying appropriate tension r. Thus, 
one can define four different ensembles [124] : 

(i) {A, Ap)-ensemble: the membrane is isolated and framed; 

(ii) (ju, Ap)-ensemble: the membrane is open and framed; 

(iii) {A, r)-ensemble: the membrane is isolated and unframed; 

(iv) (^, r)-ensemble: the membrane is open and unframed. 

All the models considered in this chapter describe framed membranes. Since the area A diverges 
in the UDSOS model with JT" = 0 in the limit t — ^ oo, irrespective of Ap, this corresponds, 
effectively, to an open framed membrane. On the other hand, in a realistic situation amount of 
amphiphiles available is likely to be finite and, therefore, it motivates us to impose restrictions 
on the allowed fluctuations of the membrane area A by restricting N to remain finite. 

Next, we describe the MC results for the RDSOS model. We have monitored the width of 
the interface W{L,t), for fixed L, as a function of time t \n d = 1 and d = 2, starting from 
an initially flat configuration, following the algorithm of the MC simulation described earlier, and 
then repeated the computation with different values of L. We carried out these computations for 
several different sets of values of the parameters J, K and N, namely, (i) J = l,iC = 0,N = 1, 
(ii) j = 0, AT = 1, iV = 1, (iii) J = l,iC = 0, Af = 2, (iv) J = 0, K = l, N = 2. 

The data for W{L, t) in the case (i) is identical to that for the case (ii) as there is no 
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difference between the effects of J and K if |Aki| = 0. 1, as already explained above. Comparing 

the data for the case (iii) with those in the cases (i) and (ii) we find that, for any given f, W{L,t) 

is certainly larger in the case (iii) as the height restrictions are less severe. 

From the best fit to the curves W{L,t) \we find that the exponent /3 is approximately 
0.25 not only in both the cases (i) and (ii) but also in the cases (iii)-(iv) thereby suggesting that 
~ 0.25 is independent not only of N but also of J and K so long as N is finite unless both J 
and K vanish simultaneously. This estimate is to be contrasted with /3 ~ 0.5 for UDSOS model. 

We illustrate our method of analysis with the data for the parameter values J — 0,K = 
= 2. The width W{L,t) of the RDSOS model for five different values of L are plotted in 
fig.3.8. The corresponding height-height correlation function C{1) of the RDSOS model for the 
same values of the parameters, namely, J = 0,ir = 1, = 2 has been plotted in fig.3.9. From 

the best fits in figs. (3. 8) and (3.9) we find /? ~ 0.25 and o: — 0.5. In fig. 3.10 we have plotted 
W{L,t)/L‘^ vs tjU using the data of fig.3.8 and the values of the exponents a, (3 estimated 
above. We find that all the curves of fig. 3.8 collapse onto a single curve in fig.3.10. thus 
supporting the Family- Vicsek scaling relation [123] 

W{L,t)sL‘‘f{^) (3.16) 


where z = aj^. Thus, we conclude, the exponent /? change its value from P = 0.5 to 0 = 0.25 
when constraints are imposed restricting the maximum allowed height difference Ahi to finite N. 
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Figure 3.8: The time evolution of the width W of the interface in the RDSOS model in d = 1 
for = 2 (i.e., \Ahi\ = 0, 1 and 2). The different symbols correspond to different system sizes 
L: L = lOO(x), T = 200(*), L = 400(n), L = 800 (■) and L — lOOOO(o). The dashed line has 
slope (3 = 0.25. The parameter values are J = 0, and K — 1. 
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Figure 3.10: The data of figure 3.9 rescaled according to equation (3.16) The different symbols 
correspond to different system sizes L: L = 100(x),L = 200(+),T = 400(n),I/ = 800 (■) and 
L = lOOOO(o). The slope of the dashed line is /3 = 0.25. 
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3.5 Comparison with Other Models 

We can now compare our UDSOS model with the model of Kahng et al [99]. Kahng et al 
considered the following SOS hamiltonian that is equivalent to Eq. (1.1) 

NN OnN LNNN 

Hsos — {j + 4M) ^ [hr — /ir'l + 2M ^ \hr — hr'l + M ^2 \^r ~ ^r'\ (3-17) 

(r,r') (r,r') (r,r') 

while the RSOS Hamiltonian is given by 

NN DNN LNNN 

Hrsos = Jj2\^>'r-hr'\+2MY: \hr - + M J2 - ^r'\, (3.18) 

(r,r') (r,r') (r,r') 

where NN, DNN, and LNNN mean nearest-neighbor columns, diagonal-nearest-neighbor, and 
linearly next-nearest-neighbor columns respectively. For both these models they calculated the 
moments of density gradient and the surface density as defined in chapter 2. The special SOS 
models considered by Kahng et al. do not include any bending energy of the type included in our 
model. Thus, interactions in their SOS models are similar to the first term on the right hand side 
of our equation (3.1). But, in addition to the interaction between the nearest-neighbour columns, 
they also allow interactions between farther-neighbour pairs of columns; the latter was motivated 
by the similar farther-neighbour interactions between the spin-pairs in the Widom model [35] The 
roughening of the interface in the Widom model has been described in chapter 2. 

One common feature of our SOS model and the Widom model [35] is that there are bending 
energies associated with bendings over one lattice spacing. Another lattice model, which asso- 
ciates bending energies over bendings of the interface over single lattice spacings, has been found 
to mimic reliably the energies associated with bendings of interfaces over truly molecular dimen- 
sions [128]. However, the reliability of such bending energies in mimicing the energies associated 
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with bendings over longer length scales remains to be explored in detail. 

We now compare the dynamics of evolution in our UDSOS and RDSOS models, from fiat 
initial configurations, with some surface growth models. The crucial difference between our SOS 
models and the surface growth models is that the former describe relaxation to the thermal equi- 
librium while the latter are never in equilibrium. Therefore, the exponent a describes equilibrium 
roughening of the UDSOS and RDSOS models studied here whereas it describes kinetic roughen- 
ing in the case of the surface growth models. The process hi hi + 1 in the SOS models studied 
here is the analogue of deposition of a particle on the column of particles rooted at the z-th site 
whereas the reverse process hi hi — 1 \n our SOS models may be interpreted as evaporation 
of a particle from the existing column at the same site. 

in the UDSOS model, except for an initial transient period starting from the flat conforma- 
tion, there are large number of sites where the increase of height does not cost any extra energy 
if vj = 0, )C ^ 0. However, in the UDSOS model, there is a fraction of sites where energy change 
is possible when height is increased even if J" = 0. Therefore, although the process hi-> hf + l 
in our SOS models may be interpreted as random deposition of particles the deposition proba- 
bility is not independent of the site. Moreover, no evaporation of particles is allowed in random 
deposition model of surface growth [123] In spite of these differences, the exponents a and 3 for 
our UDSOS model have the same values as the corresponding ones for the random deposition 
model of surface growth; for both /3 = 0.5, a = oo [116, 123] . 

The exponents a and (3 for the RDSOS model are different from the corresponding values 
for the restricted SOS model of surface growth [125] but are identical to the corresponding 
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exponents for the random deposition model with surface relaxation [125, 127] which is known to 
belong to the same universality class as that of the Edwards- Wilkinson model [123]. Mapping 
the restricted solid-on-solid model to a spin model and using the Master equation treatment for 
the spin systems, Plischke et al. [127] obtained /? = 0.25. This is identical to the value of /5 in 
our RDSOS model in spite of /C = 0 in ref. [127] because of the fact that K, does not affect the 
exponent P of the RDSOS model when N is finite. 


3.6 Conclusion 

In this chapter we have proposed a discrete model of an interface with a specific form of the 
bending energy; this bending energy is associated with bending of the interface over as short a 
length scale as the size of a molecule. In the case of UDSOS model when the out-of-plane thermal 
fluctuations of this interface are controlled only by this bending stiffness it is crumpled but the 
nature of the crumpling is very different from that observed in the continuum models. We have 
shown that although the conformation of the interface in this case satisfies one particular definition 
of a crumpled manifold, it does not satisfy two other widely used definitions of crumpling although 
all these three definitions are known to be mutually consistent when applied to the continuum 
models. We have studied some dynamical properties of a UDSOS model in the special case 
J = 0,/C ^ 0 in d = 1 by using a Master equation approach and the width of the interface W 
is found to vary with time t zsW = \/t. By identifying an appropriate crossover length, we have 
also explained why the roughening exponent a of the UDSOS model in the presence of non-zero 
J is independent of whether or not it has non-vanishing bending stiffness K,. 
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Introduction of a new RDSOS model is motivated through the examination of four different 
ensembles. By carrying out MC simulation of this RDSOS model, we compute an exponent _3 
which characterizes the time-dependence of its width during the very early stages of its dynamical 
evolution from a flat initial configuration. 

Finally, we point out some analogies between the relaxation from flat initial configuration 
to thermally roughened equilibrium conformation of these SOS models and the nonequilibrium 
growth of some kinetically roughened model surfaces. 



Chapter 4 


Oil - Water Interface in the 
2-dimensional Larson type of model 


4.1 Introduction 

, In ail the models discussed so far in this thesis bending energy of the amphiphilic mono- 
layer is introduced as a parameter. In real amphiphilic systems bending should arise from the 
microscopic interaction between different molecules and the conformational degrees of freedom 
of the amphiphilic molecules. Larson model discussed in chapter 1 explicitly takes into account 
the conformational degrees of freedom of the amphiphilic chain molecules and captures several 
essential features of these complex systems [3]. Here the system is modelled as a simple cubic 
lattice of size L^xLyX L,. Each of the molecules of water and oil can occupy a single lattice 
site. A smgle-choiu amphiphile can be described by the symbol TmMyU, where T denotes tail, 
369^199?)^ K. Maiti, Phys. Lett. A, 230, 
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% denotes head and M denotes the ’liaison or neutral part of the amphiphiies. m, p and q 
are integers denoting the lengths of the tail, neutral region and head, respectively, in the units 
of lattice sites. Thus, each amphiphile is a self-avoiding chain of non-overlapping sites of length 
i = [m + p + q). We refer to each unit of an amphiphile, regardless of whether it is a part of 
the head, tail or liaison, as a monomer. 

Jan, Stauffer and collaborators [45-47] reformulated the Larson model in terms of Ising-like 
variables, in the same spirit in which a large number of simpler lattice models had been formulated 
earlier for the convenience of calculations. The "water-loving” head group in this lattice model is 
assumed to be "water-like” and, similarly, the "oil-loving” tail group is assumed to be "oil-like”. 
In this model an Ising spin variable S is assigned to each lattice site; S', = 1 (—1) if the f-th 
lattice site is occupied by a water (oil) molecule. An amphiphile occupies several lattice sites 
each successive pairs of which are connected by rigid nearest-neighbour bond. If the j-th site is 
occupied by a monomer belonging to an amphiphile then Sj = 1, —1, 0 depending on whether 
the monomer at the jth site belongs to head, tail or neutral part. This version of the model takes 
into account more realistic features such as variation in the strength of the interactions between 
units, neutral parts of the chain, conformational interactions and repulsion between 'ionic' head 
groups [47]. 

The inter-molecular as well as intra-molecular inter-monomer interactions are taken into 
account through the interaction between the corresponding pair of Ising spins. Two spins in- 
teract provided the spins are located on the nearest- neighbour sites on the lattice. Thus, the 
Hamiltonian for the system is given by the standard form for an Ising spin system with only 
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nearest-neighbour interaction: 

h = -jY, SiSi. ( 4 . 1 ) 

<ij> 

We measure the temperature of the system in the units in which the monomer-monomer inter- 
action J = 1.0 and the Boltzmann constant Ub = 1-0. The models were simpler than that 
developed by Larson in that only nearest neighbour interactions were allowed. 

The moves allowed for the amphiphiles in the model are as follows: (i) reptation, (ii) Kink 
movement, (Hi) spontaneous chain-buckling: a portion in the middle of the chain is allowed to 
buckle spontaneously (figure 4.1(a)). (iv) pull: an amphiphile which is not fully extended because 
of the existence of a buckle is allowed to be pulled out so as to wash out the buckle (figure 4.1(b)). 
(v) wag: the end of the tail of the amphiphile are allowed to ‘wag' in a manner shown in figure 
4.1(c), where the new position of the end monomer is chosen randomly from among the nearest- 
neighbour sites, which is occupied by oil/water. The moves allowed for the amphiphilic molecules 
in the algorithm are related to similar moves introduced earlier in the context of polymers [138, 
139]; the "kink” move in our algorithm is identical to the "normal-bead motion" in ref [138]. 
Moreover, the ’ spontaneous chain-buckling” in our model not only generates new "buckles" or 
enhances the hump of existing buckles but also leads to the "crankshaft motion” [138] of existing 
buckles 

Standard "Metropolis algorithm” is implemented. Each amphiphile is allowed to try each of 
the above mentioned moves once during each MC step. Time is measured in the units of the MC 
steps per molecule. 


With the above model Stauffer, Jan and Pandey [45] studied the interfacial energy of an 




Buckling Pull Wag 

® Head o ^ Tail # Water or Oil 

Figure 4.1: New moves introduced in reformulated version of Larson model 


oil-water interface as a function of amphiphile concentration. They used amphiphiles of structure 
'HiNiTi with the spin of the head S = +1 and the tail 5 = — 1 and looked at both two and three 
dimensional interfaces. Glauber Kinetics and only the reptation movement for the amphiphile 
was used. The interfacial energy was found to decrease as the concentration of surfactant was 
increased and was close to zero at a concentration of (pa = 20%. It was also observed that 
the spontaneous magnetization Y^iSi, at a fixed temperature below the critical temperature 
Tc = 4.51 J of the pure Ising model, was destroyed by the presence of amphiphile showing that 
the amphiphiles made oil and water miscible below th&ir normal miscibility temperature, explicitly 
demonstrating three-phase co-existence. With the above reformulation of the Larson model a 
lattice model of double-chain amphiphiles and gemini surfactants can be developed which has 
been described in chapters 6-7. 

In this chapter we study the roughening properties of the interface in the above formulation 
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of the Larson model. We also investigate the physical origin of the lowering of oil-water interfacial 
tension by the single-chain surfactants in the reformulated Larson model. For these purpose we 
have computed the interfacial width in this model for different system sizes, temperatures and 
amphiphile concentrations by carrying out Monte Carlo simulation. 


4.2 The Model 

The system is modelled as a 2-dimensional square lattice of size x Ly. Each lattice 
site is occupied by one of the three molecular species, namely, either oil or water or a monomer 
belonging to an amphiphile. For simplicity we take the amphiphiles to be 'H\T\ i.e. amphiphiles 
molecule consist of a single head unit and a single tail unit. This model is similar to the model 
proposed by Kawaktsu and Kawasaki [129] for the study of phase separation process in immiscible 
binary mixture with a surfactant. A somewhat more realistic spin model of this type has been 
used to study the equilibrium properties of microemulsions by Chen et al [50]. 

The initial states are always constructed in such a way that water occupied the upper part 
while oil occupied the lower part of the lattice. The amphiphiles are put randomly at the interface. 
Periodic boundary condition is applied in the A^-direction whereas a fixed boundary condition is 
applied in the 1 -direction. The molecules in the two uppermost layers remain water while those 
in the two lowermost layers remain oil throughout the simulation. This initial condition helps us 
creating an interface in equilibrium. Liverpool and Bernardes [57], in their study of the behaviour 
of membranes with surfactants of type and UzTz without head-head repulsion, also 

found that a monolayer at the oil water interface could be stabilized by confining walls. During 
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the updating of the states of the system in our MC simulation we used Glauber dynamics (flip 
dynamics) for water and oil and the moves mentioned in section 4.1 for the amphiphiles. We 
measure the temperature of the system in the unit of J and the Boltzmann constant ks = 1-0. 
We follow the standard Metropolis algorithm, explained chapter 1. 

The interface profile was determined at each selected time as a density profile p{y) - which 
is the fraction of water (oil) molecules in a particular lattice plane (i.e. fixed y). The width of 
the interface was determined through [115, 130, 101] 

- j y‘^w{y)dy - ( J yw{y)dyf (4.2) 

with the statistical weight given by the density gradient 

^iv) = ~ ^ 2 ) (4.3) 

Here pi and p 2 are the densities in the up (water rich) and down (oil rich) phase far away from 
the interface and the buffer planes at the ends of the lattice. All these integrals run from the 
bottom to the top of the lattice. In our discrete lattice, the integrals are of course replaced by 
sums. 

Once the equilibrium is achieved, after every 1000 MC steps per site width of the interface 
is calculated and then averaged over 15 such runs to get the thickness W. 

4.3 Results 

First we have calculated the the width of the interface when there is no amphiphiles present 
in the system. This corresponds to the case of two dimensional Ising model. Width has been 
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Figure 4.2: The width of the interface of the two dimensional Ising model is plotted against 
for two different values of J : J = 1.5(*) and J = 1.8 (□). 

calculated for two different temperatures J = 1.5 and J — 1.8. We show in figure 4.2 that our 
data are indeed consistent with the known result that the width (H^) of the interface of the two 
dimensional Ising model is proportional to v^. Moreover note that for the same value of Lx, 
the smaller the value of J , the larger the width of the interface as small J corresponds to higher 
temperature. 

Next in Figure 4.3. we plot the width of the interface as a function of same 

values of J as in Figure 4.2. but now for different non-zero concentration of amphiphiles 
Clearly the interface is still rough in the presence of amphiphiles. Moreover note that for a given 
Lx and J, the interface is wider for higher concentration of amphiphiles = 80%) compared 
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Figure 4.3: The width of the interface of the Larson model for different concentration {<pa) of 
amphiphiles for different values of J is plotted against y/IZ: J = 1.5 and (j)a = 50% {*), J = 1.5 
and (pa = 80% (□), J = 1.8 and (f)a = 50% (■) and J = 1.8 and (f)a = 80% (o). 
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Figure 4.4; (a) Picture of the interface when no amphiphiles are present, (b) interface after a 
single spin flip, a denotes water molecule and . denotes oil molecule. 


to the interface width for lower concentration of amphiphiles (pa = 50%. 

The increase of interface width with increasing concentration of amphiphiles is a consequence 
of the lowering of the oil-water interfacial tension by the amphiphiles. This reduction in surface 
tension has also been observed by Stauffer et. al [45] in a ternary mixture of oil, water and 
amphiphiles. In their simulation they used amphiphiles of chain length 6. However, as is evident 
from our simulation, this lowering of interfacial tension and consequent widening of the interface 
width can be captured with amphiphiles of shortest possible length, namely, two monomers. 
Therefore it is the amphiphilicity of the surfactant which is essential for the reduction of interfacial 
tension. 


We shall now argue that the lowering of the interfacial tension by the amphiphiles in the 
model studied here is of entropic origin. 

Consider Figure 4.4(a) and Figure 4.4(b). In Figure 4.4(a) we have created an interface 
without any amphiphiles. Figure 4.4(b) gives the picture of the interface after a single spin flip 
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Figure 4.5: (a) Picture of the interface when one single amphiphile is present, (b) interface after 
the amphiphile has made a reptation movement, a denotes water molecule and . denotes oil 
molecule. A and - are the head and tail of the amphiphile molecule respectively. 

and also the value of p{y). The weight factor (Equation 4.3) in discrete form can be written as 


w{y) 


p(.y - 1 ) - pjy + 1 ) 

2(pi — P 2 ) 


(4.4) 


Here pi = 1 and p 2 = -1. The width of the interface given by Equation. 4.1 can be written for 
the discrete case as follows 


COL COL 

~ 'y (l/m yboundary) ^ (l/m y boundary)'^ iym)') • 


(4.5) 


m=l m=l 

where yboundary givcs the position of the interface and is given by COL/2. Here yboundary = 2. 
Substituting the values of w{ym) (using Equation 4.4) in Equation 4.5 we find width of the 
interface to be 0.1 Now consider Figure 4.5(a) and Figure 4.5(b). In Figure 4.5(a) we have 
created an interface with a single amphiphile. Figure 4.5(b) gives the picture of the interface 
after the amphiphile has made reptation movement in one of the possible directions and also the 
value of piy). Reptation movement in two other possible directions also give the same value 
of p(y) Substituting the values of w{ym.) (using Equation 4.4) in Equation 4.5 we find width of 
the interface to be 0.2. So we see that the width of the interface has increased when there is 
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amphiphile movement. Amphiphile movement gives rise to larger variation in p{y) as compared 
to the single spin flip. This is equivalent to the statement that the the presence of amphiphiles 
reduces the tensions of the oil-water interface and the decrease in interfacial tension allows a larger 
width of the interface. This reduction in the tension is due to the fact that as the interfacial 
area is increased the translational entropy of the amphiphile molecule is increased (for a fixed 
number of amphiphile molecules). This increase in entropy lowers the system free energy and 
thus reduces surface tension from its bare value [142]. 

4.4 Conclusion 

In summary we have shown that the width of the oil-water interface in the Larson model 
varies as Moreover, the wider interface in the presence of amphiphiles in this model is a 

consequence of the lowering of the oil-water interfacial tension by the surfactants; we have argued 
that this reduction of surface tension in this model is of entropic origin. 



Chapter 5 


Out-of-plane phase segregation and 
in-plane clustering in a binary 
mixture of amphiphiles at air-water 
interface' 

5.1 Introduction 

After studying the properties of the amphiphilic monolayer at the oil-water interface in 
chapters 2,3 and 4 we now investigate the air-water interface in the presence of single chain am- 
phiphiles. The spatio-temporal organization of amphiphiles in a system where water is separated 

paper based on this chapter has already been published by us (Debashish Chowdhury, Prabal K. 
Maiti, S. Sabhapandit and P. Taneja, Phys. Rev. E, in press 
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from air above it by the air-water interface has been studied by Chowdhury [131]. This is different 
from the oil-water interface in a microemulsion, because, unlike oil and water in a microemul- 
sion, here the air and water are not allowed to get dispersed in each other. Consequently it is 
sometimes referred to as a pseudo-ternary system. 

The Hamiltonian is the same as Equation 4.1 though air was represented by Si = -1. 
Systems of sizes Ly x were studied. The system consisted of Lair layers of air, lying on 
top of Lz — Lair layers of water, with the amphiphiles initially distributed at the interface. The 
amphiphiles were of the type '^. 2 % where n = 5 or 15. The moves of the molecules of water, 
air and amphiphiles are slightly modified so as to reflect the the fact that air and water were not 
allowed to exchange position as dispersions of air and water inside each other is not possible. 

In addition to the three movements described earlier, an extra movement of lateral diffusion of 
chains along the interface was allowed. Chowdhury studied the behaviour of a small concentration 
of chains placed at the air-water interface. The largest system size studied was 30 x 30 x 100. 
At low temperatures, they all remained at the interface but as the temperature was increased 
more and more of the amphiphiles moved to the bulk signifying an increase in the solubility with 
temperature. There was a dramatic increase in solubility over a narrow range of temperature. 
This trend is consistent with the earlier experimental observations. It was also observed that at 
any given temperature, the solubility was lower for amphiphiles with longer hydrophobic chains 
than for those with shorter hydrophobic chains, also in agreement with experiment. 

Chowdhury et al. [86] also studied the case where the polar heads of amphiphiles are either 
grafted or adsorbed on a solid substrate. When they are grafted, they are fixed to a particular 
point on the surface and can only move by chain conformational rearrangements. If the chains 
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are adsorbed, they cannot move away from the substrate, but they can diffuse laterally on the 
surface. Experimentally this has been extensively studied in the context of Langmuir-BIodgett 
films [85]. 

Chowdhury et al. studied several properties such as the number of monomers at depth z, the 
average vertical extent of amphiphiles, information on the average orientation of the amphiphilic 
molecules and that of the bonds connecting the monomers. Chowdhury et al. [86] observed, for 
QTCift&d and ddsovh&d amphiphiles, that the orientational order increased with surface density in 
agreement with experiment [85], theory [87] and Molecular dynamics simulation [88]. 

In this chapter we are interested in binary mixtures of single chain amphiphiles in water 
and at air-water interface. We introduce a microscopic lattice model of a binary mixture of 
amphiphilic molecules, of two different lengths, in a system where water is separated from the air 
above it by a sharp well defined interface. We demonstrate an entropy-driven phase segregation 
in a direction perpendicular to the air-water interface when the initial total surface density of the 
amphiphiles is sufFiciently high. We also investigate (a) the conformations of the amphiphiles, 
(b) the distribution of the sizes of the clusters of monomers belonging to the long amphiphiles 
as well as of those belonging to the short amphiphiles in planes parallel to the interface, (c) the 
effects of varying the lengths, total concentration and the ratio of the numbers of the two types 
of the amphiphilic molecules, the strength of the inter-monomer interactions and temperatures. 

We shall compare the phenomena observed in our computer simulation with some similar 
phenomena in grafted polymers and compare the underlying physical mechanisms in the two 
cases. 

The chapter is organized as follows: The model and the characteristic quantities of interest 



Out-of-piane phase segregation and in-piane clustering in a binary mixture of amphiphiles at 
84 air-water interface 


are defined in the sections 5.2 and 5.3, respectively. The concentration profiles of chemically 
identical amphiphiles of two different lengths are presented in section 5.4. The in-plane clustering 
of the amphiphile monomers are studied in section 5.5. Conformations of the amphiphiles are 
investigated in section 5.6. Results for chemically different amphiphiles are summarized in section 
5.7. In section 5.8, we compare our results with laboratory experiments and binary mixture of 
grafted polymers. Finally, in section 5.9, we present a summary of our results. 


5.2 The Model 

The model developed is based on the reformulated Larson model of microemulsions discussed 
in chapter 4. The system is modelled as a simple cubic lattice of size Lx x Ly x L~ where the 
horizontal air-water interface is parallel to the A'l'-plane and the vertically direction 

is chosen as the -f^’-axis. 

Since we are studying a binary mixture of the amphiphiles, let us represent them by the two 
symbols Hq^Tnn for the long and short amphiphiles, respectively, so that the lengths 

of the long amphiphiles is Li = qc + rrii and that of the short amphiphiles is L, = q, + rn,. All 
our investigations have been carried out for a fixed value of qi, = g, = 2. 

Notice that the monomers of the same chain as well as different chains are not allowed 
to occupy the same lattice site; this represents a hard-core intra-chain as well as inter-chain 
repulsion for chain separations smaller than one lattice spacing. Moreover, at any non-vanishing 
temperature, during the out-of-line thermal fluctuations of the chains, the hard-core repulsion 
leads to stenc repulsion between the chains. Furthermore, in the present formulation, there are 
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no potential energies associated with the torsion of isolated amphiphiles in water. 

Since both types of amphiphiles are represented by the Ising spin variables the amphiphiles 
may be regarded as chemically identical if the strength of the interaction J is also assumed to 
be identical for both types of amphiphiles. However, chemically different types of amphiphiles 
can be modelled by taking different values of the interactions Jss and Jls, where Jll is the 
interaction between two monomers both belonging to the long amphiphiles, Jss is that between 
two monomers both belonging to short amphiphiles while Ji,s is that between two monomers one 
of which belongs to a long amphiphile and the other to a short amphiphile. 

In our model here, chemically different amphiphiles can be described by first associating a 
"strength" Ci with the z-th site (z = 1, 2, ..., N). If the z-th site is occupied by water or air then 
Ci = 1 whereas Ci = Ci{> 0) if the z-th site is occupied by a monomer belonging to a long 
amphiphile while Ci = Cs(> 0) if the z-th site is occupied by a monomer belonging to a short 
amphiphile. In this case, the Hamiltonian for the system is given by 


H = -J'£ CiCjSiSj (5.1) 

<ij> 

Here we study the binary mixtures of chemically identical as well different amphiphiles. 

In each of the computer experiments, the initial state is constructed in such a way that air 
occupies the upper part while water occupies the lower part of the lattice. We begin with a 
mixture of Ng, long amphiphiles and Ns short amphiphiles; the heads of all the amphiphiles were 
initially put randomly on the sites in a horizontal lattice plane immediately below the air-water 
interface in a lattice of size Lx x Ly x Lz 3nd each of the amphiphiles was fully extended in the 
vertical direction. The system is then equilibrated at temperature T using the following dynamics. 
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during the updating of the states of the system in our MC simulation, air and water molecules 
were not allowed to exchange position as dispersion of air and water inside each other is not 
possible in our model. 

The moves allowed for the amphiphiles in our model are as follows :(i) reptation (ii) spon- 
taneous chain-buckling, (iii) Kink movement and (iv) Lateral diffusion at the interface: those 
amphiphiles whose heads are located no deeper than the molecular layer at the interface are 
allowed to move laterally where one of the four possible directions is chosen randomly and each 
of the units of the amphiphile is moved in that direction by one lattice constant. Each of these 
moves is possible only if the new positions are not occupied by any other amphiphile. 

We follow the standard "Metropolis algorithm": each of the above mentioned moves takes 
place certainly if AE < 0 and with a probability proportional to exp{—AE/T) if AE > 0 , where 
AE is the change in energy that would be caused by the proposed move of the amphiphile under 
consideration. Each amphiphile is allowed to try each of the above mentioned moves once during 
each MC step. Time is measured in the units of the MC steps per molecule. 

The largest systems simulated were 60 x 60 x 100. However, the data for smaller systems, 
e.g., 30 X 30 X 100 do not differ significantly from those for these large systems. Therefore, 
most of the production runs were made for 30 x 30 x 100 systems. The data were averaged over 
sufficient number of different initial configurations; typically ten configurations were found to be 


adequate. 
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5.3 Definitions of characteristics quantities of interest 


Now we define all the quantities we intend to compute for these models through MC sim- 
ulation. We compute all these characteristic quantities for the long amphiphiles and the short 
amphiphiles separately. 

The gross features of the spatial organization of the constituent molecules can be expressed 
through the equilibrium profiles of concentrations of the amphiphiles in the Z-direction, i.e., in 
the direction perpendicular to the air-water interface. So far as the concentration profiles of the 
amphiphiles are concerned, one can calculate two different quantities: at each layer one can count 
the number of monomers or one can count just the number of hydrophilic heads in that layer. 

It is straightforward to see that at any non-vanishing temperature, the amphiphiles cannot 
remain fully extended. One of the quantities that characterize the average size of the amphiphiles 
is their average radius of gyration [Rg] which is defined by 

( 5 . 2 ) 

1 

Rem = 7-E^ 

i=l 

where Ri{j) is the position of z~th monomer of the j-th amphiphile and Rem the corresponding 
center of mass of a single amphiphile. The symbol [....] shows the ensemble averaging over 
different configurations. The symbols La 3nd Na denote, respectively, the length and the number 
of the amphiphiles of the particular species under consideration; for the computation of the 
radius of gyration of the long amphiphiles we take La = Li, Na = Ni and for that of the short 
amphiphiles we take La = Ls, Na = Ng. 
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Next, we introduce a measure of the average projection of the amphiphiles on the ^-axis 
at a given instant of time t as follows; we calculate AZ{t), the difference of the instantaneous 
values of the ^-coordinates of the highest and the lowest monomers (i.e., the monomers with the 
largest and the smallest values of z-coordinates), on each amphiphile and average it over all the 
amphiphiles of the same species present in the system thereby getting < AZ{t) >. The data 
are, then, averaged over a sufficiently large number of runs thereby getting the configurationally- 
averaged value [< AZ{t) >]. In order to get a more detailed information on the distribution of 
vertical extent of the amphiphiles we simultaneously also calculate P{AZ), the fraction of the 
amphiphiles with vertical extent AZ at equilibrium. 

Another related quantity of interest is [< Rz(t) >], the average instantaneous projection 
of the vector R{t) on the Z-axis where R{t) for a given amphiphile is the vector from the first 
monomer on its head to the last monomer on its tail. Note that, in equilibrium, the first monomer 
on the head and the last monomer on the tail of an amphiphile do not necessarily have the largest 
and smallest z-coordinates. Therefore, in general, [< R^ >] is smaller than [< AZ >]. Moreover, 
AZ is always positive whereas, in general, for amphiphiles in water, [131], can be both 

positive and negative. Therefore, at every temperature T the equilibrium value of [< R,{t) >] has 
been used as the order parameter for measuring the average orientational order (^P) exhibited 
by the "head-to-tail” vectors at the temperature T. 

A detailed information on the orientation of the individual bonds between the monomers of 
the amphiphiles is contained in the so-called ’’bond orientational order parameter”, defined as 

$(p) = (1/A^,) ;^(3 < (2,-p+i - > _i)/2 (p = 

i=l 


(5.3) 
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where z = 1, 2, N, {N^ = for the long amphiphiles and iV, = N, for the short amphiphiles) 
refer to the amphiphiles and refers to the z-coordinate of the p-th monomer on the i - th 
amphiphile. Thus, is the difference in the heights of the two end points of the 

p-th bond connecting the p-th and p + 1-th monomers of the i — th amphiphile. For example, 
<|>(p) = -0.5 describes a situation where the p-th bond, on the average, lies in a plane parallel 
to the air-water interface whereas #(p) = 1 corresponds to a vertical p-th bond. The factors 3 
in the numerator and 2 in the denominator in equation (4) have been chosen so that ^{p) = 0 
in case of complete randomness in the orientation of the corresponding bond. In this chapter we 
shall study the effects of varying the parameters L( - L^, — Nt/N^ T, Q and Cj etc. 

5.4 Concentration profiles for chemically identical am- 
phiphiles of different lengths 

We begin with a mixture of 405 long amphiphiles and 405 short amphiphiles, i.e., Ni = 405 = 
Ns and the size of the lattice is 30 x 30 x 100. The system was equilibrated at a temperature 
T = 2.5 using the dynamics described in section 2. The concentration profile of the hydrophilic 
heads of the long as well as the short amphiphiles are shown in fig.5.1. The two profiles are clearly 
shifted with respect to each other thereby establishing a vertical phase segregation of the long 
and the short amphiphiles. We shall refer to this computer experiment as the first experiment. 

Next we carry out further experiments to establish the physical reason for this phase seg- 
regation. We reduce the concentrations of both types of the amphiphiles by a factor of 9, but 
keeping the relative concentration ratio of the two types 1 : 1, as before, without changing the 
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Figure 5.1; The equilibrium concentration profile of the hydrophilic heads perpendicular to the 
air-water interface. The system consists of 30x30x100 lattice sites and contains a mixture of 405 
long amphiphiles of length 17 and 405 short amphiphiles of length ( equilibrated at temperature 

T = 2.5. The solid lines correspond to the long amphiphiles whereas the dashed lines correspond 
to the short amphiphiles. 
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Figure 5.2; Same as in fig.5.1 except that there are 
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size of the system. In other words, we now repeat the first experiment with 45 long amphiphiles 
and 45 short amphiphiles in a system of size 30 x 30 x 100 at T = 2.5. Interestingly, unlike 
fig.5.1, no vertical phase segregation is now observed in the equilibrium profiles plotted in fig.5.2. 
We refer to this experiment as the second experiment. 


Note that the crucial difference between the first and the second experiments lies in the 
different initial surface densities of the amphiphiles: the total number of amphiphiles in the first 
experiment is 810 whereas that in the second experiment is 90 although the total number of sites 
available in the plane immediately below the air-water interface is 900 in both the experiments. In 
other words, the initial surface density of the amphiphiles in the first experiment is 90% whereas 
that in the second experiment is only 10%. Does it imply that the vertical phase segregation 
observed in the first experiment occurs only if the surface-density of the amphiphiles is sufficiently 
high, in order to check this possibility we carry out our third experiment. 


In the third experiment we keep the numbers of amphiphiles of both varieties fixed at 45 
but the lattice size is taken to be 10 x 10 x 100. Thus, the numbers of amphiphiles of both 
types in this experiment are identical to those in the second experiment while the initial surface 
density of the amphiphiles in this experiment is identical to that in the first experiment. In the 
third experiment we observe a vertical phase segregation (see fig.5.3) which is very similar to that 
observed earlier in the first experiment! Thus, we conclude that the vertical phase segregation 
takes place only if the initial surface density of the amphiphiles is sufficiently high. 
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Figure 5.4: Same as in fig.5.1 except that the lengths of the long and the short amphiphiles are 
17 and 13, respectively. 

5.4.1 Dependence on the lengths of the hydrophobic tails 

In this subsection we shall examine the effects of varying the lengths rig and n., of the 
hydrophobic tails by keeping the size of the heads of both types of amphiphiles fixed at = 
= 2 in a binary mixture of equal numbers of the two types of the amphiphiles. 

We have repeated the first experiment with amphiphiles of lengths 17 and 13, instead of 
17 and 7. The vertical phase segregation observed in this experiment (see fig.5.4) is qualitatively 
very similar to that observed in the first experiment. However, the phase segregation is not as 
complete as it was in fig.5.1. Similarly, again repeating the first experiment with amphiphiles of 
lengths 10 and 7 we observed vertical phase segregation of the long and short amphiphiles, as 
before (fig.5.5). Thus, in a binary mixture of equal numbers of chemically identical amphiphiles 
of two different lengths, the phase segregation is more complete when the longer amphiphiles are 
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Figure 5.5; Same as in fig.5.1 except that the lengths of the long and the short amphiphiles are 
10 and 7, respectively. 


very long and the shorter amphiphiles are very short. 


Since in this section we are considering chemically identical amphiphiles, a binary mixture of 
405 long amphiphiles and 405 short amphiphiles reduces to a system of 810 effectively identical 
amphiphiles in the limit of vanishing length difference. Even in such a situation, a fraction of the 
amphiphiles are forced out of the monolayer in the equilibrium state provided the amphiphiles are 
sufficiently long (fig.5.6). 
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Figure 5.6: The equilibrium concentration profile of the hydrophilic heads perpendicular to the 
air-water interface. The system consists of 30x30x100 lattice sites and contains 810 amphiphiles 
each of length 17. The system has been equilibrated at temperature T — 2.5. 
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5.4.2 Dependence on the ratio of the initial surface densities of 
the two components 

So far, we have taken equal initial surface densities of the long and short amphiphiles. In 
this subsection we study the effects of varying the ratio of the initial surface densities of the two 
types of amphiphiles. Since in all our earlier computer experiments we did not find any significant 
difference in the results (when appropriately normalized) for lattices of cross- section 30 x 30 
and 20 x 20, all the experiments to be discussed in this subsection have been carried out for 
20 X 20 X 100 lattices. Moreover, since the extent of unmixing is largest for the amphiphiles of 
lengths 17 and 7, the computer experiments in this subsection have been carried out with the 
lengths 17 and 7 of the long and the short amphiphiles, respectively. 

The concentration profiles of the hydrophilic heads of both the long and short amphiphiles 
are shown in figs.5.7(a)-(d) for four different ratios of the initial surface-densities of the two types 
of the amphiphiles. If the ratio Rd = N^/Ns is small (figs.5.7(b) and 5.7(d)) very few of the 
longer amphiphiles require to be sterically expelled from the long-rich layers. With the increase 
of the ratio Rd not only the population of shorter amphiphiles in the long-rich layers decrease, 
but an increasingly larger number of the longer amphiphiles are sterically expelled. 


5.4.3 Dependence on the temperature 

The values of the parameters in figs.5.8(a)-(d) are identical to those in figs.5.7(a)-(d) ex- 
cept that the temperature in fig. 5. 8 is T = 1.0 in contrast to T = 2.5 in fig.5.7. No vertical 
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Figure 5.7: The equilibrium concentration profile of the hydrophilic heads at T = 2.5 perpendic- 
ular to the air-water interface in a system consisting of 20 x 20 x 100 lattice sites. The numbers 
of long and short amphiphiles are, respectively. 300 and 60 in (a), 60 and 300 in (b), 240 and 
120 in (c), 120 and 240 in (d); the continuous line corresponds to the long amphiphiles and the 
dashed line corresponds to the short amphiphiles. 
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Figure 5.8: Same as in fig. 5. 7, except that the temperature of the system is T — 1.0. 
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Figure 5.9: Same as in fig. 5.4, except that the bending stiffness of each of the amphiphiles is 
K = 2\n contrast to K = 0 in fig.5.4. 

phase segregation takes place at the lower of the two temperatures. Clearly, temperature plays 
a crucial role in driving the out-of-plane phase segregation phenomenon. This is consistent with 
our interpretation that the phase segregation of the long and the short amphiphiles is entropy- 
driven; at low temperatures the likelihood of spontaneous buckling of the chains is very small 
and, hence, the steric interaction is suppressed. 


5.4.4 Effect of bending stiffness of the amphiphiles 

Because of the discreteness of the lattice only bends by 90« of an amphiphile at a lattice site 
is possible. Every bend of an anrphiphile is assumed to cost an energy K. This is the bending 
stiffness of the individual amphiphiles, in contrast to the bending stiffness of the amphiphilic 
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membranes considered in chapters 2 and 3. In fig.5.9 we show the equilibrium concentration 
profile for a mixture of 40-5 amphiphiles of length 17 and an *equal number of amphiphiles of 
length 13 in a 30 x 30 x 100 system, where each bending of the amphiphiles of either variety 
costs an energy of A' = 2 in the units of J. Clearly, there is no significant difference between 
the profile in fig.5.9 and the corresponding profile for A = 0 shown in fig.5.4. Thus, bending 
stiffness of the amphiphiles do not have significant effect on this phase segregation phenomenon 
unless K is very large; in the limit of A' -+ oo a total of x Ly number of amphiphiles can 
be accommodated in the monolayer at the air-water interface irrespective of the lengths of each 
of the individual amphiphiles because each amphiphile remains straight and fully extended in this 
limit. In other words, the extent of unmixing keeps decreasing as the value of K keeps increasing 
far beyond J. 


5.5 In-plane clustering of the amphiphiles 

So far we have focussed our attention on the concentration profiles of the amphiphiles in 
a direction perpendicular to the air-water interface. How are the monomers belonging to the 
long and short amphiphiles arranged in each of the planes parallel to the interface? In order 
to answer this question we have computed the distributions of the sizes of the clusters formed 
by the monomers belonging to the long and short amphiphiles in each of the lattice planes. 
A Hosen-Kopelman cluster counting algorithm [72] calculates the cluster size distribution. A 
cluster of monomers belonging to long (short) amphiphiles in a given plane is defined as the 
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set of those nearest-neighbour lattice sites which are occupied by monomers belonging to long 
(short) amphiphiles only [132]. 

We have analyzed the equilibrium state of a 30 x 30 x 100 lattice containing 405 long ' 
and 405 short amphiphiles. Let N{S) be the number of clusters of size S in a given lattice 
plane. The quantity N{S) x 5 is plotted against S in each of the figures 5.10(a)-(d); the 
continuous lines correspond to the long amphiphiles whereas the dashed lines correspond to the 
short amphiphiles. Clearly, as we move from the top of the system across long-rich layers to 
short-rich layers (and eventually deeper inside water), the nature of the clustering changes. In 
the 20th layer (fig.S.lOa) the monomers belonging to long amphiphiles form a single large cluster 
coexisting with a distribution of smaller clusters of monomers belonging to short amphiphiles. 
The size of the largest cluster formed by the monomers belonging to short amphiphiles in layer 
21 is larger than the largest cluster formed by the monomers belonging to the long amphiphiles 
(fig.S.lOb). As we move to the layer 22 we find that the clusters of the monomers of long 
amphiphiles are much smaller than in the preceding layer (fig.S.lOc). Beyond this layer even the 
size of the clusters of monomers of short amphiphiles rapidly decreases as we move deeper into 
water (fig.S.lOd). 

This behaviour is in sharp contrast with the corresponding cluster size distributions when 
the initial surface density of the amphiphiles is small. We have repeated the computation with 
90 long and 90 short amphiphiles keeping all the other parameters identical to those in fig.5.10. 
Since no phase segregation takes place in this case, almost all the monomers are confined to layers 
at and above the air-water interface. N{S) x 5 in the layers 14 to 17 are shown, respectively, in 

gs 11(a), (b), (c) and (d), the largest cluster in each of these layers consists of monomers 
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Figure 5.10; The distribution of the sizes of the clusters formed by the monomers belonging to 
the long and the short amphiphiles in the layers (a) 20, (b) 21, (c) 22 and (d) 23, respectively. 
The values of all the parameters are identical to those in fig.5.1. The continuous and the dashed 
lines correspond to the long and short amphiphiles, respectively. 
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Figure 5.11: Same as in fig.5.10, except that the number of both the long and short amphiphiles 
is 90 and figures (a), (b), (c), (d) correspond to the layers 14, 15, 16 and 17, respectively. 
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belonging to the longer amphiphiles. 

Liverpool and Bernardes [58] investigated the in-plane clustering of amphiphiles at the oil- 
water interface in a lattice model based on a Larson-type prescription for modelling the ternary 
system. An oil-water interface was constructed in the initial configuration where the lower half 
of the system was occupied by water (just as in our model) whereas the upper half was occupied 
by oil (instead of air) and the amphiphiles were uniformly distributed over the entire system. As 
the system began evolving with time, the amphiphiles moved quickly to the interface forming 
clusters. However, since the molecules of oil and water were allowed to exchange positions, the 
planar interface itself was subsequently destroyed spontaneously. On the other hand, the air-water 
interface in our model is stable as air and water are not allowed to exchange their positions. 


5.6 Conformation of amphiphiles 

Variations of the Z-components of the position of the center of mass, Zcm, of the amphiphiles 
with time (Monte carlo steps per molecule) are shown in fig. 5.12; the positions are measured 
from a reference plane in the air. The corresponding radius of gyration of the amphiphiles are 
plotted in fig. 5. 13. The figures (a) and (c) in both correspond to a high density for which vertical 
phase segregation was observed (see figs 5.1 and 5.3). On the other hand, the figures (b) and 
(d) in both correspond to a low density for which no vertical phase segregation was exhibited (see 
fig 5.2). 

In all the four figures (a)-(d). the center of mass of the long amphiphiles lie at a higher level. 
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Figure 5.12: The Z component of center of mass (Zcm) of the amphiphiles averaged over all the 
amphiphiles is plotted against MC steps; (a) Ne = N, = 405 and T = 2.5, (b) Nf = A", = 45 
and T = 2.5, (c) = 405 and T = 1.0, (d) Ni = Nf, = 45 and T — 1.0. The continuous 

line and white square correspond to the long amphiphiles and the dashed line and black square 
correspond to the short amphiphiles. 
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as expected. Moreover, comparing the figure 5.12(a) with 12(c) and 12(b) with 12(d) we find that 
at higher temperatures the equilibration is much quicker. Furthermore, at higher temperatures, 
the center of mass of the long amphiphiles lies at a significantly lower level compared to that at 
lower temperatures; this is a consequence of the "crumpling" of the amphiphiles caused by the 
increase of temperature. Most interesting observation is that there is significantly larger lowering 
of the center of mass of the short amphiphiles in fig.5.12(a) as compared to that in fig.5. 12(c); 
thermal crumpling alone is inadequate to explain such large lowering of the center of mass. This 
is consistent with our earlier observation of the vertical phase segregation and the interpretation 
in terms of entropically- induced steric repulsion. 

Since, at low densities, both the long and short amphiphiles get enough space to access all 
possible conformations, no significant "swelling" takes place by increasing the temperature from 
T = 1.0 to T = 2.5. in contrast, at high densities, the radius of gyration of the short amphiphiles 
at T = 2.5 is almost double of that at T = 1.0 (see fig.5. 13). 

In fig. 5.14 we plot the average vertical extension [< /S.Z >] of the amphiphiles as a 
function of time for different temperature and for different surface densities of amphiphiles. The 
corresponding distribution of the vertical extensions P{AZ) of the amphiphiles for different 
values of surface densities and temperatures are shown in Fig. 5.15. At low densities, there is 
practically no change in [< AZ >] (see figs. 5. 14(b) and (d)) as well as in P{AZ) of amphiphiles 
of both types when the temperature is increased from T = 1.0 to T = 2.5. At high densities, 
the data for [< AZ >] of the amphiphiles at T = 2.5 are almost identical to those observed at 
low densities. In contrast, both the long and short amphiphiles remain much more extended (i.e., 
much larger) at the same high density when the temperature is T = 1.0. 
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Figure 5.15; The probability distribution of the vertical extensions P(Aif) The parameter values 
and the notations for figs, (a), (b), (c), (d) are same as in fig.5.12. 
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Figure 5.16: Average orientational order parameter ^ for different temperatures and surface 
densities of amphiphiles. The parameter values and notations for figures (a), (b), (c), (d) are 


same as in fig.5.12. 
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Fig. 5.16 shows the variation of the average orientational order parameter with concen- 
tration and temperature. At low densities, there is practically no variation in 'i' by decreasing the 
temperature from T = 2.5 to T = 1.0. On the other hand, at higher concentrations, significant 
increase in the orientational ordering takes place by decreasing the temperature. A higher value 
of ^ indicates a higher position of the hydrophilic head compared to the location of the tip of 
the hydrophobic tail. 

The bond-orientational order parameters (f)p of the amphiphiles are plotted against the bond 
number p in fig. 5.17, where p is counted from the head. It demonstrates that, when the 
concentration is high, the bonds near the hydrophilic head are more likely to be oriented vertically 
than in any other direction; the closer is this bond to the head the stronger is this tendency. Even 
for high densities at low temperatures (fig.5. 17(c)) the last bond is equally likely to be oriented 
in any of the six allowed directions. 


5.7 Chemically different amphiphiles 

The strengths of the inter-monomer interactions were varied by varying C( and C,. When 
Ci — 10 = Cs both the long and short amphiphiles spontaneously entered into water (see 
fig.5. 18a). However, when Q = 10 and = 1 the longer amphiphiles move deeper into water 
although both the varieties move into water (fig.5.18b). On the other hand, when C, = 10 and 
Cf = 1 it is the shorter amphiphiles which move deeper into the water (fig.5. 18c). 
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Figure 5.18: The equilibrium concentration profile of the hydrophilic heads at T =- 2.5 perpendic- 
ular to the air-water interface in a system consisting of 30 x 30 x 100 lattice sites. The numbers 
of both long and short amphiphiles are 405; (a) = C, = 10, (b) Q = 10 and C., = 1, (c) 

Ci = 1 and Cs = 10. The continuous line corresponds to the long amphiphiles and the dashed 
line corresponds to the short amphiphiles. 
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Why is the majority of the short amphiphiles in the binary mixture so selectively pushed | 

out of the original monolayer? The answer to the this question is very closely related to the j 

physical mechanism that leads to the phenomenon of vertical phase separation. If the surface- ] 

,'t 

density is high, the tails of the amphiphiles encounter stronger steric repulsion. If the short I 

I 

chains move down by a few molecular layers, that leaves more empty space for the unconstrained i 

movement of the tails of the neighbouring long amphiphiles. The same would not be possible if | 

the long amphiphiles, instead of the short ones, were selectively pushed into water by the same i 

amount. It has been argued earlier in a different context [133] that the short amphiphiles, by | 

placing themselves in between long amphiphiles, create more space for the movement of the long | 

I 

amphiphiles. Entropy-driven phase separation has been observed recently in a wide variety of I 

I 

) 

systems [136], [137]. Very recently it has been demonstrated through computer simulations [135] j 

that entropy-driven phase separation is possible in polymer blends characterized by vanishing | 

Flory X where oc ‘Icy/ — f n ~ <^ 22 . being the energy of interaction between two nearest- | 

1 

neighbour segments of type i and j. | 

I 

i 

The system studied in this chapter have some similarities with binary mixtures of polymer j 

f 

I 

chains grafted on solid surfaces [140, 141]. A direct comparison of such polymer brushes with our 
results would require a study of mixture of grafted polymers of two drastically different lengths 
when the surface coverage is very high and temperature is not too low. | 
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5.9 Summary and conclusion 
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In this chapter we have investigated the cnt-of-plane phase segregation of the amphiphiles 
binary mixture at the air-water interface. We have also studied the nature of the in-plane 
clustering of the monomers belonging to the long and the short amphiphiles. The entropy-driven 
phase segregation can be observed in laboratory experiments with chemically identical amphiphiles 
of two different lengths. For complete unmixing of the two components the long amphiphiles 
should be as long as possible and the difference in the lengths of the two types of the amphiphiles 


should also be large. 



Chapter 6 


Self assemblies of single chain and 
double chain amphiphiles 


6.1 Introduction 

All the previous chapters deal almost exclusively with the interfacial properties except in 
chapter 5 where we have studied both the interfacial as well as bulk properties. In the remaining 
part of the thesis we will study the bulk properties of the self assembling amphiphilic systems. 
In chapters 4 and 5 only single-chain amphiphiles have been considered although the general 
conclusions drawn from our MC studies are valid also for double-chain amphiphiles. In this 
chapter we shall study the bulk properties of the single chain as well as double chain amphiphiles. 
This chapter bridges the gap between chapter 7 where we shall develop a model for gemini 
surfactants which are double chain amphiphiles and chapters 4 and 5 where we have considered 
only single chain amphiphiles. 
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In this chapter we compare and contrast the self-assemblies of single-chain and double-chain 
amphiphiles, respectively, in water. We have studied binary mixture of amphiphiles at the air- 
jvater interface in chapter 5. Here we study the self-assemblies of binary mixtures of single-chain 
and double-chain amphiphiles and compare these with those formed by the corresponding single 
species of amphiphiles 


6.2 Micelle formation 

For single chain amphiphiles micellar aggregates form at low concentration of amphiphiles, 
i.e. at concentration just above CMC. Stauffer, Jan and Pandey [45], using type studied 

the behaviour of amphiphiles in water beginning from a disordered distribution of chains and 
found that they always formed clusters (which they called micelles) whose average size, denoted 
by a peak in the size distribution, increased with time. Stauffer et al. [46] performed a more 
detailed study of the properties of the model for structure HiHaTa. The length of the chains 
and the strength of the hydrophilic and hydrophobic interaction were varied. They used Glauber 
Kinetics and only the reptation movement. Simulations were performed on lattices of sizes up 
to 173 X 173 X 173. For the same type amphiphile as that studied in [45] at a concentration 
of (pa = 0.1, they observed that between temperature of t = 2.1 and t = 4.1, a distribution of 
micelles with a stable mean size (weight average) was formed. Below t = 2.1, the micelles grow 
indefinitely and above t = 4.1 a mixed phase was observed. A simulation of the T-LiMo with the 
spin of the head 5 = -j-2 showed the growth of the mean cluster size to be slowed down^ while 
one of UiN^Ti with head S = ^2 and tail S = -2 gave similar behaviour as the with 
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head S' = +1 and tail S = -1. 

Jan, Stauffer and collaborators [46, 47] extended the model further to describe single-chain 
surfactants with ionic heads. According to their formulation, the monomers belonging to the 
ionic heads have Ising spin 4-2 to mimic the presence of electric charge. The repulsive interaction 
between a pair of ionic heads is taken into account through an (antiferromagnetic) interaction 
J = -1 between pairs of nearest neighbour sites, both of which carry spins +2; however, the 
interaction between ail other pairs of nearest-neighbour spins is assumed to be J = 1. By 
restricting the range of the repulsive (antiferromagnetic) interaction between the "charged” heads 
to only one lattice spacing one is effectively, assuming very strong screening of the Coulomb 
repulsion between ionic heads by the counterions. 

When head-head repulsion was included, the mean micelle size remained constant more or 
less independent of temperature. The repulsion was thus conjectured to be important for the 
stability of micelles [46]. Jan and Stauffer [47] also studied surfactants with structure 

4-2, -fl, 0, — 1, — 1, — 1 also with head-head repulsion and found stable micelles. However, in 
any of the above mentioned work the morphology of the micellar aggregates was not studied. 
Smit et al. [60] showed the formation of spherical micelles with Molecular dynamics simulation. 
Here we report the formation of stable spherical micelles by looking at the actual morphology of 
the aggregates formed during Monte Carlo simulation; to my knowledge, this is the first direct 
evidence of the formation of spherical micelles in the Larson model of single-chain surfactants. 

We have performed simulations for much longer amphiphiles in water. We allowed 

all the moves for amphiphiles described in chapter 4 :(i) reptation, (ii) kink, (iii) spontaneous 
buckling and (iv) pull. The simulations were performed on lattice of size 100 x 100 x 100 at 
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in amphiphile concentration Oa — 0.009 and temperature t — 2.2. Starting with amphiphiles 
jniformly distributed in the cubic lattice, spherical micelles formed at the end of the simulation 
see figure 6.1). 

Simulation with amphiphiles also gives spherical micelles; figure 6.2 shows the 

snapshot of the 3-dimensional micelles. 

Next we performed simulation for a binary mixture of single chain amphiphiles in water. Since 
A/e are studying a binary mixture of the amphiphiles, let us represent them by the two symbols 
Hg^J^iTmi and 'Hq^N'iTm, for the long and short amphiphiles, respectively, so that the lengths of 
the long amphiphiles is Le = qc + mi + l;and that of the short amphiphiles is Lg = g., 4- -f 1. 
Ml our investigations have been carried out for a fixed value of qe = Qs = 2, Lc — 12 and Lg = 6. 
The ratio of the volume fraction of the long and short amphiphils was kept at 7 : 3. In this case 
also we get a distribution of spherical micelles. Figure 6.3 shows one snapshot of such micellar 
aggregates. 


6.3 Vesicle formation 

In water, amphiphiles can form "bilayers” which consist of two monolayers of oppositely 
Driented amphiphiles such that the hydrophobic tails of amphiphiles in the two monolayers are 
bidden inside thereby minimizing contact with water. A bilayer can close itself to avoid contact 
aetween the hydrocarbon tails and water at the edges of the membrane and form what is called 
vesicles. Thus vesicles separates an inner region of an aqueous fluid from another aqueous fluid 
n its exterior. In recent years there has been an increasing interest in vesicles mainly because of 
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Figure 6.2; Same as in figure ().l(.)) cxi epf tli.it mav in !) 

their wide application in biology and medicine as model tell tiiemlH.utes, .is well as their strong 
potential as drug carriers and others encapsulating agents of iiuliistti.il televatue [143]. However, 
whereas the vesicles are widespread and form spontaneously in nature, the artificial formation of 
vesicles requires the input of some form of energy, for example sonifit at ion. .Safran et al-[142] have 
shown that mixture of two surfactants can lead to spontaneous vesicles formation. In contrast 
to the situation for single amphiphiles, where large vesicles are usu.illy not ener|',ctically stable 
in comparison with fiat bilayers, vesicles composed of two spiu ies of amphiphik'S can have lower 
curvature energies than flat bilayers and hence more stable. Vesi< les h.ive been found to form 
spontaneously in some aqueous surfactant systems, iriduding solutions (.ontaining (i) mixtures of 
lecithin and lysolecithin [144], (ii) mixtures of long and short c hain lecithins [14.6], (iii) mixtures 
of cataionic and anionic surfactants [90]. 
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Recently a microscopic lattice model of double-chain amphiphiles with a single head was 
developed by Bernardes by extending the Larson model of single-chain amphiphiles in ternary 
microemulsions and its later generalizations. A lattice model of double-chain amphiphiles with 
a single hydrophilic head can be described by the symbol TmJ^p'HqMpTm- rn, p and q are 
integers denoting the lengths of the tail, neutral region, and head, respectively, in the units 
of lattice sites. Extensive MC simulation showed the spontaneous formation of vesicles. The 
simulations were performed on amphiphile of length £ = 12 in a binary solution. 

These correspond to double chain amphiphiles which have been conjectured to be more suited 
to forming biiayer membranes [33]. Simulations were performed on a 40 x 40 x 50 system at 
a temperature of t = 2.2 with 533 amphiphiles for up to 4 x 10^ MC steps. Beginning with 
a uniform distribution of amphiphiles throughout the whole lattice, a large vesicle is formed at 
the end of the simulation. It should be noted that at this temperature the system is below 
its condensation transition temperature and consequently all the amphiphiles tend to aggregate 
in one large cluster. Simulation with shorter double chain amphiphiles ^ = 10 has also been 
performed yielding vesicle formation. 

We extend Bernardes model to include binary mixture of double chain amphiphiles and report 
the spontaneous formation of vesicles of different sizes. The simulations were performed on a 
mixture of long amphiphiles TaUiU^MiTa of length Lt = 12 and short amphiphiles 
of length Ls = 6. The ratio of the long and short amphiphiles was 7 : 3, Initially both types of 
amphiphiles were uniformly distributed in the cubic lattice of size LxLxL^. Periodic boundary 
conditions were adopted in the x and y directions. Fixed boundary conditions were applied in the 
z direction where the states of the two uppermost and two lowermost planes were kept "frozen” 
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Figure 6.4: (a) Snapshots of the vesicles formed by the binary mixture of double chain am- 
phiphiles. Black sphere represent head and grey sphere represent tail. 


in the corresponding initial states, and no movement of the amphiphiles through these planes was 
allowed. Amphiphiles can move in the lattice changing position with water molecules through 
different types of movements discussed in chapter 4; namely reptation, kink and buckling. After 
10^ MC steps we get vesicles of different sizes. Figure 6.4(a) shows different sizes of the vesicles 
formed in our simulation. For comparison, the vesicles formed by the single component double 
chain ampphiphiles has been shown in Figure 6i|(b). 
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Figure 6.4: (b) Snapshots of the vesicles formed by the single component double chain am- 
phiphiles. Black sphere represent head and grey sphere represent tail. 
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Figure 6.4: (c) A typical cross-section of the vesicle Black sphere represent 


head and grey sphere 


represent tail. 
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6.4 Conclusion 

In this chapter we have compared and contrasted the self assemblies of single chain and 
double chain amphiphiles. In binary mixture of water and single chain amphiphiles^ when the 
concentration and temperature is sufficiently low^ we get spherical micelle-like aggregates. For a 
mixture of long and short single chain amphiphiles in water also we get spherical micelles. For 
double-chain amphiphiles when only one species are present in wate^we get one large single 
vesicle. On the other hand, the amphiphiles form simultaneously more than one vesicle (of 
relatively smaller sizes) when a fraction of the long amphiphiles are replaced by shorter double- 


chain amphiphiles. 



Chapter 7 


A Microscopic Model of Gemini 
Surfactants: Self-assemblies in Water 
and at Air- Water Interface’ 

7.1 Introduction 

It IS well known that micelles are formed when the concentration of the surfactants in water 
exceeds what is known as critical micellar concentration (CMC) [10]. In reality. CMC is not a 
single concentration but a very narrow range of concentration (it is more appropriate to call it 
characteristic micellar concentration [46, 47]). Experimentally, for example, the CMC is found 

to (i) increase with the increase in the ionic strength of the heads in case of Ionic surfactants, 

^The following papers are based on this chapter: 

^ ^eSon! tCond-mat/9706052. submitted for publication , P. K. Maiti and D. Chowdhury, 
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(ii) decrease with increasing tail length [1], (iii) decrease with decreasing temperature [22]. (iv) 
increase with the increase in the size of head in case of non-ionic surfactants [54], The CMC 
calculated from computer simulation with the Larson-type models excellently agree with the above 
experimental results [46, 47, 54, 52], 

In this chapter we develop a microscopic model of ^emmz surfactants and, carrying out Monte 
Carlo (MC) computer simulations, investigate how the shapes and sizes of these molecules as well 
as their mutual interactions and their interactions with the molecules of water give rise to some 
unusual aggregation phenomena. We simulate gemini surfactants with (a) hydrophobic spacers 
and also those with (b) hydrophilic spacers. For each of these two types of gemini surfactants 
we consider both ionic and non-ionic (but polar) hydrophilic heads. Another aim of this chapter 
is to report the results of a complimentary MC study of the spatial organization of these model 
gemini surfactants (particularly, their tails and spacers) at the air-water interface in order to 
answer some of the fundamental questions raised on this point and speculations made in the 
literature [146]. A microscopic model for single-chain surfactants at the air-water interface was 
developed in chapter 5. In this chapter we replace the single-chain surfactants by the model gemini 
surfactants developed here, thereby getting the desired microscopic model of gemini surfactants 
at the air-water interface. 

In contrast, two unusual features of the CMC of gemini surfactants with ionic heads and 
hydrophobic spacer are: 

(i) for a given fixed length of each of the two tails,' the CMC increases with the length of the 
spacer till it reaches a maximum beyond which CMC decreases wth further increase of the spacer 
length [15, 152-154]; 



A Microscopic Model of Gemini SurfBctsnts: Self-sssemblies in Water and at Air-Water 
jQ Interface 

i) for a given length of the spacer, the CMC increases with increasing tail length [12, 13]. 
loreover, the micellar aggregates formed by the gemini surfactants with short spacers even at low 

oncentrations just above the CMC are " long, thread-like and entangled [16, 155], in contrast to 

ne spherical shapes of the micelles formed by single-chain surfactants at such low concentrations, 
urthermore, the CMC of gemini surfactants with ionic heads and hydrophilic spacer decrease 
nonotonically with the increase of the length of the spacer [156]. Our aim is to understand 
he physical origin of these unusual properties of gemini surfactants. We also make some new 
)redictions on the morphology of the micellar aggregates of gemini surfactants with long tails 
ind long spacer. 

Synthesis of gemini surfactants with non-ionic (polar) heads in laboratory experiments re- 
nains one of the challenging open problems. But, our computer experiments on model gemini 
lurfactants with non-ionic heads enable us to predict their morphologies and the variation of their 
ZMC with the lengths of tails and spacers. 

In the presence of the air-water interface, where do the tails and the spacer of an isolated 
jemini surfactant find themselves- do they lie inside water or outside (i.e., in the air), do they 
jet crowded close to the interface or do they spread out as far away from the interface as 
30ssible? How does the effective area of cross-section of an isolated gemini surfactant at the 
air-water interface vary with the increase of the length of the spacer when the spacer is (a) 
lydrophobic, (b) hydrophilic? How do the conformations of the gemini surfactants and spatial 
organizations of their tails and spacers vary with the increase of the density which gives rise to 
unavoidable interactions (both direct and entropic) among the surfactants? We try to answer 
these fundamental questions by carrying out computer experiments on a microscopic model that 
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we propose here. 

The model and the characteristic quantities of interest are defined in section 2. The results 
on the micellar aggregates, in bulk water, formed by the gemini surfactants with hydrophobic and 
hydrophilic spacers are reported in two subsections of section 3. The results of the investigations 
on the spatial organization of the gemini surfactants with hydrophobic and hydrophilic spacers 
at the air-water interface are given in section 4. Finally, a summary of the main results and the 
conclusions drawn from these are given in section 5. 


7.2 The Model and The Characteristic Quantities of 
Interest 

7.2.1 General aspects 

Now we propose a microscopic lattice model of gemini surfactants. In terms of the symbols 
used earlier to denote the primary "structure" of the microscopic lattice model of single-chain 
surfactants, Bernardes’ lattice model of double-chain surfactants, with a single hydrophilic head, 
can be described by the symbol Tm.J^p'K-q-^pTm- In terms of the same symbols, the microscopic 
lattice model of a gemini surfactant, which we propose here, can be represented by the symbol 
where n is the number of lattice sites constituting the spacer represented 
by the symbol <S (see fig.7.1). 

In our model, an Ising spin variable S is assigned to each lattice site; Sj = 1 (—I) the 
i-th lattice site is occupied by a water (oil) molecule. If the j-th site is occupied by a monomer 
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Figure 7.1: Larson-type models of single-chain, double-chain and gemini surfactants. 
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belonging to a surfactant then Sj = depending on whether the monomer at the jth 

site belongs to head, tail or neutral part. A monomer belonging the spacer is assigned a spin 
value -1 if the spacer is hydrophobic and +1 if the spacer is hydrophilic. The monomer-monomer 
interactions are taken into account through the interaction between the corresponding pair of Ising 
spins which is assumed to be non-zero provided the spins are located on the nearest-neighbour 
sites on the lattice. Thus, the Hamiltonian for the system is given by the standard form 


H = -J Y. (7.1) 

<ij> 

where attractive interaction (analogue of the ferromagnetic interaction in Ising magnets) corre- 
sponds to .7 > 0 and repulsive interaction (analogue of antiferromagnetic interaction) corresponds 
to J < 0 [46, 47]. The temperature T of the system is measured in the units of J/ks where ks 
is the Boltzmann constant. 

To our knowledge, all the gemini surfactants synthesized so far have ionic heads. Therefore, 
we incorporate the effects of the ionic heads following the same way as has been described in 
chapter 6. The monomers belonging to the ionic heads have Ising spin +2 to mimic the presence 
of electric charge. The repulsive interaction between a pair of ionic heads is taken into account 
through an (antiferromagnetic) interaction (J = —1). 

In order to predict the properties of gemini surfactants with non-ionic (polar) heads and to 
investigate which of the aggregation phenomena exhibited by the ionic gemini surfactants arise 
from the electric charge on their ionic heads, we have also considered a model of gemini surfactants 
with non-ionic polar heads which is obtained from the model of ionic gemini surfactants by 
replacing all the -f-2 ising spin variables by Ising spin +1 (and, accordingly, the interactions 
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Figure 7.2: Schematic representation of the moves attempted by every surfactant at every MC 
step. 

— 1 between the heads on nearest-neighbour sites are replaced by H-l). Moreover, in order to 
investigate the role of the chain stiffness we have used a chain bending energy [81]; every bend 
of a tail or a spacer, by a right angle at a lattice site, is assumed to cost an extra amount of 
energy K{> 0). 

Starting from an initial state (which will be described in the subsections 7.2.2 and 7.2.3). the 
system is allowed to evolve following the standard Metropolis algorithm: each of the attempts 
to move a surfactant takes place certainly if AE < 0 and with a probability proportional to 
exp{—AB/T) if AE > 0. where AE is the change in energy that would be caused by the 
proposed move of the surfactant under consideration. 
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Next, we specify the allowed moves of the surfactants for the appropriate sampling of the 
states of the system in a MC simulation. We have generalized the moves for single chain surfac- 
tants (discussed in chapter 4) appropriately for the gemini surfactants and schematic represen- 
tation of these moves are shown in fig.7.2: although the moves are illustrated using one of two 
tails, each of the moves (except reptation, which involves the entire surfactant) in our algorithm 
is equally likely to be applied on the two tails as well as the spacer. The moves allowed for the 
surfactants in our model are as follows: 

(i) reptation: this is identical to the reptation move for single-chain surfactants described earlier; 

(ii) spontaneous chain buckling: a portion in the middle of one of the two tails or the spacer 
is randomly picked up and allowed to buckle with the probability mentioned above; (iii) kink 
movement: a kink formed by the buckling or reptation is allowed to move to a new position 
with the appropriate probability calculated according to the prescription mentioned above; (iv) 
pull move.: this is the reverse of spontaneous chain buckling; a buckled part of one of the two 
tails or the spacer is pulled so as to make it more extended. In addition to these moves, each of 
the surfactants is allowed to move laterally in one of the six possible directions, which is chosen 
randomly with probability 1/6 and each monomer of the surfactant is moved in that direction by 
one lattice spacing. Each of these moves is possible only if the new positions of all the monomers 
are not occupied simultaneously by monomers belonging to other surfactants. Each surfactant is 
allowed to try each of the above mentioned moves once during each MC step. The moves of the 
molecules of water and air are described in the subsections 7.2.2 and 7.2.3 below. 

In principle, one can study the aggregation of gemini surfactants deep inside bulk water 
and their spatial organization close to air-water interface by MC simulation of a single system 
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ere the lower part of the lattice representing water is sufficiently large and the density of 
'factants is also sufficiently high so that a large fraction of these can be found deep inside 
iter in the form of micellar aggregates. However, for the convenience of computation, we study 
ese two aspects of the problem separately: in the first part we investigate only the phenomenon 
aggregation of gemini surfactants in bulk water and in the second part we investigate only 
e spatial organization of the tails and spacers of the gemini surfactants at (and near) the air- 
ater interface when the total volume fraction occupied by the surfactants is quite small. Both 
lese parts of our computations are based on the general model described above and the specific 
fferences involved in these two are explained separately in the next two subsections. 


.2.2 Model of Gemini Surfactants in Bulk Water 


In order to investigate the spontaneous formation of micellar aggregates and their morphol- 
gy, model gemini surfactants are initially dispersed randomly in a x Ly x system which 
ontains only water and surfactants. Periodic boundary condition is applied in all the three di- 
jctions thereby mimicking bulk water which is infinite in all the directions. So far as the moves 
f the molecules of water are concerned, each molecule of water is allowed to exchange position 
nth a monomer belonging to a surfactant, provided that is necessary for the implementation of 
n attempted move of a surfactant. 
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7.2.3 Model of Gemini Surfactants at Air- Water Interface 

Just as in the preceding subsection, the system under investigation is modeled as a simple 
cubic lattice of size x Ly x T,. However, in contrast to simulating a bulk of water, which 
is infinite in all directions, we now simulate a semi-infinite vertical column of air separated from 
a semi-infinite vertical column of water below by a sharp horizontal air-water interface. In the 
Cartesian coordinate system the horizontal air-water interface is parallel to the XT-plane and the 
vertically dowmtHmi direction is chosen as the -l-T-axis. Each of the molecules of water and air 
can occupy a single lattice site. A classical Ising spin variable 5 is assigned to each lattice site; 
Si = 1 ( — I ) if the /-th lattice site is occupied by a water molecule (air molecule or empty). Our 
prescription for assigning the Ising spin variables to the sites occupied by the monomers of the 
amphiphiles is identical to that given in the preceding subsection. Periodic boundary conditions 
are applied along the A' and V' directions. The lattice sites in the uppermost and the lowermost 
layers are occupied by "down" and "up” spins respectively, which were not updated during the 
computer simulation. These boundary conditions mimic the physical situation mentioned above 
which we intend to simulate. 

In the initial state the surfactants are so arranged that their spacers lie flat, and fully ex- 
tended, horizontally in the first layer of water immediately below the air-water interface and their 
tails are fully extended vertically into air. The system is then allowed to evolve towards equilib- 
rium following the Metropolis algorithm explained earlier. So far as the moves of the individual 
molecules are concerned, air and water are not allowed to exchange positions, as dispersion of 
air and water inside each other is not possible in our model. However, if some monomers of a 
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ictant come out of water the vacant sites are occupied by inserting water molecules; this is 
iistent with our assumption that the water column is semi-infinite in the ^-direction. More- 
, we impose an additional constraint that none of the heads of the surfactant molecules can 
le out of water. 

5.4 Characteristic Quantities of Interest 

The most direct approach to investigate the morphology of the micellar aggregates and the 
tial organization of the different parts of individual surfactants is to look at the snapshots of 
system after equilibration. For studying the variations of CMC with the lengths of the tails 
I spacers one has to use a well-defined prescription for computing the CMC; this is a subtle 
nt as the CMC is not a unique single concentration, as mentioned before. We follow the 
scription proposed, and used successfully in the case of single-chain surfactants [46, 47]; we 
ntify CMC as the amphiphile concentration where half of the surfactants are in the form of 
ated chains and the other half in the form of clusters consisting of more than one neighbouring 
phiphiie. 

We have introduced a quantitative measure of the effective cross-sectional area A of the 
nini surfactants projected onto the air-water interface. We compute and |A?/|m which 

the maximum differences in the X- and 5^-coordinates, respectively, of the monomers and 
fne A as 

A = [(|Aa:U)2 + (lAyP^jV^ (7.2) 


The vertical extension < Z > is defined as the difference in the Z- coordinates of the 
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highest and lowest monomers (i.e. monomers with highest and lowest value of ^-coordinates) 
of a single surfactant. A quantitative measure of the gross features of the spatial organization 
of the tails and spacers of the gemini surfactants at the air-water interface is the equilibrium 
profiles of the concentrations of the corresponding monomers in the ^-direction, i.e., in the 
direction perpendicular to the air-water interface. More precisely, at each molecular layer, we 
count separately the number of monomers belonging to the tails and the spacers (and also the 
heads and neutral parts) in that particular layer and average the data over sufficiently large 
number of configurations after equilibration of the system. 


We have carried out MC simulations of the model TmMp'HgSn'HqM^pTm of gemini surfactants 
for p = q = I and for three different values of the tail length, namely, m = 5, 15 and 25. In our 
simulation of the surfactants at the air-water interface we do not find any observable difference 
in the concentration profiles obtained in single runs for 100 x 100 x 100 systems and for larger 
systems containing identical surface-density of surfactants, all the profiles reported in this chapter 
have been generated for system sizes Lx = Ly — = L = 100 by averaging over sufficiently 

large (10-25) number of runs. The same size of the system was also found to be large enough to 
avoid severe finite-size effects on the CMC data; each of the data points for the CMC is obtained 
by averaging over typically 10 runs. For a given m we have computed the CMC for spacer lengths 


2 < n < 20. 
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3 Micellar Aggregates of Gemini Surfactants 

3.1 Aggregates of Gemini Surfactants: Results for Hydropho- 
bic Spacers 

The CMC of ionic gemini surfactants with hydrophobic spacers are plotted against the 
acer length for two different lengths of the tail, namely, m = 5 and m = 15, in figs. 7.3 
) and 7.3(b), respectively. The non-monotonic variation of CMC with the spacer length, in 
;s. 7.3(a) and 7.3(b), is in qualitative agreement with the experimental observations [152-155]. 
awever, this is in sharp contrast to the monotonic decrease of the CMC with the length of the 
'drophobic tail of single-chain model surfactants of the type TmMpUq [46, 47, 54]. Moreover, 
r a given length of the hydrophobic spacer, the CMC of this type of gemini surfactants increases 
hen the bending stiffness K of the hydrophobic chains is switched on (see figs. 7.3(a) and (b)). 
jtthermore, we observe that, for a given length of the hydrophobic spacer, the CMC of ionic 
amini surfactants increase with the increase of the tail length (see fig.7.4); this trend of variation 
also consistent with the corresponding experimental observations [12, 13]. 

For a given tail length (see fig.7.5 for m = 5 and m = 15), the CMC of model gemini 
Lirfactants with hydrophobic spacers decreases monotonically with the increase in the spacer 
;ngth when the polar head group is non-ionic. This is in sharp contrast to the non-monotonic 
ariation observed for ionic gemini surfactants. However, for a given spacer length, the trend 
)f the variation of CMC of non-ionic gemini surfactants with the tail length is similar to that 
)bserved for ionic gemini surfactants, i.e., CMC increases with the increase of the length of the 
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Figure 7.3: (a)Variation of CMC of ionic geminis with hydrophobic spacer length; m — 15, 
T = 2.2. The symbols □ and x correspond to iC = 0 and iC = 2, respectively, (b) Same as (a), 
except that m = 5. The symbols A and * correspond to iC = 0 and K = 2, respectively. The 


continuous curves are merely guides to the eye. 
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igure 7.4: Variation of CMC of ionic geminis with tail length at T = 2.2 for three different 
:ngths of the hydrophobic spacer, namely, n = 6.8,16. The straight lines connecting the 
uccessive data points are merely guides to the eye. 
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Figure 7.5: Variation of CMC of non-ionic geminis with the length of the hydrophobic spacer; 
m = 15 (□) and m = o (A) both with K = 0 and at T = 2.2. The continuous curves are 


merely guides to the eye. 




Figure /.b; Snapshots of the micellar aggregates formed by ionic geminis with hydrophobic 
spacer; m - 15, n 2 and K - 0 at T = 2.2 when the surfactant density is 0.007. The symbols 
black spheres, dark grey spheres and light grey spheres represent monomers belonging to head, 
tail and spacer, respectively. 

tail. 

Instantaneous snapshots of the micellar aggregates formed by long tailed {m = 15) gemini 
surfactants with ionic heads and hydrophobic spacer are shown for spacer lengths n = 2 (fig. 
7.6) and n = 16 (fig. 7. 7). The morphology of the aggregates in fig.7.6 are similar to the 
long, thread-like and entangled” micelles observed in laboratory experiments [16] and in MD 
simulations [148] on gemini surfactants with short hydrophobic spacers. Moreover, our data in 
fig-7.7 suggest that rod-like ("columnar”) micelles are formed by gemini surfactants with ionic 
head and long tail (rti = 15) when the length of the hydrophobic spacer Is also long (n = 16). 
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Figure 7.7: Same as in fig.6, except that n = 16 and the density is 0.005. 

'he morphologies of the aggregates in fig.7.6 and 7.7 are in sharp contrast with the spherical 
hape of the micelles (see fig.7.8) formed by single-chain ionic surfactants of comparable tail size 
ven at concentrations somewhat higher than those in the figures 7.6 and 7.7. 

There is no significant difference between the morphologies of the micellar aggregates of 
onic and non-ionic single-chain model surfactants represented by the symbol TmMpHq [98]. 
iimilarly, we do not observe any significant difference also in the shapes of the aggregates of 
Dnic surfactants (figs.7.6 and 7.7) and those of non-ionic gemini surfactants (see figs. 7. 9 and 
.10) with hydrophobic spacers, for given values of m, n and comparable concentration, in spite 
if qualitatively different trends of variation of their CMCs with spacer lengths. 

Note that the rod-like micelles, shown in fig.7.7, are formed when the length of the hydropho- 
be spacer and the combined length of the tail and the neutral part of the gemini surfactants are 




Figure 7.8: Snapshots of micellar aggregates formed by single-chain ionic surfactants with m 14 
and the density 0.0 L. The symbols black spheres and grey spheres represent monomers belonging 


to head and tail, respectively. 
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figure 7.12: Variation of CMC of ionic geminis with hydrophilic spacer length; m. = 15, 

= 2.2. The symbols □ and x correspond to iv = 0 and K = 2, respectively. 

)Oth equal to 16. Does this imply that rod-like micelles are formed whenever the hydrophobic 
pacer and the tail are equal (or comparable) in length? In order to answer this question we have 
Iso looked at the snapshots of the micellar aggregates of similar gemini surfactants with shorter 
ails and spacers; a typical example, shown in fig.7.11, corresponds to rn = 5, n = 6. The fact 
hat these micelles are also "long, thread-like and entangled”, like those in fig. 7. 6, in contrast 
0 the rod-like micelles of fig.7.7, suggests that the morphology of the ionic gemini surfactants 
nth hydrophobic spacers is dominantly determined by the length of the spacer; long, thread-like 
licelles are formed if the spacer is short and rod-like micelles are formed if the spacer is long. 

.3.2 Aggregates of Gemini Surfactants; Results for Hydrophilic 


Spacers: 



7.3 Micellar Aggregates of Gemini Surfactants 


149 



f'igurr! 7.1.3; Same as fig. 7.6 except that the spacer is hydrophilic. 

In fig. 7. 12 we plot the CMC against the length of the hydrophilic spacer for gemini surfac- 
tants with io'nic head and tail length m = 15 (the qualitative features of the corresponding curve 
for m = 5 are very similar and, therefore, not shown). In contrast to the non-monotonic variation 
of CMC observed earlier with the variation of the length of hydrophobic spacers, now we find 
a monotonic decrease of CMC with the increase of the length of the hydrophilic spacer. This 
trend of variation is in qualitative agreement with the corresponding experimental observation 
[156]. Moreover, for given lengths of the hydrophobic spacer and the tail, n and m, the CMC 
for the bending energy K = 2 is lower than that for iC = 0 (see fig.7.12); this trend of variation 
is exactly opposite to the corresponding trend observed earlier in the case of gemini surfactants 
with hydrophobic spacers. 

The snapshots of the micellar aggregates formed by the gemini surfactants with ionic heads 
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Figure 7.14: Same as fig. 7.7 except that the spacer is hydrophilic. 

and hydrophilic spacer are shown for spacer lengths n = 2 (fig. 7.13) and for ;= Ki (fig.7.14) 
for densities which are identical to those in the figs.7.6 and 7.7, respectively. Comparing the 
morphologies of the aggregates in fig.7.6 and fig.7.13 we find that the gemini surfactants with 
hydrophobic spacers form coarser (albeit fewer in number) aggregates compared to the corre- 
sponding geminis with hydrophilic spacers; this is also consistent with the fact that the CMC of 
the gemini surfactants with spacer length n = 2 are higher when the spacers are hydrophilic as 
compared to that for hydrophobic spacers. 

The difference in the morphologies of ionic geminis with hydrophobic and hydrophilic spacers 
s much more striking when the spacer is longer (n = 16) (compare the fig.7.7 with fig.7.14)- 
:he micelles are more or less spherical when the spacers are hydrophilic! 

An important difference between the micellar aggregates of gemini surfactants with hy- 
drophobic spacers and those with hydrophilic spacers is that more spacer monomers are found on 
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the outer surface of the aggregate (i.e., in contact with water) when the spacer is hydrophilic. 
This is consistent with one’s intuitive expectation because the hydrophilic spacers like to be in 
contact with water. 


The snapshots of the micellar aggregates of non-ionic gemini surfactants with hydrophilic 
spacers are very similar to those for the corresponding ionic gemini surfactants (and, therefore, 
not shown in any figure). 

Hydrophilic spacers gain energy by remaining surrounded by water. On the other hand, 
hydrophobic spacers as well as tails try to avoid contact with water by hiding inside micellar 
aggregates. That is why in the snapshots of micellar aggregates we see that a larger number of 
monomers belonging to the spacers are in contact with water, when the spacers are hydrophilic, 
than those when the spacers are hydrophobic. And this is prominent particularly for long spacers. 


7.4 Spatial Organization of Gemini Surfactants at Air- 
Water Interface: Results for Hydrophobic and 
Hydrophilic Spacers 

7.4.1 Dilute regime 

First let us consider the dilute regime where the concentration of the surfactants is so low 
that not only all of them remain, almost certainly, at the air-water interface but every surfactant 
may be regarded as, effectively, isolated from each other. In this limit the cross-sectional area A 
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Figure 7.15: Variation of cross-sectional area of isolated individual gemini surfactant with spacer 
spacer length. The solid line for hydrophobic spacer and broken line for hydrophilic spacer. To 
give an indication of the accuracy of the data points the error bar of only one point has been 


shown. 




Spt;" Hyd.p.„.c 

— — ^ 

of the molecules is determined by only intra-molecular interactions, which is dominated by the 

steric (entropic) interactions among the tails and the spacer. We plot the cross-sectional area .4 

of isolated individual gemini surfactants as a function of the length of the spacer in fig.7.15, in 

both the cases of (a) hydrophobic and (b) hydrophilic spacers. The spacer is very stiff when its 

length is n = 2 as no wiggle can form. The area ^ for n = 4 is smaller than that forn = 2 

irrespective of the nature of the spacer (i.e., hydrophobic or hydrophilic); this is caused by the 

formation of wiggle on the spacer which brings the two heads closer. Further increase of the 

spacer length gives rise to a linear increase of the area A. However, a sharper increase in A takes 

place when the length of the spacer becomes equal to that of the tails; on both sides of this 

regime of sharp rise, the rate of increase of A with n is practically identical. 

Because of its stiffness against wiggle formation, the spacer of length n = 2 can buckle 
neither towards air nor towards water and remains parallel (like a rigid rod) to the air-water 
interface. Therefore, if n = 2, the cross-sectional area A of isolated gemini surfactants with 
hydrophilic spacers is practically identical to that of gemini surfactants with hydrophobic spacers. 
However, for all larger values of n, A is smaller if the spacer is hydrophilic; a hydrophilic spacer 
buckles into water thereby leaving most of the space in the air above the heads available for 
occupation by the tails. On the other hand, the hydrophobic tails take up a substantial amount 
of available space in a cap-like volume in the air just above the heads thereby forcing the tails to 
spread out radially outward and, hence, increasing the effective area A. 


Evidence in support of this scenario emerges also from the plots of vertical extension < Z > 
of the isolated gemini surfactants against the length of the spacer (see fig.7.16); a larger < Z > 
of gemini surfactants with hydrophilic spacer, as compared to those of gemini surfactants with 
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Figure 7.16: Variation of vertical extension of individual gemini surfactant with spacer length. 
The solid line for hydrophobic spacer and broken line for hydrophilic spacer. 




7.4 Spatial Organization of Gemini Surfactants at Air-Water lni-Prf:,rfl 
and Hydrophilic Spacers: 

hydrophobic spacer of identical length, arises from the fact that the 


: Results for Hydrophobic 
155 

hydrophilic spacers buckle 


into water while their tails remain outside water. 


7.4.2 High Surface density regime 

The concentration profiles of the tails of the gemini surfactants with hydrophobic spacers 

are shown in fig.7.17(a) and the corresponding concentration profiles for gemini surfactants with 

hydrophilic spacers are shown in fig.7.17(b). In the case of gemini surfactants with hydrophobic 

spacer, the spacers minimize contact with water by arranging themselves just outside the water, 

but do not venture out too far from the interface. On the other hand, the hydrophilic spacers 

gain energy by moving inside water thereby leaving more space just outside water which become 

available for occupation by the tails; consequently, one would have naively expected, the tails of 

the gemini surfactants with hydrophilic spacers are likely to be found closer to the interface than 

those of the gemini surfactants with hydrophobic spacers. However, what we observe in reality 

in fig.7. 17(b) is much more dramatic- a significantly large fraction of the monomers belonging to 

# 

the tails are pulled into water along with the heads (see fig.7.18) to which they are attached! The 
loss of energy due to the increase in the area of contact between the hydrophobic tails and water 
is compensated by the gain of energy from the increase of contact between hydrophilic spacers 
and water as well as the gain of conformational entropy of the system arising from the larger 
amount of space available to those chains which remain at the interface. This interpretation is 
supported by our observation that this effect is more prominent at higher densities of surfactants. 
Some other manifestations of the entropic effect have been observed earlier [80-82]. 

The conclusions drawn from averaged concentration profiles are supported by the instan- 
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Figure 7.17: Concentration profiles for tail monomers for three different cases when the number 
of amphiphiles present in the systems are (i) 500 (solid line), (ii) 100 (broken line) and (iii) 10 
(dashed line). (a) hydrophobic spacer (b) hydrophilic spacer. 
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Figure 7.18: Concentration profiles for heads when the number of amphiphiles is 500. Solid line 
is for hydrophobic spacer and broken line for hydrophilic spacer. 
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taneous snapshots of the surfactants. Figure 7.19(a) 7.20(a) gives the snapshot of the gemini 
surfactants with hydrophobic spacer at the interface. For clarity we have shown only heads and 
spacers in fig. 7.19(a) and heads, tails and liaisons in fig. 7.20(a). This can be contrasted with 
the snapshots for hydrophilic spacer (fig 7.19(b) and 7.20(b)). In fig. 7.19(b) we have shown only 
heads and spacers and in fig. 7.20(b) we have shown heads, tails and liaisons. For hydrophilic 
spacer we find that the spacers and heads prefer to go inside water leaving space for tails at the 
interface. 


7.5 Summary and Conclusion 

In this chapter we have developed models of both ionic and non-ionic gemini surfactants with 
lydrophobic spacers as well as those with hydrophilic spacers. We have investigated the mor- 
jhologies of the micellar aggregates of these gemini surfactants and computed the corresponding 
IMCs by carrying out MC simulations. 

The main features of the aggregation of gemini surfactants with hydrophobic spacers can 
le summarized as follows; (i) the micelles are far from spherical- for short spacers these are long 
thread-like" and for long spacers these are "rod-like"; (ii) the CMC varies non-monotonically 
'ith increasing spacer length; (iii) the CMC increases with the increase of the bending stiffness 
F the tails and spacers. 

The main features of the aggregation of gemini surfactants with hydrophilic spacers can 
e summarized as follows: (i) the micelles are more or less spherical; (ii) the CMC decreases 
lonotonically with increasing spacer length; (iii) the CMC decreases with the increase of the 
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Figure 7.19: Snapshots of the gemini surfactants at the interface (a) with hydrophobic spacer 
and (b) with hydrophilic spacer. The black sphere represent monomers belonging to heads and 
grey sphere represent monomers belonging to spacers. 
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igure 7.20: Same as in figure 7.19 except that now the black sphere represents monomers 
elonging to heads, grey sphere represent monomers belonging to tails and white sphere represent 


lonomers belonging to liaisons. 
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bending stiffness of the tails and spacers. 

In contrast to the case of single chain surfactants the CMC increases with the hydrocarbon 
tail length for both the ionic and non-ionic gemini surfactants irrespective of whether the spacer is 
hydrophobic or hydrophilic. However like the case of single chain surfactants the morphologies of 
the ionic gemini surfactants are identical to that of the corresponding non-ionic gemini surfactants 
both for hydrophobic as well as hydrophilic spacer. 

Therefore, we conclude that (i) the shapes of aggregates are dominantly determined by the 
geometric shape and size of the molecules and whether the spacer is hydrophobic or hydrophilic, 
whereas (ii) the variation of CMC with spacer length is strongly influenced by the ionic charge 
and, again, whether the spacer is hydrophobic or hydrophilic. 

In the case gemini surfactants at the air-water interface for dilute regime, the cross-sectional 
area for single isolated gemini surfactant increases with the spacer length both for hydrophobic 
and hydrophilic spacer. However beyond a certain length of the spacer the cross-sectional area 
is larger for the hydrophobic spacer as compared to that for hydrophilic spacer. These trends 
are consistent with the variation of vertical extension < Z > with spacer length; a larger value 
of < Z > for hydrophilic spacer as compared to the hydrophobic spacer of identical length is 
observed in our simulations. 

For extremely ligh surface density of surfactants at the air-water interface we have demon- 
strated qualitatively the spatial organization of the gemini surfactants for both the case of hy- 
dropho’''c and hydrophilic spacers. 

In view of the above observations, it seems that the main effects of introducing the spacer 
is to impose an additional geometrical constraint on the packing of surfactant molecotes and. 
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therefore, to influence their aggregate shape and other properties. 

Molecular dynamics (MD) simulations of a similar molecular model of gemini surfactants has 
been carried out by Karaborni et al. [148]. In their model, particles of water interact mutually 
via a truncated Lennard-Jones (LJ) potential with sufficiently long cut-off to incorporate both 
the short-range repulsion and long-range attraction. The mutual interactions between the pairs 
of particles belonging to the tail were also similar. But, the cut-off range of the tail-water and 
head-head interactions were so short that no attraction was possible. However, the chains and 
spacers simulated by Karaborni et al. were much smaller than those investigated here. Besides, 
Karaborni et al. neither investigated the CMC and its variations with lengths of the tails and 
spacer nor considered any model of gemini surfactants with hydrophilic spacers. One should 
also try to develop more efficient MD algorithms to repeat our computations with more realistic 
interaction potentials on a continuum to check if any of the morphologies observed in this chapter 
have been influenced significantly by the discrete lattice. 

It would be interesting to investigate the effects of weakening of the screening (i.e., increasing 
the range) of the repulsive Coulomb interaction between the ionic heads on the results reported 
here; but, such a MC study will require much larger computational resources. 



Chapter 8 


Monte Carlo simulation of 
amphiphiles in complex fluids 
confined in porous medium 

8.1 Introduction 

Effects of confining surfaces on simple fluids have been investigated for several years, par- 
ticularly in the context of fluids in porous media [158]. However, the statistical mechanics of 
confined complex fluids has received very little attention so far. Investigations of these systems 
will not only contribute to the basic understanding of phase ordering in confined systems, but may 
also find application in oil extraction technology [160], because oil, brine (a solution of inorganic 
salts in water) and surfactants in earth’s crust form a complex fluid like microemulsions. 

Chowdhury et al.[159] have studied the structure of ternary microemulsions confined in model 
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pores through Monte Carlo (MC) simulation of Widom model. However, they have considered 
only single ’pores’ in the form of ’slits’. The real porous materials are known to consist of 
an interconnected network of pores of different shapes and sizes. Moreover the model they 
consider, the Widom model, has many drawbacks in describing the phase diagram of ternary 
mixture. Recently Liverpool and Bernardes [58] have studied the aggregation of amphiphiles in 
ternary mixtures confined by walls (which can be considered as single pore). They have used 
reformulated Larson model discussed in chapter 4. The main result of their study is that at low 
temperature almost all the amphiphiles form a spanning cluster. However at high temperature, 
the cluster breaks up into small clusters. In this chapter we propose a methodology for studying 
the equilibrium structure of ternary microemulsions confined in a porous medium (such as Vycor 
glass). We illustrate the use of our proposed method by applying it to Larson model of ternary 
microemulsions confined in a model porous medium with neutral confining walls. 

The chapter is organized as follows; in section 2 we describe our model.,and in section 3 we 
describe how a model vycor glass (which is a porous medium) is prepared in a computer exper- 
iment. In the subsequent sections we illustrate the methods we propose here for the computer 
simulation of microemulsions in porous media. 


3.2 The Model 

Our model is based on the Larson model of microemulsions [41, 42, 44] and its later 
leneralizations discussed in chapters 1 and 4 [46, 47], 

A detailed analysis of the distribution of the sizes of clusters of amphiphiles can be performed 
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using techniques developed for the analysis of cluster statistics in percolation theory [132], An 
amphiphile belongs to a cluster if any of its monomers is a nearest-neighbour of any of the 
monomers of amphiphiles already counted as part of the cluster. A Hoshen-Kopelmann cluster 
counting algorithm [72] calculates the cluster size distribution ti(s) of clusters with s monomers 
(units). 

8.3 Preparation of the Model porous medium 

The model porous medium used in our simulation is constructed by first allowing a binary 
mixture to phase segregate using the cell dynamics scheme (CDS) [162] and, then, etching 
out one of the two components. This method was suggested by Chakrabarti [161] for modelling 
Vycor-like porous media. The porosity of the medium, which is defined as the volume fraction 
occupied by the pores, was kept r)0(±l)%. In the following we outline the cell dynamics scheme 
[162] for generating the model porous media. The single cell dynamics is described by 

il}{t + l,n) = f{ip{t,n)) (8.1) 

where is the value of the order parameter in the cell n at time t. For the nonconserved 

order parametei^ the CDS model gives 

+ 1, n) = /('(/>(i, n)) -f D[« ij{t, n) » n)] (8-2) 

where D is a positive constant proportional to the phenomenological diffusion constant, << 
* >> — * is isotropized discrete Laplacian. We use the following definition of << * >> on the 
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three dimensional cubic lattice 


«ib{t,n)» = 


+ 

+ 


80 

80 

J_ 

80 


in the nearest - lu'ighlnnn- cells) 
in the next - lu^arest - ueighhour c;(dLs) 
^(■0 in the 3rd nearest - neighbour cells) 


(8.3) 


Since the choice of the injection function / does not affect the result^we choose the following 
form for the function / for our simulation 


f{x) = A tankx 


(8.4) 


where A\s3 parameter which is greater than 1. In our simulation we start from a x Ly x L, 
cubic lattice consisting of a 50% - 50% random mixture of two components. We use periodic 
boundary condition for the simulation. The values of the parameter A and J) were chosen to be 
1.3 and 0.5. The system was then allowed to evolve according to the cell dynamics equation (8.3) 
. A porous medium of desired pore thickness can be created by stopping the phase separation 
process after an appropriate number of iteration steps and etching out one of the two phases. 
A model porous media has been shown in fig 8.1(a), F igure 8.1(b) shows a cross section of 
this model porous medium. The porosity, which is defined as the fraction of the vacancies on the 
lattice, of the model porous medium is about 50%. 


8.4 Results and Discussion 

Once the model porous medium is created, oil, water and amphiphiles molecules were put 
inside the pore following Larson-like prescription, as explained before. If a site is occupied by the 
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Figure 8.1; (a) A model porous medium (b) A typical cross section. Dark grey represent rocks 


and light grey represents pores 



Figure 8.2; Cluster size distribution amphiphiles inside (a) single pore (c) vycor glass. Cluster 
size distribution for oil inside (b) single pore and (d) vycor glass, I lie solid litu,. is for temperature 
1 and dashed line is for temperature 4. 

rock It IS represented by S, = 0 (for non interacting rocks). If tire rocks rirc interacting then 
depeniiing on the nature of interaction we put S, =: I if the rocks .are "w.iler- Irjving" and S, = -I 
If It IS "oil-loeing". The amphiphiles were of the type first we simulated the system 

when oil-water-amphiphiies are confined by walls from the above and below (this corresponds to 
g e pore). We have studied the cluster size distribution of the amphiphiles for 2.5'‘ lattice where 
the total volume fraction occupied by the amphiphiles is lf,% (number of amphiphiles is 390). 
The ratio of oil and water was kept at 50% - 50%. At low temperatures (e.g.,7' == 1 .0) almost 
all the amphiphiles (331) form a huge cluster (figure 8.2(a)) [58]. At high enough temperatures 
(=.g..r = 4.0) the large cluster breaks into smaller clusters fig 8.2(a). When the simple confining 
replaced by the interconnected tortuous porous structure (vycor glass) the amphiphiles 
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Figiin' 8.3: (a) A cross section of bicontinuous microemulsion inside single pore, w denotes water 
molecule, o denotes oil molecule and H and T denote head and tail of the amphiphile molecule. 

cannot form a large cluster since the movements of the amphiphiles are now constrained by the 
presence of rocks (Fig. 8.2(c)). For low concentration (10%) of oil (water) droplets of oil in 
water (water in oil) are formed. Figures 8.2(b) and 8.2(d) give the cluster size distribution of oil 
in single pore and vycor glass respectively. 

When the concentration of oil and water are comparable we get bicontinuous structure for 
both single pore and vycor like porous meduim. Figures 8.3(a) and 8.3(b) show the bicontinuous 
microemulsions inside single pore and porous medium respectively. 

For low concentration (say 10%)) of either oil or water , micelles are formed at sufficiently 
low temperatures. Cross section of the miceller aggregates have been shown in figures 8.4(a) an 
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Figure 8.3: 


(b) A cross section of bicontinuous microemulsion inside porous medium. * denotes 


rocks. 
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Imriut 8.4: (a) Micellar aggregates (oil droplet) inside single pore. 


8.4(b). 


8.5 Summary and conclusion 

In this chapter we have proposed a method of computer simulation of the ternary microemul- 
sion confined in porous medium (such as vycor glass). The emphasis was on the methodology 
and not on the results. To illustrate the method we have calculated the cluster size distribution 
of amphiphiles and oil inside the porous medium at two different temperatures. We have also 
looked at the structure of the microemulsion for different concentration ratios of oil and water. 
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Figure 8.4: (b) Micellar aggregates (oil droplet) inside porous medium. 



Chapter 9 


Conclusion 


In this tlicsis wc have studied the bulk as well as interfacial behaviour of the amphiphilic 
systems from the microscopic point of view using statistical mechanical approaches, particularly, 
the techniques of MC simulation. We have considered two different approaches to study the am- 
phiphilic systems, namely, interface/membrane model and microscopic lattice model. However, so 
far as the membrane approach is concerned, we have studied only the out-of-plane thermal fluc- 
tuations of some model membrane. We have not studied any self-assemblies and bulk properties 
using membrane model approach. We have investigated the self-assemblies exclusively through 
MC simulation of microscopic lattice models. 

In chapter 2 of the thesis we have studied the nature of the amphiphilic monolayer membrane 
at the oil-water interface in the Widom model which is the simplest microscopic lattice model of 
ternary microemulsions. The Widom model is a classical Ising spin model with nearest-neighbour 
as well as farther-neighbour interactions on an appropriate lattice. The nearest-neighbour inter 
actions gives rise to interfacial tension and the farther-neighbour interactions contribute to the 



Conclusion 

bending stiffness of the amphiphilic monolayer; these farther- neighbour interactions give rise to 
a richer bulk phase diagram of the Widom model, as compared to that of the nearest-neighbour 
Ising model. We have found that these farther-neighbour interactions do not affect the interface 
roughening exponent i.e. the membrane is rough but not crumpled. Nevertheless, the oil-water 
interface in the Widom model is wider than that in the ising model with only nearest-neighbour 
interactions. 

The absence of the crumpling transition in the Widom model is attrit)uted to the fact 
that in the Widom model we can not switch off the interfacial tension and study the effect of 
bending stiffness only. On the other hand liquid-like amphiphlic membranes whos(' out-of-plane 
thermaf fluctuations are governed by the bending stiffness alone are known to be crumpled. To 
make a clear distinction between the interfacial tension and bending stiffnes.s and to get a possible 
crumpling transition we have proposed a discrete model of an interface (botfi unrestricted discrete 
SOS (UDSOS) and restricted discrete SOS (RDSOS)) with a specific form of the l)onding energy; 
this bending energy is associated with bending of the interface at right angles. These have been 
discussed in chapter 3. When the out-of-plane thermal fluctuations of the interface in the UDSOS 
model are controlled only by this bending stiffness it is crumpled but the nature of the crumpling 
is very different from that observed in the continuum models. We have shown that although 
the conformation of the interface in this case satisfies one particular definition of a crumpled 
manifold, it does not satisfy two other widely used definitions of crumpling although all these 
three definitions are known to be mutually consistent when applied to the continuum models. 
We have studied some dynamical properties of a UDSOS model in d =. 1 in the special case of 
vanishing interfacial tension but non-zero bending stiffness by using a Master equation approach 
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and the width of the interface U’ is found to vary with time as v^. By identifying an appropriate 
crossover length, we liave also explained why the roughening exponent a of the UDSOS mode! in 
the presence of non-zero interfacial tension is independent of whether or not it has non-vanishing 
bending stiffness. 

Introduction of a new RDSOS model is motivated through the examination of four different 
ensembles. By carrying out MC simulation of this RDSOS model, we compute an exponent 3 
which characterizes the time-dependence of its width during the very early stages oi \ts dynamkzl 
evolution from a flat initial configuration. 

So far the bending rigidity has been introduced as a parameter in all the models of mem- 
branes. In real amphiphilic systems this arises naturally from the microscopic interaction among 
the constituent molecules and from their conformational degrees of freedom. So it is natural to 
consider a more realistic model where the conformational degrees of freedom of the amphiphilic 
molecules are taken into account. In chapter 4 we have studied the roughening properties of the 
amphiphilic membrane at the oil-water interface in the Larson model which explicitly incorporates 
the conformational degrees of freedom of the amphiphilic molecules. We have shown that the 
width of the interface varies as v/L;,:, and hence the oil-water interface in the Larson model is 
rough but not crumpled. Moreover, the wider interface in the presence of amphiphiles in this 
model is a consequence of the lowering of the oil-water interfacial tension by the surfactants, 
we have argued that this reduction of surface tension in this model is of entropic origin. After 
studying the nature of amphiphilic monolayer at the oil-water interface we shift our emphasis 
towards air-water interface and, then, study some of the bulk properties the amphiphilic systems. 

In chapter 5 we have investigated the out-of-plane phase segregation of the amphiphiles in 
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a binary mixture, of two species which are chemically identical but has different lengths, at the 
air-water interface. We have demonstrated an entropy-driven [)has(' segregatioti in a direction 
perpendicular to the air-water interface when the initial total surface density of the amphiphiles 
is sufficiently high. The steric repulsion between the tails of the longer amphiphiles is argued 
to be the driving force behind this phase segregation phenomena. In other words, this phase 
segregation is driven by the gain of conformational entropy, contrary to tlie phenomenon of 
energy-driven phase segregations, e.g., in a binary alloy. We have <ilso studied the nature of the 
in-plane clustering of the monomers belonging to the long and the short amphiphiles. 

So far we have studied only the interfacial properties of the systems cunt .lining only single 
chain amphiphiles although in chapter 5 we have studied sonic of the bulk firoperties. In the 
subsequent chapters we have focussed attention on both the interfacial as well as bulk properties, 
in chapter 6 we have looked at the morphologies of the aggregates formed by single chain and 
double chain amphiphiles. Our computer simulations have shown the formation of stable spherical 
micelles for single chain amphiphiles both for single species as well as for binary mixture of single 
chain amphiphiles. On the other hand double chain amphiphiles have been found to aggregate 
forming vesicles. 

Gemini surfactants which are complex double chain amphiphiles are thought to be the next 
generation surfactants for their unusual characteristics, such as very low critical micellar concen- 
tration (CMC) and a high efficiency in reducing the oil-water interfacial tensions in comparison 
with the conventional single chain surfactants. We have proposed a microscopic lattice model 
of these gemini surfactants in chapter 7. By extensive MC simulation we have studied the the 
^SS'’sg3t>on behaviour of these gemini surfactants in water and the spatial organization of these 
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gemini surfactants in a system where water is separated from the air above it by a sharp well 
defined interface. Wo have also calculated the critical micellar concentration (CMC) by varying 
different parameters. The results of our numerical investigations are summarized as follows: (i) 
For hydrophobic spacer these model geminis form long, thread-like and entangled micellar aggre- 
gates even at low concentrations, where single-chain model surfactants form spherical micelles, 
provided the spacer is not too long. On the other hand, rod-like micellar aggregates form when 
the hydrophobic: spacer is sufficiently long. When the spacer is hydrophilic we get more or less 
spherical micelles. There is no ciualitative difference between the morphologies of the aggregates 
formed by ionic gemini surfactants and those formed by the corresponding nonionic gemini sur- 
factants, irrespective of whether the spacer is hydrophobic or hydrophilic. 

(ii) For ionic gemini surfactants with hydrophobic spzcer the CMC varies non-monotonically with 
the increase in spacer lengtfi, whereas it decreases monotonically with spacer length for non-ionic 
gemini surfactants. On the other hand, when the spacer is hydrophilic the CMC decreases mono- 
tonically with spacer length both for ionic as well as non-ionic gemini. (iii) The CMC is found 
to increase with the increase of the bending stiffness of the tails and spacers if the spacer is hy- 
drophobic. in contrast, the CMC dcrnascs with the increase of the bending stiffness of the tails 
and spacers when the spacer is hydrophilic, (iv) In contrast to the single chain surfactants the 
CMC of gemini surfactants increases with the tail length (both for ionic and non-ionic) irrespective 
of whether the spacer is hydrophobic or hydrophilic, (v) In case of the gemini surfactants at the 
air-water interface for extremely dilute regime the cross-sectional area of the isolated individual 
gemini surfactants increases with the spacer length for both hydrophobic as well as hydrophilic 
spacer. We also find the similar trend of variation for the vertical extension. However when the 
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spacer is long, for a given value of spacer length, the cross-sectional area is smaller if the spacer 
is hydrophilic, whereas the vertical extension is larger for the hydropliilic spacer as compared to 
the case when the spacer is hydrophobic. For the case of high surface density regime we have 
shown the spatial organization of the gemini surfactants qualitatively for both the case when the 
spacer is hydrophobic as well as for the case when the spacer is hydrophilic. Some of the trends 
of variation are counter-intuitive but are in excellent agreement with the available experimental 
results. 

From the observations in chapter 7 we infer that (i) the shapes of aggregates are dominantly 
determined by the geometric shape and size of the molecules and whether the spacer is hydropho- 
bic or hydrophilic, whereas (ii) the variation of CMC with spacer lenglli is strong, ly influenced by 
the ionic charge and, again, whether the spacer is hydrophobic or hyclrophilu. 

In chapter 8 we have presented a numerical scheme for computer simulation of ternary 
microemulsion inside porous medium, e.g., vycor glass. 


9.1 Future direction of research 

It is important to calculate the gain in conformational entropy over the energy loss quantita- 
tively and to show that it is the driving force behind the phase segregation phenomena described 
described in chapter 5. Future research also should try to explain quantitatively (i) the non- 
monotonic variation of CMC with spacer length and (ii) increase of CMC with the increase of 
hydrophobic tail length for gemini surfactants. To explain the unusual morphologies of the ag- 
gregates formed by the gemini surfactants, both in water and at air-water interface, is another 
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challenging problem. 

The - B block copolymer serves as a surfactant in a X/B homopolymer blend, A binary 
polymer blend containing an amphiphilic block copolymer shows similar phase separated struc- 
tures like microemulsions [163], It would be interesting to apply our methodology of modelling to 
such ternary polymer mixtures and study the phase behaviour and to see whether any interesting 
new phenomena can arise. Recently there has been lot of interest on 'living polymers’. There is 
substantial experimental evidence that for low surfactant concentration certain surfactant solu- 
tions contain large one dimensional aggregates which can be considered as ‘living polymer’ [164]. 
So it would be interesting to study the structure and dynamics of the "living polymers” using the 
lattice models of the same type as studied extensively in this thesis. 
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SPIN GLASS THEORY 
IN NEURAL 
NETWORKS 


8.1 Introduction 

8.1.1 Generalities 

Bioiogicai computers - the brain aind the nervous system have existed for 
millions of years and they are marvellously effective in processing sensory 
information and controlling the interactions of the living systems with their 
environment. Seeing something, recognizing a face or recalling from mem- 
ories takes an enormous amoimt of computation. For example, to “see" 
something, first the images received in our eyes are coded by our nervous 
system in fantastic richness of details, then processed and restored in our 
condousness, corrected for changes of illmnination, mood and any other 
circumstances. The fact that bioiogicai computation is so effective suggests 
that it may be possible to attain similar capabilities in artifidal devices 
or at least it is possible to model the brain where the peculiarities of its 
performances emerge naturally. 

Although our brain functions as an information processing computo:, 
the basic prindple on which its functioning is based is quite differenc from 
that of a conventional computer. A conventional computer stores informa- 
tions by assigning addresses, which identify the physical locations where the 
data will be stored in the hardware. When the central processor requires 
a piece of data, it issues an instruction to read. the data at a particular 
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A spin-rocationaily-symmetric patk-integrai formulation, with a vector order parameter is pre- 
sented for a generalized A' -orbital Hubbard model which jdeids the correct Hartrcc Fock result in 

the saddle-point approximation, and also the spin-wave mode representing transverse Gaussian fluc- 
tuations. A systematic expansion of the action in powers of inverse .V, wherein spin-rotational 
svxnmetrv and hence the Goldstone mode are preserved order by order, is developed, and the first 
order 0(l/-\f) quantum corrections to the ground-state energy and subiattice magnetization ob- 
tained. Setting = 1 yields a quantitatively correct description of the antiferromagnetic state of 
the one-orbital Hubbard model at half- fill ing. 


L INTRODUCTION 

Recently there has been interest in developing a spin- 
rotationaily-symmecric path-integral formulation of the 
Hubbard model. Up until recently the situation has 
been that the process of taking the saddle-point ap- 
proximation has led to inconsistencies in the resulting 
theor}'.^ Thus if the Hubbard interaction term is decom- 
posed in terms of charge-density and spin-density terms. 
TtitTiii = ~ {riif — resulting the- 

ory' yields correct HF energies, however, the presence of 
a scalar magnetic order parameter (n,-t — nn) in the the- 
ory leads to the absence of Goldstone modes. On the 
other hand, the — |St -S,- representation introduces a vec- 
tor order parameter (S,*), and hence Goldstone modes in 
the theory, however, because of the | factor, the single- 
panicle spectrum obtained at the Hanree-Fock level is 
incorrect. 

In the present paper we show that, using a generalized 
^V-orbital Hubbard modei.'^*^ it is possible to develop a 
path-integral formulation wherein the effective action is 
both spin-rotationaily invariant and exhibits a correct 
Hartree-Fock result at the saddle point. For the Hub- 
bard model there is no formal expansion parameter in the 
problem to control the perturbaci\''e expansion, whereas 
in a perturbative treatment it is imponant to ensure that 
adl conser\'ation laws and symmetries are obeyed system- 
atically. La particular the Goldstone mode which fol- 
lows from spin-rotational symmetry must be preserved 
in perturbation theory. In the path-integral formulation 
of the generalized -V-orbitai Hubbard model acts as 
an expansion parameter and spin-rotational symmetry 
is guaranteed order by order in the expansion becauuse 
the perturbative piece of the Hamiltonian is manifestly 
rotationaily symmetric. Ln the saddle-point approxima- 
tion, which rieids the classical action of 0(1), we re- 
cover the correct HF result, and this is also exact in the 
limit -T oc. The quantum fluctuation effects appear 
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at higher orders and we consider Gaussian flue 
around the saddle point yielding an 0(1/ j\f) com 
the action, from which we discuss evaluation of t 
turn corrections to the ground-state energy and st 
magnetization. . 

We consider the following Hamiltonian for a 
ized j\f orbital Hubbard model with j\f orbitals 
site:*^’^ 

-t 

(ij)cra 



Here a.p are orbital indices. The factor ^ i 
in the interaction term to render energy density 
the j\f CO limit. The two interaction terms are 
density and exchange type wiihxespect to orbital 
and together can be expressed as^ 

^ V ( ^ Ita Tot ^ \xB i/3 "h 

ai5 


where “La). S,- = is t 

spin-density operator, and n,* = Yla ^ 

averaged charge-density operator at site L Tlie 
tion term, expressed as above with the expansiot 
eter -y in front, is manifestly spin-rotationally sy 
and hence the Goldstone mode is preserved ord< 
der in a systematic ^ expansion. 


IL PATH-INTEGRAL FORMULATIO 

To set up the path-integral formulation of the 1 
model we go to the Grassmann-'variable represeni 


0163-1829/94/49(9)/6078(5)/S06.00 


6078 


© 1994 The American Physic 



;r aniisymmerry of the many-body fermioii 



Now to get the HF result we obtain the action in 
the saddle-point approximation, wherein the action is 
minimized with respect to variations in the field, i.e.. 
= 0, where ^ = {x,y,z). At the saddle point we 

obtain 


( 3 ) 

.d ^ia column vectors containing Grass- 
)les for the two spins, corresponding to the 
rators and Here we have dropped 

which simply shifts the position of the Fermi 
>w the partition function 2 = Tre*^^, in 
id at zero temperature, can be expressed in 
egral formalism as a functional integral over 
variables: 


where we have defined the single-particle Green's func- 
tion in the two-spin basis as 


Gocih ri ^ ) 



( 11 ) 


We define the local magnetic moment due to an electron 
in orbital a by 




.agrangian density is given by 

{ij)a 



(Siot(f)) — iTror ^Got(x, t; I, t) (12) 

and hence at the saddle point the vector field Is noth- 
ing but the local magnetization 

4>i^{t)^±{Sio.{t)). (13) 


( 3 ) 

III. HARTREE-FOCK RESULT 


dth the quartic interaction term, we make use 
Dard-Stratanovich transformation, and intro- 
vector field Oi{t) = Yla each point 

d time. 


= const J 

(6) 

i>{{t) couples with the total spin at site i 
ver all orbitals). The interaction term is now 
ladratic and it is now possible to integrate out 
dc field to get an effective action in terms of 
(t). Integrating over the Grassmamn variables 
for the partition function,'^ 

(7) 


effective action is given by 
U 


-I 


(it 


i 

^ TrlnuSt - Sacicj>)] 


( 8 ) 


operator Hai<p) is the effective one-particle 
tn for orbital a and in the two-spin basis is 


Now we show that within the saddle-point approxima- 
tion we get the correct HF result for the half-filled system. 
For half-filling we expect an antiferromagnetic ground 
state, which is characterized by an ordering wherein the 
local magnetization has equal magnitude but opposite di- 
rection on sites of the two sublattices. We choose the c 
direction as the direction of ordering without any loss of 
generality. Thus, we have = 0, and 

<(>L = ( 14 ) 

where m = |nT^ — the staggered or subiattice mag- 
netization is the order parameter in the AF state. The 
sum over orbitals in the total field <pf in Eq. (9) yields a 
factor AG and therefore, at the saddle point, the effective 
one-particie Hamiltonian may be written, again in the 
two-spin basis, as 

23(|x)(ii +H.c.)l + 23(-1 )'Ao-", (15) 

m i 


where A = \[J\fm)U/J^ = XjTnjl is the gap parameter. 
The corresponding Hamiltonian in momentum space, and 
now expressed in the two-sublattice basis, can be written, 
for spin-ij electrons, as^ 



€k 

crA 


(16) 




(ij) 



where e\c = — 2t cos(A:^a) is the free-particie energy' 

on the lattice. T he eigenva lues yield the quasiparticie en- 
ergies, Ek = ±-y/A2 -h €|, and electronic amplitudes on 
sites of the two sublattices are obtained from the eigen- 
(9) vectors (ufecr hka)- The sublattice magnetization, which 
is the absolute difference in electronic densities on sites 



of the nvo sublactices. is given by 


=-?r73^. (17) 

u ;c V - er 


Here rhe sum is over all occupied states lying in the lower 
Hubbard band, and A is the number of lattice sites. With 
2A = mC. this leads to the seif-consistency condition 


_L-.iL ^ 

U N 4- 


(18) 


from which the (HF) sublattice magnetization may be 
determined as a function of the interaction strength UJ 


IV. SPIN WAVES 

We now study the spin-wave modes, which are low- 
lying, collective excitations in the antiferromagnet, rep- 
resenting transverse fluctuations around the broken- 
symmetry direction. We do this by allowing for small 
transverse fluctuations in the vector field around 

the saddle-point value. 

4>ia = (19) 

Since 4>i = (pict, summed over aU orbitals, we have, 
for substitution in Eq. (8), 

0 ? = 4 - 20 -^^ • ^ ocbia +• ^ S(f>ia • S<i>{0, 

ai3 

( 20 ) 

Correspondingly Hcxi<P) being linear in 0 can be written 
as 


= +SH^, (21) 

where represents the effective one-particle Hamilto- 
nian [Eq. (15)] at the saddle point, and in the two-spin 
basis. SHct is given by 

I 


U 

d'Ha = ^ S(pi0 - cr|f)(zi. 




ia 


By making the above substitution in Eq. (8j 
panding in powers of we obtain an expansio 
effective action S^q{o] around the saddle point 
ers of 6(j>i{t). Within the Gaussian approximati 
keeping terms only up to second order, we haxe 

where is the (classical) action evaluated at the 
point: 


/*[-f i;.v7»r)= 




which is j\f times the classical action for the one-t 
model, and in terms of the Green's fnnctioi 
(11)] at the saddle point, we have 



dt 


j\f ^ 

iaS 


. u 
+1 ^ 


T=i 


Since we are interested in spin waves, we cor 
only transverse fluctuations, i.e.. d'of = 0. No 
is convenient to work with bdr- which are denm 
(Sd)^ ±iS6^)/V2 and z: rT^)/2, instead c 

X and y components. The term in SH is ther 

expressed as y/2{Sdf^cr^ +S0'^cr'^). Substituting for 
from Eq. (22) and performing the sum over the ini 
space and time coordinates coming from the trace i 
second term, we can write this as 


ij aJ • 


where the superscripts s and s refer to the two signs, s = — s — h, and the matrix elements oi D are gi’ven be 
D being given analogously, 

V 

7=1 

= SijS{t - t') - iUGi°\it,jt')Gf\jt'.it). 

Fourier transformation to frequencv space and momentum space (the latter within a two-subiattice basis as tr 
lationai symmetry is present only within this basis) yields the foUowing expression for the quantum correction tp 

action: 




;race is over the twosubiattice eiements. and 
= 1 - Ux^iQ.^) = -0). (29) 

Q.f>) is the zeroth-order. antiparailel-spin 

2) = iY^ f ^G'°\k.u^)Gl°Hk-0.^-.n}. 

k 

(30) 

D~[Q, n) is now rendered in a form familiar 
mdom-phase approximation (RPA) expression 
isverse spin snsceptibilit}^ Poles of £?“(<?, f2) 
pin-wave energies at the RPA level.' 
discuss the order of relative to that of 
al action Assuming no correlation be- 

nations in different orbitals, i.e., ) = 

)”), we see that the sums over a./3 contribute 
:tor of j\f. which cancels the j\f in the denom- 
Iding a contribution which is 0[l/\') relative 


to the classical action. In fact, generalizing to the nth- 
order term in the expansion for the action, with n being 
even for the averaging over fluctuations to suni^^e, we 
have n factors oiM in the denominator, coming from the 
explicit factor of l/M in 5H. and n sums over orbitals. 
Again, the absence of correlation between fluctuations 
in different orbitals requires the n = 2m orbitals to be 
paired off into m pairs, leaving jn independent sums over 
orbitals, yielding a factor of The net order of the 

nth-order term is therefore {1/jV)'^^^, so that the expan- 
sion of the action, into quadratic, quartic. etc., terms, is 
formally seen as a perturbative expansion in powers of 
l/AT: 


5 = 


Al 


Af^ 




(31) 


Therefore, up to the (Gaussian) level of small trans- 
verse fluctuations aroimd the saddle point, we have for 
the partition function 




V(56^)V{50'^) exp 


Q 5==‘ 


dn 

2tt 




I over the bosonic fluctuations, we obtain, 


xp 


u\r + ^E / ^Tr;lnZ?^-(Q,n) 


(321 


(33) 


trace is over the two-sublattice eiements within 


GROUND-STATE ENERGY AND 
UBLATTICE MAGNETIZATION 


)ve correction to the action leads to following 
to the -ground-state energ}': 

E / (34) 


Af) (quantum) correction to ground-state en- 
0 transverse Gaussian fluctuations around the 
at is. in a diagrammatic language, nothing but 
cion due to transverse spin fluctuations at the 
. This correction has been evaluated in the 
y tractable strong-coupling limit. ^ Up to this 
fluctuations lower the ground-state energy rel- 
ae HP \aiue due to mixing between the HF 
the spin-wave states. Up to the Oil/Af) level 
i-state energy per site is^ 


-J 


^ ■ ,\rE 


(35) 


1 

where iQ ^ ^ The first-order quantxun cor- 

rection to. the ground-state energy is about 16% for the 
square-lattice case. 

We now consider the quantum correction to the sub- 
lattice magnetization due to the Gaussian fluctuations 
around the saddle point. Now, the ground-state en- 
ergy correction discussed above corresponds, in a dia- 
grammatic language, to the ladder diagram consisting of 
ladders of interaction lines between fermion loops of op- 
posite spin. Cutting fermion lines of both spins at the 
same ladder step yields the ladder-sum (RPA) contri- 
bution to the transverse spin susceptibiity, and cutting 
just the i-spin (say) line yields the j^A spin-fluctuation 
correction to the t“Spin fermion propagator. The cor- 
rection to subiattice magnetization due to this (RPA) 
spin-fluctuation correction has been evaluated again in 
the strong-coupling limit.® At the 0(l/j\f) level there is 
a 40% reduction for the square-lattice case, yielding a 
sublattice magnetization of 0,6. Due to an exact cancel- 
lation of the 0(1/A/)^ contribution, this result actually 
holds to the two-loop level.® 

These 0[1/Af) corrections yield results which are iden- 
tical to the 0(1/25) corrections obtained within the spin- 
wave theory for the spin- 1/2 quantum Heisenberg AF-^ 
This is not surprising as in the strong-coupiing limit 
the A/"-orbital Hubbard model is equivalent to the spin 
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S = A'/2 Heisenberg AF. The systematic expansion in 
powers of 1/A' is perfectly equi\-aient. order by order, to 
the 1/25 expansion. 

In conclusion, we have presented a spin-rotationally- 
symmetric path-integral formulation with a vector order 
parameter for a generalized AAorbital Hubbard model. 
This procedure yields the correct Hartree-Fock result 
in the saddle-point approximation, and also the correct 
spectrum of spin-wave modes. A systematic expansion 
in powers of inverse j\f is developed, and the first-order 
quantum corrections to ground-state energy and sublat- 
tice magnetization are obtained. As mentioned in the 
Introduction, the 1/A/* expansion serves to systematize 
the fluctuation corrections about the HF state, and set- 
ting A* = 1 yields a quantitatively correct description of 
the AF state of the one-orbital Hubbaird model at half- 


filling. 

In view of the above obser%-ation it m uartint 
quire whether the ground state of the doped sp 
the large-A' limit might also <iescril)e the A" - I 
whether there are featuress ms.sornjued with thr one 
model which are missing from the largr-A'' modi 
analysis of the doped system in the Iarge«A'' limit 
presented elsewhere. We simply note here that 
havior of holes in the ground state of the large-A 
is quite different from that in the A/* = 1 case, 
ticular. the sequence of hopping events initiated 
added holes, leaving strings of upturned spins ab 
AF arrangement, is absent in the huge- A’’ modi 
inhomogeneous magnetic: grounci states like the 
mensurate domain- wail states are therefore re 
ground states of the doped larg^e'A^* mociel, 
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IMPROVED PERFORMANCE OF THE HOPFIELD AND 
LITTLE NEURAL NETWORK MODELS WITH I'lME 
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We report the results of simulation of neural lu'twork iiuxIpIs with tin* fiyuaplu nuuiet 
tions constructed using the IIc4;)h’s rule anti the <lynaini<.*i det<'ruun<‘<l by ih*' ujtftnal 
field, which has a weighted contril)ution from the time ilrbiyod ragmvhi \VV i 
both the asynchronous (or Glaubtu’; nopli<4<l) and syni hronoufi (bitth*| dyiunuM i hu 
numerical results and the finite siz<‘ variation stmly (for si 7 .«‘r) S withm ihr 
< N < 4000) support the pr<‘vious imlh ation |S<*n and < !haktab,it<i, /Vit/.s A 

327 (1992)] of improved performaiu<‘ in the re( all and overlap priipei tter^ m (he ihnmo 
dynamic limit. It is identified that the time delayetl term m the dynamn ri alhavr* ib** 
network to come out of the spurious valh^y.s in tin* 'b‘nergv landra ape'* plelmed withmu 
the delay term; Hopfield mo<l<‘l). In an a(>proximate analytn al ,»itudy ol ran h modefij 
in the extreme dilution limit, tin* roh* of tin* time dehtyed term (o .suppierei the ('apm 
glass-like) noise is also indirat(Hl. 

1. Introduction 

The properties of associative iiHuaory models like the Hopfield ami the Little models 
are now well studied. ’ In .such models, (*a<di two .state (McCullouf.h and Ihti:.) 

neuron is represented by the Ising spin variabU's a, ■; 11; i 1, A'. The 

symmetric synaptic interactioms betwwrn the mniron.s aiv cou.structed following the 
Hebb’s rule of learning, which say.s that for p nund)er of patterns the syiinpHr 
strength Jij for the pair of neuron.s i and j i.s 

EC <;■, (II 

/i I 

where {Cfi}, i = 1, 2, ... , V r<!l>r<^slmt.s the ;/th p.attem to be learned by the 
network. Each {Cf } can takts valmis .tl. N is the total number of nemom., eaeh 
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connected to all other neurons, and p = aN denotes the total number of patterns 
to be learned. Any arbitrary pattern {Si} representing the state of the neurons at 
time t then evolves following the dynamics 

Si(f + 1) = sgn {hi (t)) , (2) 

where hi is the internal field on the neuron i due to all other neurons and is given 
by 

hi{t) = ^ ^ Jij Sj{t) . (3) 

j 

The “associative learning” of (p = aN) number of random patterns (represented by 
{,^'‘} n = l,2, ... ,p) corresponds to the creation of that many attractors or fixed 
points (may also be limit cycles in the Little model) of the dynamics. In the case 
of sequential updating in (2), or the zero temperature Glauber dynamics as in the 
case of Hopfield model, one can define an energy-like function 

N 

n{t) = -'£jijSiit)Sj{t). (4) 

i>j 

This energy function can be shown to decrease monotonically with time t, guaran- 
teeing a fixed point of the dynamics such that after a finite number of iterations 
(t*), the network dynamics stabilizes for all i and the fixed point equation is given 
by 

Si{e + i):=Si{e). (5) 

For parallel updating in (2), as in the case of the Little model, we can define a 
Lyapunov function^ 

N 

Cit) = -Y,JijSi{t)Sjit-l), (6) 

i>j 

which can be shown to decrease monotonically with time t, guaranteeing a limit 
cycle or fixed point of the dynamics such that after a finite number of iterations 
(t*) the network settles to the state 


1 ) = Si{t^ ^ 1 ) 


( 7 ) 


for all L 

Because of the possibility of defining the energy function for the Hopfield model 
with sequential or Glauber dynamics, the static properties of the model are well 
known and well understood from statistical physics.^’^ In particular, the memory 
loading capacity (with significant overlap of the fixed point states with 
the learned patterns Cii starting from any distorted pattern with initial overlap 
< 1 within the domain of attraction) shows a discontinuous transition at Oc — 
0.138 (m^ > 0.97 for a <ac and = 0 for o: > ac for any fi). The dynamics shows 
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an unusual critical slowing down near In particular, the average relaxation 
or recall time r (average of t* over various patterns, for a fixed amount of initial 
distortion) shows an unusual form of growth near a = ac, although there is no real 
(static)phase transition (the recall or overlap states are metastable for a > 0.05, 
where the minimum energy are the spin glass or spurious and the metastable overlap 
states disappear after a 0.14).^’^ Although no analytic theory could be developed 
(for finite a), the static and dynamic properties of the fixed points (recall states) of 
the Little model are almost identical to that of the Hopfield model.^’^ 

In the recent literature, there has been a number of indications^’® that dynami- 
cally defined networks , with asymmetric synaptic connections Jij may have better 
recall performance because of the suppression of spin glass-like states. The exactly 
solvable model of Derrida et al^ indeed uses the same Hebb’s rule for connection 
strength Jij but with extreme dilution (inducing asymmetry). This gives better 
recall properties. There were also indications^’® that addition of some local mem- 
ory of each neuron, in the sense that the internal field hi{t) in (3) is determined 
by the cumulative effect of all the previous states of the neighboring aj{t - t?), 
?7 = 1, 2, . . . , t (as indeed is biologically plausible; Peretto^), gives better recall 
dynamics in some analog network models. In all these cases, no effective energy 
function can be defined (because of the asymmetry of Jij or of the local memory 
effects in hi) and the use of statistical physics of spin glass-like systems is not possi- 
ble. All the networks are thus defined dynamically. In a recent simulation study®’® 
it has been shown that a dynamics with a single time step delay with some tunable 
weight A . 


Si{t + 1) = sgn 




( 8 ) 


j 


instead of (2), improves the performance of the Hopfield-like model (with sequential 
updating) enormously. Here the same synaptic connections Jij in Eq. (1), obtained 
by using Hebb’s rule, are employed. In particular, for a network with N = 500, the 
simulations indicated® that the discontinuous transition for the overlap function at 
ac 0.14) disappears and becomes continuous. Accepting recall with final overlap 
TTif > 0.9, the effective threshold loading capacity Oc increases from 0.14 for A = 0 
to about 0.25 for A 1. It was also found that average recall time or relaxation 
time^ T (defined before) becomes minimum around A ~ 1 at any particular loading 
capacity a. 

The above interesting and somewhat surprising results for this single step time 
delayed dynamics®’® are from simulations of a rather small size (N = 500) of the 
network. Here we present the results of a detailed numerical study of this dy- 
namics, both with asynchronous (sequential as well as random; Hopfield-like) and 
synchronous (Little-like) updating (Sec. 2). We check our numerical results for a 
bigger network (for 250 ^ N ^ 4000) and study if there is any significant finite 
size effect. We find that the fixed point memory capacity improves considerably 
with A, both with Hopfield and Little dynamics (asynchronous and synchronous, 
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respectively )and that the performance saturates after A = 1. With 90% final over- 
lap (rrif > 0.90) with the learned patterns, the loading capacity increases to almost 
ac 0.25 for A > 1. The relaxation or average recall time r, for any fixed loading 
capacity a (< 0.25), is also minimum at A = 1. We could not detect any significant 
finite size effect (for 250 < AT < 4000) and we believe that this result is true in the 
thermodynamic limit. We also observe that for A = 1, both the Little dynamics 
and Hopfield dynamics give identical results (all the limit cycles of the Little model 
disappear to give fixed points). The performance of the network for memory storage 
in limit cycles for negative A has also been studied for Hopfield and Little dynamics. 
Here, the negative A term provides a damping^® and limit cycles become frequent 
but the overlap properties deteriorate significantly. We have also studied the per- 
formance of such dynamics for randomly diluted networks. In Sec. 3 we give some 
approximate analytical formulation. In particular, in the extreme dilute limit,® the 
better performance of the network is indicated by analyzing the evolution equation 
of the overlap function for the overlap states, by using local field analysis technique. 
In Sec. 4, we discuss the significance of these results. 


2. Numerical Simulation and Results 

We consider a network of N neurons (spins). The synaptic interactions between 
neurons i and j are found using the Hebb’s rule (1) for a set of p-random Monte 
Carlo generated patterns {C^}, p = 1, 2, . . . , p; z = 1, 2 . . . iV. After this ‘learning’ 
stage, each pattern is corrupted randomly by a fixed amount (typically 10% of the 
neuron states are changed randomly; with initial overlap = 0.90 for any pattern 
p) and the dynamics (updating process) following (8) starts with these corrupted 
patterns. Both sequential (Hopfield) and parallel (Little) updating are employed. 
Since our updating process requires two initial patterns (at {t - 1) and {t — 2)), 
we start with two randomly distorted patterns with the same value of distortion 
(for each p). We study the overlap of the network state with the pattern p during 
the updating process and the final overlap of the fixed points (checked for 
two successive time steps) or limit cycle patterns are noted. The averages over all 
p patterns (and over 2 to 5 sets of random number seed for the entire run) are 
taken {rrif = average of over all p). The number of iterations (t*) required to 
reach the fixed points (or limit cycles) are also noted and their average (over the 
patterns) give the average recall or relaxation time (r). Our typical results for the 
Hopfield model (dynamics) are for N = 1000 and those for the Little model are for 
N = 500. We consider mostly positive A values. Negative A values do not improve 
the performance and have been studied just for some typical investigations. 

In Fig. 1, we give the variation of the final average overlap nzf (of any pattern) 
for the fixed point (checked over two successive time steps) with the redefined A 
[Ar = (1 - A)/(l = A)] at fixed loading capacity a = p/N 0.25 for sequential 
as well as random (Glauber) updating. It may be noted that for the Hopfield 
dynamics (A = 0) the final overlap rrif < 0.30 for a = 0.25. Strictly speaking, this 
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Fig. 1. The variation of the average final overlap mf (for the fixed point obtained) with Ar(= 
(1 — A)/(l + A), positive A) at a fixed loading capacity a = 0.25, with sequential dynamics for 
networks with N = 250, 500 and 1000. The inset shows the variation of average recall or relaxation 
time r against Ar for diiferent N. 


value should be zero in the thermodynamic limit. The results shown in Fig. 1 are 
for N = 250, 500 and 1000, and the nonvanishing value of rrif for A = 0 (Ar = 1) at 
a = 0.25 is precisely due to this finite size effect as can be easily checked (decreases 
considerably as N becomes large). It is seen that the final overlap mf gradually 
improves with increasing A and it saturates beyond A = 1 (Ar = 0) to mf 0.90 
(for a = 0.25). The inset figure shows that the average recall time (or relaxation 
time) r also decreases as A increases from A = 0 and finally attains an optinal 
value (dependent on the network size N and loading capacity a) at A = 1 (Ar = 0), 
beyond which it increases again. The above results for r are for sequential updating. 
For random updating, the qualitative behavior is again the same (although the 
magnitude of r is somewhat larger). 

In Fig. 2, we again plot the final average overlap mf (for the fixed points) with 
Ar for parallel (Little-like) updating. The results for the sequential (Hopfield-like) 
updating are also given here for comparison. The results for the fixed points are 
qualitatively and quantitatively same in both the Hopfield case and the Little case. 
The results are for N = 500. In particular, we find mf 0.90 at a = 0.25 for A > 1 
('^R < 0). The inset shows the variation of average recall time r and the fraction 
of limit cycles with A. Recall time r again shows a minimum at A = 1 (Ar = 0), 
as in the Glauber (Hopfield) case. An interesting observation is that at A = 1 the 
limit cycles of the Little model (parallel dynamics) disappear and all reach to fixed 
points with identical overlap properties as in the case of the Hopfield model. 
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Fig. 2. Average final overlap rrif (for the fixed point obtained) at a fixed loading capacity A = 0.25, 
for a network with N = 500 against Ar(=: (1 ~ A)/(l + A), positive A). The results are for 
parallel updating (Little model); overlap results for the Glauber dynamics (Hopfield model) are 
also indicated. The inset shows the variation of average recall time r and the fraction of limit 
cycles fi against Ar. 


For sequential updating, the variation of the average final overlap rrif with the 
loading capacity a, for some typical positive values of A, are shown in Fig. 3. The 
results are for N = 1000. The inset here shows variation of average recall time (r) 
with a for different A. For negative values of A, even for sequential updating, one gets 
limit cycles (e.g. the fraction of limit cycles 0.48, for A = -0.8 at a = 0.25) and 
the fixed point fraction decreases further for parallel (Little) dynamics. However, 
the overlap with the learned patterns decreases considerably and we do not see any 
significant memory capability (for storage in limit cycles) for negative values of A 
(typically nif ~ 0,42 for A = 0.8 at a = 0.25, using sequential updating). 

In order to investigate the finite size effect in the improved performance of the 
network for positive A values (for sequential dynamics), we studied the variation of 
average final overlap rrif with 1/N at fixed loading capacity a (= 0.20 and 0.25) 
and fixed A (= 1). These results, for 250 < iV < 4000, are shown in Fig. 4. The 
inset there shows the same variation of rrif with 1/N at fixed a (= 0.25) for A = 0. 
One can clearly compare the finite size effects observed: for A = 1, no significant 
variation of rrif with 1/N is observed for both a = 0.20 and 0.25 (mf c:; 0.94 for 
a = 0.20, and mf ~ 0.89 for a = 0.25 at A = 1 as AT — > oo). For A = 0, however 
there is significant finite size effect. From the extrapolation we find the extrapolated 
result for mf (for large N) to be around 0.29 here at this a (= 0.25); see the inset. 
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Fig. 3. Variation of average final overlap mf against loading capacity for some typical values of A, 
with sequential dynamics for N = 1000. The inset shows the variation of the recall time r against 
a for the same values of A. 



Fig. 4. Variation of average final overlap mf against 1/iV for fixed loading capacity a — 0.20 
and a = 0.25 for A = 1. The inset shows the same for a = 0.25 at A = 0. The results are for 
250 < TV < 4000. 
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It is to be noted that for positive A, the overlap and other properties saturate 
beyond A = 1 (Ar = 0) for any a. This is true even for very large values of A, 
in the updating (8). This is somewhat tricky in the sense that apparently one 
time step is almost skipped in every updating; practically {Si{t - 2)} (together 
with only a small fraction of Si{t - 1)) determines {5i(t)} through (8) and is kept 
in a ‘buffer’ for latter use and {Si{t - 1)} from the ‘buffer’ is used to determine 
{Si(t = 1)}. We have checked independently the performance of such an updating 
(5i(t+l) = sgn[I.jJijiX' Sj{t)+XSj{t-l))], instead of (8), with A' -> 0). We find for 
A' = 0, the Hopfield model result is recovered. In fact this can be seen from Fig. 2, 
where the overlap ruf decreases rapidly for very large values of A in (8) (A > 10). 
It may be noted, however, that the effective size of the domain of attraction of the 
fixed points depends very much on the ratio A/A'. This can be seen easily in the 
limi t A -4 0 but A' finite (or A' -+ 0 but A finite), where the change in energy of the 
network per updating step depends on the value of A' (or A ), which determines the 
effective shape of the energy landscape seen by the dynamics. 

As can be easily seen, for nonzero A, the updating dynamics in (8) does not 
minimize the ‘energy’ function E{t) defined as E{t) = — Ej>i JijSi{t)Sj{t), as 
the delay term contributes in the internal field. Some typical variations in E{t) 



Fig. 5. The variation of the configurational “energy” B{t) with time (updating iteration) t, for 
the sequential updating with A = 4 for two typical cases (arbitrarily distorted patterns) before 
and after reaching the fixed points. The inset shows the same for a distorted pattern with A = 0 
dynamics. 
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with iterations (updating) are shown in Fig. 5 for a = 0.5 and A = 4. (The large 
values of a and A are chosen here to ensure larger r, so that the effect can be 
displayed over longer time range). The results clearly show the contribution of the 
A term in escaping over the spurious valleys, in the “energy” landscape. The inset 
in Figure 5 shows the same variation for E(t) in the A = 0 (Hopfield) case. Here of 
course, it clearly decreases monotonically with time. 

We have also studied the performance of this dynamics for diluted network 
(with sequential dynamics; N = 500), where a random fraction c of the synaptic 
connections Jij removed. We show in Fig. 6, some results for overlap (rrif) against A 
at a fixed loading capacity a = 0.20 for some typical values of dilution concentration 
c. The inset shows the variation of mf against a at a fixed A (= 1) for some typical 
values of c. 
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Fig. 6. The variation of average final overlap mf with Ar (= (1 — A)/(l + A), positive A) at a fixed 
capacity a = 0.20 with sequential dynamics for N = 500 at some typical values of c(= 0.10, 0.20 
and 0.25). The inset shows the variation of mf with a at A = 1, for the same values of c. 


3. Approximate Analytical Solution in Case of Strong Dilution 

In this section we present an approximate analytical solution of the asymmetrically 
diluted version of the model described in the preceding section. In the limit of strong 
dilution, where only an infinitesimal fraction of the original number of synaptic 
connections remain, this model can be solved exactly (for A = 0 in (8)) by using the 
method developed by Derrida et al^ Following Derrida et al. we define the diluted 
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version of the Jij as 

= ( 9 ) 

jl 

where Cij and Cji are independent random variables and hence Jij is no longer 
symmetric. We have chosen the normalization l/K rather than 1/iVso as to give 
sensible results; as is the number of remaining connections. The probability 
distribution for Cij is assumed to have the form 

P(Qj) = §SiCij - 1) + (^1 - I j 6{Cij) . (10) 

With the above form of Jij the local field at site i at any time t is given by 

3 

= E + ^sj{t - 1)1 . (11) 

3 

For any state {Sj} of the network, we now break up the field hi{t) in Eq. (11) into 
a term coming from a particular pattern u and a remaining crosstalk (noise) term: 

= ^Ci E w + E - 1) + vr(t - 1) , ( 12 ) 

3 3 

where 

j 

and 

- 1) =1 E E - 1) . (13) 

ixi:u j 

We see that the dynamics at zero temperature is governed by the update rule (for 
asynchronous as well as synchronous dynamics) Si{t - 1) = sgn(/ii(t)), where hi{t) 
is given by Eq.(12). 

Like the fully connected network, in this case also we define the loading capacity 
a = p/iiT, where p is the number of pattern to be stored in the network. In the very 
diluted case we can calculate critical loading capacity ac (the maximum number 
of pattern we can store in the network with retrieval quality better than 90%). In 
order to calculate the critical loading capacity ac we examine the stability of the 
patterns To do that we find the overlap of an actual stable configuration 

Si = sgn(/ii) of the network with the pattern We define the overlap as 

]v ^ Cl S'i = — ^ Cij(jSj . 

^ 3 


m 


(14) 
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Both definitions of rnf' are equivalent (Hertz et al?). Using the dynamics 5i(t4*l) = 
sgii{hi{t)) with hi{t) given by Eq.(12), we have a self-consistent equation for 
given by 

i 

= sgnIm-'Ci) + - 1) + QivUt) + - 1))] . 

i (15) 

Assuming independent fluctuations at (t) and — the last term inside the third 
bracket of the above equation becomes a random variable whose mean is zero and 
whose variance is given by 

^2 _ + 2Xm^{t)m^{t — 1)] , (16) 

and has a Gaussian distribution given by 

-PC’?) = ^/^^ exp[-7?V2o-^] . (17) 


We now replace the summation in Eq.(15) by an integration over the random Gaus- 
sian variable rj and we have the evolution equation for the overlap 


i^{t + 1) = j drjP{r))sg!i[m^ (t) H- Xm^{t - 1) + 1 ?] 


m^{t) Xmy{t — 1) 


To replace the sum by the integral we had to assume the statistical independence of 
Q and t^^’s, so that the symmetry P{7]) = P(— ry) could be used. Moreover we are 
assuming that have the Gaussian distribution P{Ti) given by (17), even though 
that was strictly valid only for Si = . It can be shown that these assumptions are 

valid if (a) dilution is independent for ij and ji (which is indeed the case since Cij 
and Cji are assumed to be independent) and (b) we are in strongly diluted regime, 
K <. N (Hertz et al?). At the steady state (i.e. when the fixed point has been 
reached) we have m^{t) = m^{t — 1) = rrif and the self-consistent equation for the 
overlap becomes 

TUf = erf 

The critical loading capacity ac beyond which the only solution is mf = 0 and below 
which there are solutions with mf ^ 0, is given by^° 

, 2 (1 + A)" 

^ TT l-h 





( 20 ) 
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This increased loading capacity owing to the suppression of fluctuation due to 
the A term (see Eqs. (18) and (19)), can also be rationalized in a different way if 
one notes® that the ‘shadow’ neurons <Xi{t) (= 5i(t - 1)), representing the states of 
the neurons at the previous time steps, effectively increases the number of neurons 
from N to 2N (assuming there is no correlation between the states of the neurons 
and shadow neurons; away from fixed points). One can then rewrite the dynamical 
Eq. (8) in a single step updating form 



Such an asymmetric synaptic matrix might be argued (eg. Ref. 5) to be responsible 
for suppressing the spin glass noise, thereby increasing the loading capacity. In 
fact the largest eigenvalue of the synaptic matrix in (21) increases by a factor A (in 
the A 0 limit). This may be compared with that implied by (19), giving higher 
capacity due to suppression of noise. 

4. Discussions 

Our detailed numerical investigation clearly indicates that the recall time (r) and 
the overlap properties (measured through the overlap function mf) of a neural net- 
work with the synaptic connections formed using the Hebb’s rule (1), are clearly 
improved with appropriately weighted (Ac::!) time delayed dynamics (8) with se- 
quential and parallel updating. These indeed support the observations reported in 
Refs. 3 and 9. Our simulation results for rather large network sizes were used for 
systematic extrapolations (see Fig. 4), which suggests the improved performance 
of the network (for A > 1) to be valid in the thermodynamic limit. Approximate 
analytical consideration (see Sec. 3) also indicates the same. 

It may be noted that with any finite A, the dynamics (even for sequential up- 
dating) is not Glauber-like, since no energy function can be defined (as in the case 
of Hopfield model for A = 0) which gets minimized during the updating process. 
Also, it may further be noted that although the A term provides a momentum term, 
allowing the network to escape shallow spurious valleys (see the results for the time 
variation of the ‘energy’ function E{t) in Sec. 2; Figure 5), the fixed point equa- 
tions of the dynamics in (8) are the same as that of the Hopfield model: namely 
Si{t* -f 1) = sgn(/ii(r)) with /ii(r) = (1 -f- A)[Ej Jij5j(r)]. This indicates that 
the same Jij^ as defined in (1), contains more overlap states (with mf > 0.9 for 
a < 0.25 with A 1^’®) than those obtained using the Hopfield (Glauber) dy- 
namics (which gives only the fixed point states with rrif > 0.97 for a < 0.14 and 
mf = 0 for a > 0.14 A = 0). The use of the dynamics given by Eq. (8) (which for 
finite A does not necessarily minimize the energy in successive iterations) helps to 
access of these additional (over those for the Hopfield model) fixed points [already 
built-in with the synaptic connection given by Hebb’s rule (Eq. (1))]. In fact, as 
already shown,® see Sec, 3, the dynamics in (8) can be rewritten in a single time 
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step update form with redefined neuronal states, but then the S 3 Tiaptic interaction 
becomes asymmetric in an extended neuron state space (of size 2N), which may be 
argued for the better performance. The essential role of the A term (for positive 
A) to suppress the spin glass-like noise has been indicated in an approximate noise 
calculation in the extreme dilute limit (in Sec. 3). 

Acknowledgment 

PKM is grateful to DST, India, for financial assistance and to SEMP, Calcutta, for 
warm hospitality. BKC is grateful to A. C. C. Coolen and D. Sherrington for some 
useful comments and suggestions. We thank P. Sen for useful discussions. 

References 

1. D. J. Amit, H. Gutfreund and H. Sompolinsky, Ann. Phys. 173, 30 (1987). 

2. P. Peretto, An Introduction to the Modelling of Neural Networks, (Cambridge Univer- 
sity Press, Cambridge, 1992); J. Hertz, A. Krogh and R. G. Palmer, Introduction to 
the Theory of Neutral Computation, (Addison- Wesley, California, 1991); P. K. Maiti, 
in Modeling and Prediction in Complex System, ed. J. K. Bhattacharjee and A. K. 
Malhk (Norasa, Delhi, 1995) (in press). 

3. B. K. Chakrabarti and P. K. Dasgupta, Physica A186, 33 (1992). 

4. M. Ghosh, A. K. Sen, B. K. Chakrabarti and G. A. Kohring, J. Stat. Phys. 61, 501 
(1990); V. Banerjee and S. Puri, Int. J. Mod. Phys. B (1994) (in press). 

5. G. Parisi, J. Phys. A19, L675 (1986). 

6. B. Derrida, E. Gardener and A. Zippehus, Europhys. Lett. 4, 167 (1987). 

7. I, Kanter and H. Sompohnsky, Phy. Rev. A35, 380 (1987). 

8. C. M. Marcus and R. M. Westerfield, Phy. Rev. A42, 2410 (1990). 

9. P. Sen and B. K. Chakrabarti, Phys. Lett. 162, 327 (1992). 

10. V. Deshpande and C. Dasgupta, J. Phys. A24, 5015 (1991). 



